试验三 描述统计

1.熟悉课本P277-279的内容。
2.使用R语言完成课本P24 第4题和第5题。；

setwd("C:/Users/zhang/Desktop")
a<-read.csv("fishlenth.csv",header=F,col.names="体长")
#2-第四题
max(a)
min(a)
range<-max(a)-min(a)
breakspoints<-seq(0,300,20)
labels<-as.character(seq(0,300,20)[-1])
labels
c<-table(cut(a$体长,breaks = breakspoints,labels = labels),exclude = NULL)
write.csv(c,"1.csv")
jpeg("1.jpg") 
hist(a$体长, breaks = breakspoints, col = 'red', xlab = '体长', main = '直方图')
dev.off()  
dev.off()  

#2-第五题
a<-c(a$体长)
mean<-mean(a)
median<-median(a)
mode<-which(table(a)==max(table(a)))
var<-var(a)
range<-max(a)-min(a)
std.dev<-sd(a)
C.V<-100*std.dev/mean
middle.quartile<-quantile(a,3/4)-quantile(a,1/4)
mean;median;mode;var;range;std.dev;C.V;middle.quartile

3.使用R语言中的grubbs.test函数分析课本P24第4题的数据中是否有异常值，如果有异常值，列出各异常值。

#3
install.packages("outliers")
library(outliers)
grubbs.test(a)
grubbs.test(a,opposite = T)
grubbs.test(a,opposite = F)

4.对自由度为20的t分布，分析：
（1）当t=10时的概率密度函数值和概率值；
（2）P=0.05和0.01时对应的t值；

#4
df <- 20
t_value <- 10
pdf_value <- dt(t_value, df)  
cdf_value <- pt(t_value, df)   
cat("当t=10时的概率密度函数值:", pdf_value, "\n")
cat("当t=10时的概率值:", cdf_value, "\n\n")
p_values <- c(0.05, 0.01)
t_values <- sapply(p_values, function(p) qt(p, df))
t_values<-qt(p_values, df)
n=0
for (p_value in p_values) {
  n=n+1
  cat(paste("P=", p_value, "时对应的t值:", t_values[n], "\n"))
}