Displaying 1 result from an estimated 1 matches for "ydatasd".
Did you mean:
datas
2009 Jun 17
0
nls with weights
...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)
yDataSd <- apply(yMat, 2, sd)
plot(yData ~ xData, ylim = c(min(yData - yDataSd), max(yData + yDataSd)))
arrows(xData, yData - yDataSd, xData, yData + yDataSd, angle = 90,
length = 0.1, code = 3)
# without weights
model <- nls(yData ~ sigmoid(xData,bottom, top, slope, logec50), start =
list(bottom...