[R] Functions in formulae ??
Jeremy Clark
jeremyclarkmel at gmail.com
Thu Nov 21 11:57:29 CET 2013
Dear All,
In the following simple case I can't seem to get an improved fit,
despite trying all of the control possibilities. As there seem to be
no examples anywhere which show use of functions such as "dnorm"
within a formula, and as I am not confident at all that my formula is
correctly configured, I would appreciate some feedback ! Many thanks
in advance.
library(minpack.lm)
gg<-c(170, 171, 172, 173, 174, 175,
176, 177, 178, 179, 180, 181, 182,
183, 184, 185, 186, 187, 188, 189,
190, 191, 192, 193, 194, 195, 196,
197, 198, 199, 200, 201, 202,
203, 204, 205, 206, 207, 208, 209)
dd<-c(1673, 1659, 1738, 1687, 1882, 2010, 2202,
2248, 2409, 2417, 2215, 2279, 2539, 2479, 2341,
2395, 2314, 2404, 2369, 2254, 2048, 1899,
1892, 1744, 1333, 1183, 982, 772, 526,
451, 335, 253, 157, 108, 67, 44,
26, 12, 10, 0)
dd100000 <- dd * 0.00001
dd100000gg <- data.frame(dd100000, gg)
## modLM4:
modLM4 <- nlsLM(dd100000 ~ dnorm(gg, mu,sigma, log = FALSE), data =
dd100000gg, start = c(mu = 190, sigma = 10), trace = TRUE)
windows()
## plot data
plot(gg, dd100000, main = "modLM4")
## plot fitted values
lines(gg, fitted(modLM4), col = 2, lwd = 2)
More information about the R-help
mailing list