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))
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