Аблайн для ГЛС не будет сюжет - PullRequest
1 голос
/ 16 апреля 2020

Я пытаюсь построить модель ниже с линией регрессии, и я могу сгенерировать график, но ни одна из команд аблайна, которые я пробовал (также ниже), кажется, не работает. R будет казаться счастливым и не будет выводить никаких сообщений об ошибках с помощью первой команды abline, которую я попробую, но на графике не появится ни одной строки. Я попытался добавить функцию abline в сам код, но это тоже не сработало. Любая помощь будет оценена.

model2<-gls(transformedlongevity~transformedproportion, data=independencedata, 
            correlation=corPagel(1,trimtree))
summary(model2)
plot(model2, xlab= "Transformed Value of Proportion of Life Spent Under Parental Care (Years)", 
     ylab="Transformed Value of Maximum Life Span (Years)", 
     main= "Time Spent Under Parental Care Relative 
     to Longevity of Passerine Species")


abline(gls(transformedproportion~transformedlongevity, independencedata), untf=T)
abline(a=0, b=coef(model2)) 
abline(coef(model2)) 

1 Ответ

0 голосов
/ 16 апреля 2020

Пожалуйста, включите пакет, использованный в вашем сообщении. Так что если ваш gls взят из nlme, то график, который вы видите, предназначен для остатков, а не для прогнозируемых значений. Используя пример, заданный в nlme:

library(nlme)
fm1 <- gls(follicles ~ Time, Ovary,
            correlation = corAR1(form = ~ 1 | Mare))
plot(fm1)

enter image description here

Таким образом, подобранная линия здесь не имеет смысла. Скорее всего, вы хотите сделать:

plot(Ovary$Time,Ovary$follicles)
abline(fm1,col="blue")
Добро пожаловать на сайт PullRequest, где вы можете задавать вопросы и получать ответы от других членов сообщества.
...