Displaying 1 result from an estimated 1 matches for "ydata2".
Did you mean:
data2
2009 Jun 17
0
nls with weights
...ent"
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)
yDataSd <- apply(yMat, 2, sd)
plot(yData ~ xData, ylim = c(min(yData - yDataS...