Programy R

######### Estymator regresyjny (LPBZ)

estreg_sr<-function(dane,N,xsrpop) {
#dane - pierwsza kolumna x, druga kolumna y
n<-nrow(dane)
xsrs<-mean(dane[,1])
ysrs<-mean(dane[,2])
eB<-cov(dane[,1],dane[,2])/cov(dane[,1],dane[,1])
r<-cor(dane[,1],dane[,2])
est<-ysrs-(xsrs-xsrpop)*eB
eD2<-(N-n)*var(dane[,2])*(1-r^2)/(N*n)
eD<-sqrt(eD2) #wartość estymatora średniego błędu szacunku
egam<-100*eD/est #wartość estymatora średniego względnego błędu szacunku
wynik<-matrix(c(est,eD2,eD,egam),ncol=1)
########### poniższe 3 linie formatujące wynik końcowy mogą być pominięte na kolokwium
wynik<-formatC(wynik,format="f")
nazwa<-list(c("estymator","estymator wariancji","estymator średniego błędu szacunku","estymator średniego względnego błędu szacunku (w %)"),c("wartość"))
dimnames(wynik)<-nazwa
###########
wynik
}