R kode som illustrerer binormalfordeling

 

 

# Den binormale fordelingen er beskrevet p? sidene 254-257 i l?reboka.

 

# Vi vil generere n=100 observasjoner fra den binormale fordelingen med E(X)=E(Y)=0 og V(X)=V(Y)=1 for ulike verdier av korrelasjonen (rho) og plotte observasjonene

library(MASS)

n=100

par(mfrow=c(2,2))

r=0.0

plot(mvrnorm(n,c(0,0),matrix(c(1,r,r,1),2,2)),xlim=c(-3,3),ylim=c(-3,3),xlab="x",ylab="y",main="rho=0")

r=0.30

plot(mvrnorm(n,c(0,0),matrix(c(1,r,r,1),2,2)),xlim=c(-3,3),ylim=c(-3,3),xlab="x",ylab="y",main="rho=0.30")

r=0.60

plot(mvrnorm(n,c(0,0),matrix(c(1,r,r,1),2,2)),xlim=c(-3,3),ylim=c(-3,3),xlab="x",ylab="y",main="rho=0.60")

r=0.90

plot(mvrnorm(n,c(0,0),matrix(c(1,r,r,1),2,2)),xlim=c(-3,3),ylim=c(-3,3),xlab="x",ylab="y",main="rho=0.90")

par(mfrow=c(1,1))

 

# R kommandoene ovenfor er mer kompliserte enn de komandoene vi har brukt f?r, s? du trenger ikke ? forst? detaljene her.