Displaying 1 result from an estimated 1 matches for "normalcovariance".
2007 Aug 13
1
simulate data from multivariate normal with pre-specified correlation matrix
...N<- 10 # number of samples
# define population correlation matrix sigma
sigma<-matrix(0,p,p) #creates a px p matrix of 0
rank<-2
for (i in 1:rank){
for (j in 1:rank){
rho<-0.75
sigma[i,j]<-rho^abs(i-j)
sigma[i,i]<-1
}
}
# create a data set of 10 observations from multivariate normalcovariance matrix sigma
library(MASS)
Xprime<-mvrnorm(N,rep(1,p), sigma )
Also, if I calculate
S<- cov(Xprime)
will S be the sample correlation matrix?
Thank you very much for your time!