Advertisement
Guest User

Untitled

a guest
Dec 11th, 2019
136
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
R 2.54 KB | None | 0 0
  1. #exercise 1
  2. z <- runif(100, min = 0, max = 1) # выборка объемом 100
  3. z
  4. x <- z[1:49] # частичная выборка объема 49
  5. x
  6. y <- z[50:100] # частичная выборка объема 51
  7. y
  8. ks.test(x,y,alternative = "two.sided") # Критерий Смирнова с альтернативной гипотезой x <> y
  9. ks.test(x,y,alternative = "less") # Критерий Смирнова с альтернативной гипотезой x < y
  10. ks.test(x,y,alternative = "greater") # Критерий Смирнова с альтернативной гипотезой x > y
  11.  
  12. #exercise 2
  13. x<-runif(100,min=2,max=3)#выборка объемом 100
  14. a<-mean(x)#мат. ожидание
  15. s<-sd(x)#среднеквадратичное отклонение
  16. a<-seq(a-s/2,a+s/2,length=99)
  17. p<-sapply(a,function(newa) t.test(x,mu=newa, alternative="two.sided")[[3]])
  18. p#значение p-value с разными значениями a
  19. pl<-sapply(a,function(newa) t.test(x,mu=newa, alternative="less")[[3]])
  20. pl#значение p-value с разными значениями a
  21. pg<-sapply(a,function(newa) t.test(x,mu=newa, alternative="greater")[[3]])
  22. pg#значение p-value с разными значениями a
  23. matplot(a,cbind(p,pl,pg),type="l",lty=c(1,2,4),xlab="a",ylab="Level of significance",
  24.         main="The dependence of the significance level on parametr a",col=c("chartreuse1","red","blue3"))
  25. legend("topleft",legend=c("equal","less","greater"),
  26.        col=c("chartreuse1","red","blue3"),lty=1:2,cex=0.7)
  27. abline(h=c(0,0.05),lty=3)#граница принятия нулевой гипотезы
  28.  
  29. #exercise 3
  30. set.seed(9)
  31. n<-200
  32. Y<-0.2
  33. x <- rexp(n,Y)
  34. breaks = c(-Inf,0.7,1.5,2.2,2.9, 4,4.7,6.2,8.2,10.2, +Inf) # границы интервалов
  35. breaks
  36. table <- table(cut(x, breaks = breaks)) #матрица
  37. table
  38. p3 <- diff(pexp(breaks, 0.2)); p3 # вектор вероятностей
  39. sum(p3)
  40. chisq.test(x = table, p = p3)
  41.  
  42. #exercise 4
  43. n<-100
  44. a<-8
  45. sd<-2
  46. x<-rnorm(n,a,sd)
  47. m<-2:20 #  Вектор числа степеней свободы
  48. v<-x-mean(x) # центрированный вектор значений выборки  
  49. alpha<-sapply(m, function(z) ks.test(v, "pt", z)[[2]])#значение p-value в критерии Смирнова при разных m
  50. alpha
  51. plot(m, alpha, type="b", pch=20, ylim=c(0, max(alpha, 0.05)),
  52.      xlab="Number of degrees of freedom",
  53.      ylab="Significance level")
  54. abline(h=0.05, lty=2)# граница принятия нулевой гипотезы
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement