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 arra
y

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.