Displaying 2 results from an estimated 2 matches for "lmix2".
Did you mean:
fmix2
2004 Aug 03
1
Using MASSv3's example from 8.7 in R?
...================
Example from MASSv3 8.7
=========================================================================
library(MASS)
data(geyser)
attach(geyser)
truehist(waiting, xlim=c(35,110), ymax=0.04, h=5)
width.SJ(waiting)
wait.dns <- density(waiting, 200, width=10.24)
lines(wait.dns, lty=2)
lmix2 <- deriv3(
~ -log(p*dnorm((x-u1)/s1)/s1 + (1-p)*dnorm((x-u2)/s2)/s2),
c("p", "u1", "s1", "u2", "s2"),
function(x, p, u1, s1, u2, s2) NULL)
p0 <- c(p=mean(waiting < 70), u1=50, s1=5, u2=80, s2=5)
tr.ms <- function(info, thet...
2001 May 06
1
deriv3 example on Venables/Ripley page 263
What does it require to run the example on page 263
of Venables/Ripley book using R 1.3.0?
I got the following error and I have no clue how to correct it.
Thank you very much. (If you reply to the list, please send a
copy to me.)
> lmix2 <- deriv3(
+ ~ -log(p*dnorm((x-u1)/s1)/s1 + (1-p)*dnorm((x-u2)/s2)/s2),
+ c("p", "u1", "s1", "u2", "s2"),
+ function(x, p, u1, s1, u2, s2) NULL)
Error in deriv3.formula(~-log(p * dnorm((x - u1)/s1)/s1 + (1 - p) * dnorm((x - :
Function...