信頼区間

Question

私のデータが2列で構成されています。以下のような時間と累積数：

time <- c(1:14) 
cum.num <- c(20, 45, 99, 195, 301, 407, 501, 582, 679, 753, 790, 861, 1011, 1441) 

は私の非線形関数である：

c1*cos(0.6731984259*time)+c2*sin(0.6731984259*time)+c3*(time)^2+c4*time+c5 

私の目標は、これをモデル化することです関数を非線形回帰を使用してRでnls()とし、信頼区間を計算します。私は、次のdonr持って、次のエラーを得た

confint(m1.fit) 

私は、このコマンドを発行すると：：

m1.fit<-nls(cum.vul~c1*cos(0.6731984259*time)+c2*sin(0.6731984259*time)+c3*(time)^2+c4*time+c5,start=list(c1=-50,c2=-60,c3=5,c4=8,c5=100)) 

信頼区間を計算しながら、私はエラーを得た、私は次のことを試してみました

Waiting for profiling to be done... 
Error in prof$getProfile() : 
step factor 0.000488281 reduced below 'minFactor' of 0.000976562 

誰もこの点で私を助けることができますか？

Answer 1

てみパッケージnlstools：

> nlstools::confint2(m1.fit) 
     2.5 %  97.5 % 
c1 -48.556270 54.959689 
c2 -175.654079 -45.216965 
c3 3.285062 9.529072 
c4 -49.254627 46.007629 
c5 -34.135835 272.864743`