Возникли проблемы при получении coefplot () для построения двух регрессий друг над другом

Question

Я пытаюсь использовать coefplot (), как описано в этой статье: http://www.r -bloggers.com / визуализация-коэффициенты-регрессии-в-r /

Однако, когда я запускаю этот точный код, я получаю только одну регрессионную диаграмму, а не 3. Вот скриншот, показывающий точный код, который я запустил, плюс график вывода.http://i.imgur.com/YtDND.png

Я действительно не уверен, что еще делать.Ваша помощь будет высоко ценится.

Answer 1

Я укажу, что coefplot2 (мое расширение coefplot, которое имеет некоторые приятные функции, но все еще требует дополнительной работы), доступное от R-forge, делает это:

library(coefplot2)
coefplot2(list(M2,M1,M3),col=c("black","red","blue"),legend=TRUE)

Answer 2

Код для coefplot больше не принимает аргумент для смещения. Его нет в документации и нет в списке официальных документов. Вы можете сделать версию, которая делает это, изменив код для coefplot:

Тип coefplot2 return . Скопируйте и вставьте функцию в командную строку и поставьте перед ней

coefplot2 <-  # the rest of the pasted function should follow

затем добавьте voffset=0 к списку формальностей и измените эту строку:

arrows(ci1, (1:k), ci2, (1:k), lty = lty[1], lwd = lwd[1], col = col,

К этому

arrows(ci1, (1:k)+voffset, ci2, (1:k)+voffset, lty = lty[1], lwd = lwd[1], col = col,

И измените линию точек на:

points(cf , (1:k)+voffset, pch = pch, col = col)

Затем нажмите Enter, и у вас должна появиться новая функция coefplot2. Тогда это должно работать

coefplot2(M2, xlim=c(-2, 6) )
par(new=TRUE)    # could not get the add=TRUE argument to work either.
coefplot2(M3, col="blue", xlim=c(-2, 6), voffset=0.4)
par(new=TRUE)
coefplot2(M1, col="red", xlim=c(-2, 6) , voffset=0.2)