se vc se refere a "eu tenho um *vetor* de dados que vem da mistura de 2 normais desconhecidas e quero identicar: A) os parametros de cada normal; B) qual observacao vem de qual distribuicao", entao vc esta' tratando de um problema de mistura (normal) de duas componentes... e eu faria: set.seed(10) dist1<-rnorm(n=400,mean=10,sd=5) dist2<-rnorm(n=400,mean=35,sd=8) dist.c <- c(dist1, dist2) library(mixtools) res <- normalmixEM(dist.c) res$mu ## compare com as medias q vc usou na geracao res$sigma ## compare com os desvios-padrao usados na geracao apply(res$posterior, 1, which.max) ## a q distribucao cada obs pertence b 2013/6/3 ASANTOS <alexandresantosbr@yahoo.com.br>:
Boa tarde pessoal,
Simulei duas distribuições normais e representei em um histograma que obviamente assumiu a forma bimodal, sendo: ### set.seed(10) dist1<-rnorm(n=400,mean=10,sd=5) dist2<-rnorm(n=400,mean=35,sd=8) dist.c<-rbind(dist1,dist2) # hist(dist.c) # Gostaria de saber se existe algum procedimento que me permita separar dist1 e dist2 com seus respectivos intervalos de confiança, imaginando que eu não conheça os parâmetros media e desvio padrão que os criaram?
Obrigado
-- ====================================================================== Alexandre dos Santos Proteção Florestal Coordenador do curso Técnico em Florestas Vice Coordenador do curso de Engenharia Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos@cas.ifmt.edu.br ======================================================================
_______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.