SMSfacenorm (R 2.13.1)

SMSfacenorm computes Flury faces for two normal samples: x~N1(0,1) and y~N2(2,1)

Download File

Tomas Marada

Multivariate Statistics: Exercises and Solutions

Tue, July 10 2012 by Dedy Dwi Prastyo

-


Description: -

Description: -


rm(list=ls(all=TRUE))
graphics.off()
install.packages("aplpack") 
library(aplpack)
q = matrix(0, ncol = 36);
q[,13]=1;                   # right face line
q[,14]=1;                   # right darkness of hair
q[,31]=1;                   # left face line
q[,32]=1;                   # left darkness of hair
x = rnorm(50);              # generates 50 standard normal distributed data
y = rnorm(50) + 2;          # generates 50 normal distributed data with mean 2
z = t(cbind(t(x),t(y)));    # puts the data together in matrix
par(ask=TRUE)
faces(z, q, face.type=1, main="random faces", col.nose=rainbow(2),col.eyes=rainbow(2, start=0.6,end=0.85),col.hair=terrain.colors(25),col.lips=rainbow(20,start=0,end=0.2),col.face=heat.colors(15),col.ears=rainbow(20,start=0,end=0.2), nrow.plot=10,ncol.plot=5,plot.faces=T);   #computes the Flury faces