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