Displaying 1 result from an estimated 1 matches for "mix2".
Did you mean:
mix
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...