Advertisement
Guest User

Untitled

a guest
May 22nd, 2019
91
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 3.22 KB | None | 0 0
  1. # Zadania
  2.  
  3. set.seed(100)
  4. #####
  5.  
  6. X <- rnorm(100,2,0.5)
  7. alpha1 <- 0.05 #poziom istotności
  8. statystyka1 <- ((mean(X) - 2.05)/0.5)*10
  9. #1a
  10. #testowanie 1 sposobem
  11. S1a <- qnorm(0.05) # w ten sposób szukamy obszaru krytycznego
  12. S1a > stat # False, zatem nie ma podstaw do odrzucenia H na rzecz K
  13.  
  14. #p - wartość
  15.  
  16. pwartosc1a <- pnorm(statystyka1)
  17. pwartosc1a < alpha1 #False, zatem nie ma podstaw do odrzucenia H na rzecz K
  18.  
  19. #1b
  20. #testowanie 1 sposobem
  21. S1ba <- qnorm(0.025)
  22. S1bb <- qnorm(0.975)
  23.  
  24. S1ba > statystyka1 || S1bb < statystyka1 #zatem nie ma podstaw do odrzucenia H na rzecz K
  25.  
  26. #p - wartość
  27.  
  28. pwartosc1b <- 2*min(pnorm(statystyka1),1-pnorm(statystyka1))
  29.  
  30. pwartosc1b < alpha1 #zatem nie ma podstaw do odrzucenia H na rzecz K
  31. #####
  32.  
  33. Y <- rnorm(100,5,2)
  34. alpha2 <- 0.01
  35. statystyka2 <- (mean(Y)-5.5)*10/sd(Y)
  36. #2a
  37. #testowanie 1 sposobem
  38. S2a <- qt(0.99, 99)
  39. S2a < statystyka2 #zatem nie ma podstaw do odrzucenia H na rzecz K
  40.  
  41. #p - wartość
  42.  
  43. pwartosc2a <- 1 - pt(statystyka2, 99)
  44. pwartosc2a < alpha2 #zatem nie ma podstaw do odrzucenia H na rzecz K
  45.  
  46. #2b
  47. #testowanie 1 sposobem
  48. S2ba <- qt(alpha2/2,99)
  49. S2bb <- qt(1-alpha2,99)
  50.  
  51. S2ba > statystyka2 || S2bb < statystyka2 #zatem odrzucamy H na rzecz K
  52.  
  53. #p - wartość
  54.  
  55. pwartosc2b <- 2*min(pt(statystyka2,99),1-pt(statystyka2,99))
  56. pwartosc2b < alpha2 #zatem odrzucamy H na rzecz K
  57. #####
  58.  
  59. Z <- rnorm(100,20,5)
  60. alpha3 <- 0.1
  61. statystyka3 <- 99*var(Z)/30
  62. #3a
  63. #testowanie 1 sposobem
  64. S3a <- qchisq(1-alpha3,99)
  65. S3a < statystyka3 #zatem nie ma podstaw do odrzucenia H na rzecz K
  66.  
  67. #p - wartość
  68.  
  69. pwartosc3a <- 1 - pchisq(statystyka3,99)
  70. pwartosc3a < alpha3 #zatem nie ma podstaw do odrzucenia H na rzecz K
  71.  
  72. #3b
  73. #testowanie 1 sposobem
  74. S3ba <- qchisq(alpha3,99)
  75. S3bb <- qchisq(1-alpha3,99)
  76.  
  77. S3ba > statystyka3 || S3bb < statystyka3 #zatem nie ma podstaw do odrzucenia H na rzecz K
  78.  
  79. #p - wartość
  80.  
  81. pwartosc3b <- 2*min(pchisq(statystyka3,99),1-pchisq(statystyka3,99))
  82. pwartosc3b < alpha3 #zatem nie ma podstaw do odrzucenia H na rzecz K
  83. #####
  84.  
  85. U <- rexp(100,0.1)
  86. alpha4 <- 0.1
  87. statystyka4 <- 50*mean(U)
  88. #4a
  89. #testowanie 1 sposobem
  90. S4a <- qchisq(alpha4,200)
  91. S4a > statystyka4 #zatem nie ma podstaw do odrzucenia H na rzecz K
  92.  
  93. #p - wartość
  94.  
  95. pwartosc4a <- pchisq(statystyka4,200)
  96. pwartosc4a < alpha4 #zatem nie ma podstaw do odrzucenia H na rzecz K
  97.  
  98. #4b
  99. #testowanie 1 sposobem
  100. S4ba <- qchisq(alpha4/2,200)
  101. S4bb <- qchisq(1-alpha4/2,200)
  102.  
  103. S4ba > statystyka4 || S4bb < statystyka4 #zatem odrzucamy H na rzecz K
  104.  
  105. #p - wartość
  106.  
  107. pwartosc4b <- 2*min(pchisq(statystyka4,200),1-pchisq(statystyka4,200))
  108. pwartosc4b < alpha4 #zatem odrzucamy H na rzecz K
  109. #####
  110.  
  111. V <- rbinom(100,1,0.6)
  112. alpha5 <- 0.05
  113. statystyka5 <- (mean(V)-0.5)*10/sqrt(mean(V)*(1-mean(V)))
  114. #5a
  115. #testowanie 1 sposobem
  116. S5a <- qnorm(1-alpha5)
  117.  
  118. S5a < statystyka5 #zatem odrzucamy H na rzecz K
  119.  
  120. # p - wartość
  121.  
  122. pwartosc5a <- 1 - pnorm(statystyka5)
  123. pwartosc5a < alpha5 #zatem odrzucamy H na rzecz K
  124.  
  125. #5b
  126. #testowanie 1 sposobem
  127. S5ba <- qnorm(alpha5/2)
  128. S5bb <- qnorm(1-alpha5/2)
  129.  
  130. S5ba > statystyka5 || S5bb < statystyka5 #zatem odrzucamy H na rzecz K
  131.  
  132. #p - wartość
  133.  
  134. pwartosc5b <- 2*min(pnorm(statystyka5),1-pnorm(statystyka5))
  135. pwartosc5b < alpha5 #zatem odrzucamy H na rzecz K
  136. #####
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement