2017-02-26 89 views
0

我正在尝试使用mixtox包执行Beckon曲线拟合。为了获得起始值,我使用了tuneFit()函数。但是,我正在收到下面提到的错误消息。我也试着直接用我能想到的初始值进行拟合,但也收到错误信息。任何想法如何解决这个问题?Mixtox tuneFit函数获取Beckon曲线拟合的初始估计值:参数的长度为零

library(mixtox) 
via_data$Concentration <- c(0.1974457, 0.3980034, 0.7921297, 1.5754719, 
          3.1546263, 6.2702489, 12.6250182, 25.2772970, 
          25.1103882, 49.9450449, 74.6800090) 
via_data$Viability <- c(100.611247, 94.621027, 96.210269, 87.530562, 
         80.562347, 62.224939, 39.119804, 
         30.806846, 30.684597, 22.860636, 2.567237) 

via_tuneFit <- tuneFit(via_data$Concentration, via_data$Viability, 
       eq="Beckon") 
Error in if (nrow(param) > highBar) param <- param[sample(nrow(param), : 
    argument is of length zero 

via_curvFit <- curveFit(x=via_data$Concentration, rspn=via_data$Viability, 
       eq="Beckon", rtype = "continuous", 
       param = c(77, 26, 1, 1, -6)) 

Processing continuous quantal dose-response data 
Error in solve.default(t(jac) %*% jac) : 
    system is computationally singular: reciprocal condition number = 2.79026e-31 

任何帮助将不胜感激。

感谢, Krina

回答

0

这个工作。

via_curvFit <- curveFit(x=via_data$Concentration, rspn=via_data$Viability, 
       eq="Beckon", rtype = "continuous", param = c(1.32, 1.26, -0.94, 1.7, 0.68))