Modelo linear de segundo grau
Sei que posso rodar dois modelos e fazer o LRT, entre outras análises, mas por algum motivo preciso da rsm! Alguém saberia o motivo deste modelo não "rodar"? x=c(1,1,1,2,2,2,3,3,3,4,4,4,5,5,5) y=c(4,8,5,9,12,9,10,13,12,12,14,13,11,11,12) plot(y~x) require(rsm) m3=rsm(y~ FO(x)+FO(x^2)) Error in dimnames(LM$B) <- list(fonm, fonm) : comprimento de 'dimnames' [1] não é igual ao tamanho do array Obg
André, O seu modelo não deveria ser especificado como: *> m3=rsm(y~ FO(x)+PQ(x))*
m3
Call: rsm(formula = y ~ FO(x) + PQ(x)) Coefficients: (Intercept) FO(x) PQ(x) -0.1333 6.7190 -0.8810 HTH -- Cesar Rabak 2018-05-21 14:45 GMT-03:00 Andre Oliveira via R-br <r-br@listas.c3sl.ufpr.br
:
Sei que posso rodar dois modelos e fazer o LRT, entre outras análises, mas por algum motivo preciso da rsm!
Alguém saberia o motivo deste modelo não "rodar"?
x=c(1,1,1,2,2,2,3,3,3,4,4,4,5,5,5) y=c(4,8,5,9,12,9,10,13,12,12,14,13,11,11,12) plot(y~x) require(rsm) m3=rsm(y~ FO(x)+FO(x^2))
*Error in dimnames(LM$B) <- list(fonm, fonm) : comprimento de 'dimnames' [1] não é igual ao tamanho do array*
Obg
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.
Verdade César, obrigado! André Oliveira Souza. Em quinta-feira, 24 de maio de 2018 08:58:04 BRT, Cesar Rabak <cesar.rabak@gmail.com> escreveu: André, O seu modelo não deveria ser especificado como:
m3=rsm(y~ FO(x)+PQ(x))> m3 Call:rsm(formula = y ~ FO(x) + PQ(x)) Coefficients:(Intercept) FO(x) PQ(x) -0.1333 6.7190 -0.8810 HTH--Cesar Rabak
2018-05-21 14:45 GMT-03:00 Andre Oliveira via R-br <r-br@listas.c3sl.ufpr.br>: Sei que posso rodar dois modelos e fazer o LRT, entre outras análises, mas por algum motivo preciso da rsm! Alguém saberia o motivo deste modelo não "rodar"? x=c(1,1,1,2,2,2,3,3,3,4,4,4,5, 5,5) y=c(4,8,5,9,12,9,10,13,12,12, 14,13,11,11,12) plot(y~x) require(rsm) m3=rsm(y~ FO(x)+FO(x^2)) Error in dimnames(LM$B) <- list(fonm, fonm) : comprimento de 'dimnames' [1] não é igual ao tamanho do array Obg ______________________________ _________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/ cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br- guia) e forneça código mínimo reproduzível.
participantes (2)
-
Andre Oliveira -
Cesar Rabak