Displaying 1 result from an estimated 1 matches for "ydata3".
Did you mean:
data3
2009 Jun 17
0
nls with weights
...t;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)
yDataSd <- apply(yMat, 2, sd)
plot(yData ~ xData, ylim = c(min(yData - yDataSd), max(yData + yDataSd)))
arrows(xData, yData - yDataSd, xData, yData + yD...