- 1、本文档共10页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
R基础课后(非参数统计)
#A.1
AGE=c(18,23,22,21,20,19,20,20,20)
first3=sort(AGE)[(length(AGE)-2):length(AGE)]
delete.1=order(AGE)[(length(AGE)-2):length(AGE)]
except=AGE[-delete.1]
c(except[1:2],19,except[3:length(except)])
except[2]=22
#A.2
a1=rep(1:3,rep(2,3))
a2=c(1,8,10,11)
a3=seq(1,30,length(a2))
a4=seq(1,5,2)
a4=c(a1,a4,rep(0,2))
a5=2:10
a6=c(a2,a3[-(1:3)],a4)
a7=c(rep(1,10),rep(0,8))
#1.3
library(MASS)
data(geyser)
a=geyser
a1=subset(a,waiting70)
a1=geyser[geyser$waiting70,]#两个a1等价
a2=subset(a,waiting70waiting!=57)
a2=geyser[geyser$waiting70geyser$waiting!=57,]#两个a2等价
a3=subset(a,waiting70,c(duration))
a4=subset(a,duration70,c(waiting))
#A.3
library(MASS)
a=geyser
attach(a)
b1=waiting[waiting70]
b2=waiting[waiting70waiting!=57]
b3=duration[waiting70]
b4=waiting[duration70]
#1.4
x=c(0,1,1,2,3,4)
num=function(x){
r=0
p=c(rep(0,length(x)-1))
q=c(rep(0,length(x)-1))
for(i in 1:(length(x)-1))
{
for(j in (i+1):length(x)){
p[i]=p[i]+I(x[i]x[j])
q[i]=q[i]+I(x[i]x[j])
}
r=r+p[i]-q[i]
}
r
}
h=rep(0,1000)
for(i in 1:length(h)){
x=runif(5,-5,5)
h[i]=num(x)
}
y=as.factor(h)
yy=levels(y)
a=rep(0,11)#记数 记录-10:10:2各种结果出现的次数
p=rep(0,11)
for(i in 1:11){
a[i]=sum(h==(-12+2*i))
p[i]=a[i]/length(h)
}
a
p
hist(h)
h=rep(0,10000)
for(i in 1:length(h)){
x=sample(1:5,5,replace=T)#-5到5中间的整数随机抽样
h[i]=num(x)
}
y=as.factor(h)
yy=levels(y)
a=rep(0,length(yy))#记数 记录-10:10各种结果出现的次数
p=rep(0,length(yy))
for(i in 1:length(yy)){
a[i]=sum(h==(-11+i))
p[i]=a[i]/length(h)
}
a
p
hist(h)
#当随机取10000次的一个结果 a=[71 321 774 1255 1637 1825 1684 1256 743 338 96]
# p=[0.0071 0.0321 0.0774 0.1255 0.1637 0.1825 0.1684 0.1256 0.0743 0.0338 0.0096]
#当随机取十万次数据的一个结果
#a=[795 3421 7553 12521 16771 18180 16538 12553 7418 3400 850]
#p=[0.00795 0.03421 0.07553 0.12521 0.16771 0.18180 0.16538 0.12553 0.07418 0.03400 0.00850]
#1.5
uniroot(f=function(x) 2*x^3-4*x^2+3*x-6, interval=c(-10,10))
f=function(x){2*x^3-4*x^2+3*x-6}
f(0)
a=-10
b=10
root=function(a,b)
{
c=(a+b)/2;
whi
文档评论(0)