r语言统计分析第五次


# 6.1
data1 <- c(27, 18, 15, 24, 36, 30)
ks.test(data1,y = 'runif')


# 6.2
data2 <- c(36,36,37,38,40,42,43,43,44,45,48,48,50,50,51,
           52,53,54,54,56,57,57,57,58,58,58,58,58,59,60,
           61,61,61,62,62,63,63,65,66,68,68,70,73,73,75)
shapiro.test(data2)


#6.3
data3 <- data.frame(nrow = c(1,1,1,2,2,2),
                    ncol = c(1,2,3,1,2,3),
                    freq = c(341,405,105,103,11,15))
mid_tab = xtabs(freq~nrow+ncol,data = data3)
mid_table<-matrix(c(341,405,105,103,11,15),nrow=2,ncol=3)
chisq.test(mid_tab)
chisq.test(mid_table)


# 6.4
X<-c(25.6, 22.2, 28.0, 29.8, 24.4, 30.0, 29.0, 27.5, 25.0, 27.7,
       23.0, 32.2, 28.8, 28.0, 31.5, 25.9, 20.6, 21.2, 22.0, 21.2)
A<-factor(rep(1:5, each=4))
data4<-data.frame(X, A)
aov.mis<-aov(X~A, data=data4)
summary(aov.mis)

plot(data4$X~data4$A)

# # summary(fit)
# library(gplots)
# plotmeans(X ~ A, data=data4)
# 
plot(TukeyHSD(aov.mis))

pairwise.t.test(X, A, p.adjust.method="none")
pairwise.t.test(X, A, p.adjust.method="bonferroni")

# p.adjust(p, method = p.adjust.methods, n = length(p))

# 6.5
group_<-rep(c("A","B","C"),each=8)
Weight_A <- c(15,13,11,12,12,16,14,17,17,16,
                    18,18,21,22,19,18,22,24,20,23,
                    25,27,30,32)
Y <-c(85,83,65,76,80,91,84,90,97,90,
                     100,95,103,106,99,94,89,91,83,
                     95,100,102,105,110)

data5<-data.frame(group_,Weight_A,Y)


res_5 <- aov(Y ~ Weight_A+group_ , data=data5)
summary(res_5)



############################  ... 
my_own <- function(x){
  res <- mean(x)
  return(res)
}
my_own(c(1:100))
my_own(c(1:100,NA,NA,4,5))


my_own2 <- function(x,...){
  res <- mean(x,...)
  return(res)
}
my_own2(c(1:100,NA,NA,4,5),na.rm =T)
r语言 统计分析 第五次

r语言统计分析第五次
