Advertisement
Guest User

estadistica aplicada - metodos estadisticos - clase 5

a guest
Apr 25th, 2017
70
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
R 1.69 KB | None | 0 0
  1. data <- read.delim("clipboard")
  2. attach(data)
  3. detach(data)
  4. args(subset)
  5. typeof(data)
  6. class(data)
  7. data$Universidad == "UNCI"
  8. unci <- subset(data, Universidad == "UNCI")
  9. attach(unci)
  10. unci.t.test <- t.test(Nota.2015.2,conf.level = 0.92)
  11. unci.t.test$conf.int[[2]]
  12.  
  13. attach(unci.t.test)
  14. detach(unci.t.test)
  15. icm <- function(x, nc, N){
  16.   alfa <- 1 - nc
  17.   xbar <- mean(x)
  18.   s <- sd(x)
  19.   n <- length(x)
  20.   ee <- s / sqrt(n)
  21.   fcpf <- (N-n)/(N-1)
  22.   fcpf.sqrt <- sqrt(fcpf)
  23.   LI <- xbar - qt(1 - (alfa/2), n - 1) * ee * fcpf.sqrt
  24.   LS <- xbar + qt(1 - (alfa/2), n - 1) * ee * fcpf.sqrt
  25.   L <- c(LI, LS)
  26.   print(L)
  27. }
  28.  
  29. icm(Nota.2015.2, 0.92, 100)
  30.  
  31.  
  32. install.packages("EnvStats")
  33. library(EnvStats)
  34.  
  35. varTest(Nota.2015.2, conf.level = 0.92)$conf
  36. sqrt(varTest(Nota.2015.2, conf.level = 0.92)$conf)
  37.  
  38. #####################
  39. ####################
  40.  
  41. t.test(Nota.2015.2, Nota.2015.1, paired = T, conf.level = 0.92)$conf
  42. t.test(Nota.2015.1, Nota.2015.2, paired = T, conf.level = 0.92)$conf
  43.  
  44. ####################
  45. ###################
  46.  
  47. install.packages("DescTools")
  48. library(DescTools)
  49.  
  50. CoefVar(Nota.2015.2,conf.level = 0.92)
  51. CoefVar(Nota.2015.1,conf.level = 0.92)
  52. ?CoefVar
  53.  
  54. ##################
  55. ##################
  56.  
  57. detach(unci)
  58. detach(data)
  59. attach(data)
  60. no_unci <- subset(data, Universidad != "UNCI")
  61.  
  62. attach(no_unci)
  63. detach(no_unci)
  64. var.test(Nota.2015.2~Universidad, conf.level = 0.92)
  65.  
  66.  
  67. ##################
  68. #################
  69.  
  70. #by default: paired = FALSE, var.equal = F
  71. t.test(Nota.2015.2~Universidad, conf.level = 0.92, paired = F, var.equal = T)$conf
  72.  
  73.  
  74. install.packages("BSDA")
  75. library(BSDA)
  76. ?z.test
  77.  
  78. z.test(x = Nota.2015.2[Universidad = "UNP"],sigma.x = 1, sigma.y = 1, conf.level = 0.92)
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement