Displaying 1 result from an estimated 1 matches for "ydata1".
Did you mean:
data1
2009 Jun 17
0
nls with weights
...n(yData)??? Can anybody explain, why it ends up with:
"singular gradient"
8< ------------------------------------------------------- >8
sigmoid <- function(x, bottom, top, slope, logec50) { bottom + (top -
bottom) / (1 + 10^((logec50-x)*slope)) }
xData <- seq(-10,10,1.2)
yData1 <- sigmoid(xData,10,1,-0.2,3) + rnorm(1:length(xData), sd = 0.5)
yData2 <- sigmoid(xData,10,1,-0.2,3) + rnorm(1:length(xData), sd = 0.2)
yData3 <- sigmoid(xData,10,1,-0.2,3) + rnorm(1:length(xData), sd = 1.1)
yMat <- rbind(yData1, yData2, yData3)
yData <- apply(yMat, 2, mean)
yDataS...