search for: ydata3

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...