4パラメーターロジスティックモデルまたは3パラメーターロジスティックモデルのいずれかに合うようにtry（またはtryCatch）を使用するには

Question

パッケージdrcを使用して、毒性データに4パラメーターロジスティック回帰モデルを適合させたいと思います。この警告に

exp.df <- as.data.frame(matrix(nrow = 14, ncol = 2)) 

exp.df[,1] <- c(3200, 9600, 25600, 25600, 25600, 1600, 800, 
       6400, 19200, 0, 12800, 1200, 400, 2400) 

exp.df[,2] <- c(0.5855615, 0.9625668, 0.4064171, 0.4973262, 
       0.4732620, 1.0000000, 0.6764706, 0.4652406, 
       0.5106952, 0.7566845, 0.5427807, 0.5106952, 
       0.5935829, 0.4759358) 

names(exp.df) <- c("Dose", "Response") 

Dose_Response <- drm(Response ~ Dose, data = exp.df, 
        fct = LL.4(), type = "binomial") 

結果：OPTIMで

エラー（startVec、opfct、ヘッセ= TRUE、方法= optMethod、 コントロール=リスト（maxit回= maxit回、これは例示的なデータセットであります非有限有限差分値 [4] が制約drmOpt（opfct、opdfct1、startVecSc、optMethod、warnValにエラーが、私は、3パラメータのロジスティックのREGを実行した場合の収束が

に失敗レジーソン、それは動作します。

Dose_Response <- drm(Response ~ Dose, data = exp.df, 
        fct = LL.3(), type = "binomial") 

私は、まず3パラメータロジスティック回帰を4パラメータロジスティックをテストしようかtryCatchを使用することはできますか？それとも別の解決策がありますか？

Answer 1

あなたのcluster(id)とは何ですか？私は偽の1を追加した後、あなたのデータはLL.3()とLL.4()の両方で解決されました。