pdf("rpd7_9r.pdf")


p=matrix(c(2320,3300,3290,3460,3770,4210,5470,5510,6120,6780),ncol=1)
y=matrix(c(2380,3190,3270,3530,3980,4390,5400,5770,6890,8320),ncol=1)


delta=p-y

nobs=nrow(y)

msep1=(t(delta) %*% delta)/nobs;
delta
msep1



msep2=((nobs-1)*var(delta)/nobs)+(mean(delta)^2)
mean(delta) 
var(delta) 
msep2

dev.off()