[R] Mixture normal distribution

aleid2001@yahoo.com aleid2001 at yahoo.com
Wed Feb 8 13:18:23 CET 2006


Dear R helper,


I hope that u can help me to sort out my problem
because I sent an E-mail last night to R-list but I
have not receive any help and at the same time I think
this problem is not so hard.

I have used the following functions before 

> K<-10
> prime<-c(2,3,5,7,11,13,17)
> UN<-seq(1:K)%*%t(sqrt(prime))
> U1<-UN-as.integer(UN)
> U<-matrix(qnorm(U1),K,7)
> U

            [,1]        [,2]        [,3]        [,4]  
     [,5]       [,6]
 [1,] -0.2167193  0.61902728 -0.71900806  0.37387492
-0.47715805  0.2677426
 [2,]  0.9479680 -0.09010569 -0.06990169 -0.54900030 
0.34047235 -0.8026015
 [3,] -0.6978334 -0.85544499  0.54814530  1.53212280 
1.64363675  0.9026864
 [4,]  0.4038929  1.46253909  1.59168198  0.20958766
-0.62343558 -0.1962555
 [5,] -1.4678844  0.41315654 -0.91407055 -0.74294809 
0.20989183 -1.9148412
 [6,] -0.0369025 -0.27331681 -0.21109168  1.14796196 
1.28012119  0.3406266
 [7,]  1.2786790 -1.15348516  0.39201340  0.05080407
-0.78449952 -0.7099779
 [8,] -0.4853656  1.06431267  1.21882042 -0.97005119 
0.08280907  1.0127500
 [9,]  0.6065405  0.22357828 -1.15223725  0.88440727 
1.03481830 -0.1257587
[10,] -1.0707736 -0.46627852 -0.35664238 -0.10670096
-0.96909898 -1.5936006
             [,7]
 [1,] -1.15960110
 [2,] -0.68646092
 [3,] -0.33366316
 [4,] -0.01899511
 [5,]  0.29375681
 [6,]  0.63913903
 [7,]  1.08816734
 [8,]  2.16601566
 [9,] -1.23750068
[10,] -0.73537275


Now I want instead of using (qnorm) the inverse of
slandered normal distribution, I want to use the
mixture of normal distribution, and also, I want to
know how to read the data to txt file.


Many thanks for your help in advance.

Al-Eid



		
___________________________________________________________ 
Win a BlackBerry device from O2 with Yahoo!. Enter now. http://www.yahoo.co.uk/blackberry




More information about the R-help mailing list