In [1]:
# rm(list=ls()) 
options(OutDec = ",")
#==============================================================================
# Descricao: representacao grafica de dados multivariados no R
#
# Note que o R entende variaveis por colunas e as amostras por linha
# Entao, e' necessario transpor a matriz de dados
#
#==============================================================================
set.seed(5485)                    # fixando a semente pseudo-aleatorio
p       <- 3                      # numero de variaveis
n       <- 10000                  # tamanho da amostra
X       <- matrix(rnorm(p*n),p,n) # definindo X

par(mfrow=c(1,1),lwd=2.0,cex.lab=1.5,cex.axis=1.5,
     lab=c(10,5,5),mar=c(0,1,0,2.5),xpd=T,cex.main=2.0)
pairs(t(X),pch=15,xlim=c(-4.5,4.5),ylim=c(-4.5,4.5))
No description has been provided for this image
In [2]:
print(apply(X,1,"mean"))  # vetor (linhas) de medias amostral
print(var(t(X)))          # matriz de covariancias amostral
print(cor(t(X)))          # matriz de correlacoes amostral
[1]  0,005208232 -0,004029804  0,003825141
           [,1]        [,2]        [,3]
[1,] 1,01385413  0,01349226  0,02016458
[2,] 0,01349226  1,00964929 -0,01230010
[3,] 0,02016458 -0,01230010  1,00726439
           [,1]        [,2]        [,3]
[1,] 1,00000000  0,01333557  0,01995398
[2,] 0,01333557  1,00000000 -0,01219696
[3,] 0,01995398 -0,01219696  1,00000000