預測更大範圍

Question

我想擬合更寬的x範圍的數據並預測y值。

讓我們假設我有'IRIS的數據集，並從該post

library(dplyr) 
cc <- iris %>% 
    group_by(Species) %>% 
    do({ 
    mod <- nlsLM(Sepal.Length ~ k*Sepal.Width/2+U, start=c(k=10,U=5), data = ., trace=F, control = nls.lm.control(maxiter=100)) 
    pred <- predict(mod, newdata =.["Sepal.Width"]) 
    data.frame(., pred) 
    }) 

使用下面的數據進行預測這是擬合圖

我想，以配合這個數據更寬Sepal width範圍使得

new.range<- data.frame(x=seq(2,10,length.out=20)) 

和修改腳本

pred <- predict(mod, newdata =new.range) 

繪製的（功能new.range擬合

library(ggplot2) 

ggplot(cc,aes(y=Sepal.Length,x=Sepal.Width ,col=factor(Species)))+ 
    geom_point()+ 
    facet_wrap(~Species)+ 
    geom_line(aes(x=new.range,y=pred),size=1) 

錯誤（...，row.names = NULL，check.rows = FALSE， check.names = TRUE，：參數意味着不同的行數： 20，150

我不明白爲什麼要這樣做錯誤。我想pred從new.range計算，所以他們應該有相同的長度？

類似的帖子

using-predict-in-nls

trouble-with-predict-function-in-r

predict-maybe-im-not-understanding-it?

Answer 1

這是東西，達到你想要的。原始問題的原因是，在您的迴歸中，預測的名稱是Sepal.width而不是x，而您的預測完全不使用您的new.range，因此您必須執行類似new.range<- data.frame(Sepal.Width=seq(2,10,length.out=50))的操作來對您的new.range進行預測。

另一個問題是，你必須使new.range的長度爲50，以便pred和new.range適合原始數據幀。

然後你可以畫出你想要的圖，注意new.range變成Sepal.Width.1。

library(dplyr) 
cc <- iris %>% 
    group_by(Species) %>% 
    do({ 
     mod <- nlsLM(Sepal.Length ~ k*Sepal.Width/2+U, start=c(k=10,U=5), data = ., trace=F, control = nls.lm.control(maxiter=100)) 
     new.range<- data.frame(Sepal.Width=seq(2,10,length.out=50)) 
     pred <- predict(mod, newdata =new.range) 
     # pred <- predict(mod, newdata =.["Sepal.Width"]) 
     data.frame(., new.range, pred) 

    }) 

library(ggplot2) 

ggplot(cc,aes(y=Sepal.Length,x=Sepal.Width ,col=factor(Species)))+ 
    geom_point()+ 
    facet_wrap(~Species)+ 
    geom_line(aes(x=Sepal.Width.1,y=pred),size=1)