Zadania 2, Studia, Stopień 2 Semestr I, Statystyczne systemy informatyczne w AD, R moje prace, Wykłady i zadania, 2. ...
[ Pobierz całość w formacie PDF ]
#ZAD1wzrost=read.csv2("wzrost.csv", header=F)summary(wzrost)quantile(wzrost[,1], 0.25)quantile(wzrost[,1], 0.75)median(wzrost[,1])sd(wzrost[,1])(r=max(wzrost)-min(wzrost))library(e1071)#kurtozakurtosis(wzrost[,1],type=2) #nieobiciazony bo type2 i e1071k=n*(n+1)/(n-1)/(n-2)/(n-3)*sum((wzrost[,1]-mean(wzrost[,1]))^4)/(sd(wzrost[,1]))^4-3*(n-1)^2/(n-2)/(n-3))#Skosnoscn=length(wzrost[,1])s=n/(n-1)/(n-2)*sum((wzrost[,1]-mean(wzrost[,1]))^3)/(sd(wzrost[,1]))^3s=e1071::skewness(wzrost[,1],type=2) #trzeba przez e1071 obliczyc#przedział ufności dla wartosci oczekiwaneja=0.1qt(a,n-1,lower.tail=T)(m=mean(wzrost[,1])-qt(1-a/2,n-1,lower.tail=T)*sd(wzrost[,1])/(n)^(1/2))(n=mean(wzrost[,1])+qt(1-a/2,n-1,lower.tail=T)*sd(wzrost[,1])/(n)^(1/2))cat("przedział ufnośći dla wartosci oczekiwanej wynosi:","[",m,",",n, "]")#Przedział ufnosci dla wariancjib=n*(sd(wzrost[,1])^2)/(qchisq(1-a/2,n-1,lower.tail=T))c=n*(sd(wzrost[,1])^2)/(qchisq(a/2,n-1,lower.tail=T))cat("przedział ufnośći dla wariancji wynosi:", "[",b,",",c, "]")#ZAD2zakupy=read.csv2("zakupy.csv")N=length(zakupy[,3])zakupy[,2]w=(zakupy[,2]=="K")zakupy$WYDATEK[w]kobiety=zakupy$WYDATEK[zakupy[,2]=="K"]man=zakupy$WYDATEK[zakupy[,2]=="M"]mean(kobiety)mean(man)var(kobiety)var(man)median(kobiety)median(man)quantile(kobiety,0.25)quantile(kobiety,0.75)quantile(man,0.25)quantile(man,0.75)skewness(kobiety,type=2)skewness(man,type=2)kurtosis(kobiety,type=2)kurtosis(man,type=2)#ZAD3X=c(1,2,3,4,5)P=c(7/30, 1/6, 4/15, 2/15, 1/5)wartocz=function(P,X){EX=0for (i in 1:5) {EX=EX+X[i]*P[i]}return(EX)}wartocz(P,X)EXX=function(P,X) {EXX=0for(i in 1:5){EXX=EXX+X[i]^2*P[i]}return(EXX)}EXX(P,X)Wariancja=function(P,X) {VarX=EXX(P,X)-wartocz(P,X)^2return(VarX) }Wariancja(P,X)#ZAD4pu_test=function(n,mi,sigma,poz){w=c(rnorm(n,mean=mi,sd=sigma))m=mean(w)-qt(1-(1-poz)/2,n-1,lower.tail=T)*((sigma/n)^(1/2))n=mean(w)+qt(1-(1-poz)/2,n-1,lower.tail=T)*((sigma/n)^(1/2))if(mi>=m && mi<=n)y=TRUEelsey=FALSEy}pu_test(100,0,1,0.9)ileliczb=function(n,mi,sigma,poz) {T=0F=0for(i in 1:1000) {if(pu_test(n,mi,sigma,poz)==TRUE)T=T+1elseF=F+1}cat(" TRUE otrzymano:",T,"razy\n","FALSE otrzymano:",F,"razy\n")}ileliczb(100,0,1,0.9) [ Pobierz całość w formacie PDF ]