Displaying 4 results from an estimated 4 matches for "verosimiglianza".
2004 Sep 20
3
montecarlo simulation
Hy!
I would like to know how run a montecarlo simulation with R.
Thank you!!!!
Francesca Matalucci
__________________________________________________________________
Accesso Internet Gratis per utenti Excite! Attivalo subito!
http://www.excite.it/hitech/accesso
Il Mio Excite. Personalizza la tua Home page Excite come vuoi tu!
http://www.excite.it
AAA/Relazioni. Sfoglia gli annunci e trova la
2008 Jun 11
1
Some problem with the OLS
? stato filtrato un testo allegato il cui set di caratteri non era
indicato...
Nome: non disponibile
URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20080611/2da35c97/attachment.pl>
2009 May 22
0
EM algorithm mixture of multivariate
...rm(y, mu01, sd01, log=FALSE)
f2<-dmvnorm(y, mu02, sd02, log=FALSE)
p=alpha0*f1+(1-alpha0)*f2
scatterplot3d(y[,1], y[,2], p , highlight.3d=TRUE, col.axis="blue",
col.grid="lightblue",main="val iniz mistura normali multivariata",
angle=120, pch=20)
#verosimiglianza iniziale
l0<-sum(log(p))
l1<-l0
alpha<-alpha0
mu1<-mu01
mu2<-mu02
sd1<-sd01
sd2<-sd02
for (iter in 1:itermax)
{
#passo E
for (i in 1:n) {
tau[i,1]<-(alpha*f1[i])/p[i]
tau[i,2]<-((1-alpha)*f2[i])/p[i]
}
#passo M
alpha= mean(tau[,1])
mu1=colSums(tau[,1]*y)/sum(tau[,...
2009 May 22
0
EM algorithm mixture of multivariate gaussian
...rm(y, mu01, sd01, log=FALSE)
f2<-dmvnorm(y, mu02, sd02, log=FALSE)
p=alpha0*f1+(1-alpha0)*f2
scatterplot3d(y[,1], y[,2], p , highlight.3d=TRUE, col.axis="blue",
col.grid="lightblue",main="val iniz mistura normali multivariata",
angle=120, pch=20)
#verosimiglianza iniziale
l0<-sum(log(p))
l1<-l0
alpha<-alpha0
mu1<-mu01
mu2<-mu02
sd1<-sd01
sd2<-sd02
for (iter in 1:itermax)
{
#passo E
for (i in 1:n) {
tau[i,1]<-(alpha*f1[i])/p[i]
tau[i,2]<-((1-alpha)*f2[i])/p[i]
}
#passo M
alpha= mean(tau[,1])
mu1=colSums(tau[,1]*y)/sum(tau[,...