总结和更正

#-------------------第一套--------------------------#
###第一题######
## 编写函数计算h（x,n）=1+x+x^2
#####第一题############
fun1 = function(x,n){s = 0 # 定义一个变量接受求和for (i in 0:n) {s = s + x^i}message('前n项目和为：',s) # 打印结果
}
fun1(2,4) # 调用函数######第二题###############
## 编写函数计算15，21，39，45，17的乘积，分别用for循环和while循环
######第二题###############
data1 = c(15,21,39,45,17)  #定义数据
## for 循环
fun2 = function(data1){m = 1 # 这里设为1来装连乘for (i in 1:length(data1)) {m =  m*data1[i]}message('最后乘积为：',m) # 打印结果
}
fun2(data1)###########分割线################################# while 仅此条件不同而已
fun3 = function(data1){m = 1i = 1 # 这里需要先定义iwhile (i<=length(data1)) {m = m * data1[i]i =i + 1 # 需要更新变量}message('最后乘积为：',m)
}#######第三题#########
#用循环函数求出1-20之间能被4整除的数的和，并输出结果。
n = 20 # 20结尾的数
s = 0
a = c()# 定义空列装符合条件的值
for (i in 1:n){if(i %% 4 == 0){a = c(a,i)
}
}
print(sum(a))
########第四题#########
# x向量的值为100个5以内（包括5）的随机整数，
# 产生一个长度与x相等的y向量，对于x中值为3的元素，
# y的对应元素赋值为0，其余赋值为1。计算y中有多少1。
set.seed(0)
x = round(runif(100,0,5.1)) # 产生100个0-5的整数
y = c(1:length(x)) # 产生空向量
x
y[which(x!=3)]=1
y
y[which(x==3)]=0
y############第五题############
## （1）读入文件5.txt中的数据，建立y关于x1 x2 x3 x4的线性回归方程。
# 读取数据
data5 = read.csv('./5.txt',sep = '\t',header = TRUE)
# 建立模型
ress = lm(y~x1+x2+x3+x4,data = data5)
summary(ress)
## （2）采用逐步回归法建立y关于x1 x2 x3 x4的线性回归方程，并对方程和回归系数进行显著性检验。
data5.lm1 = step(ress,direction = 'both')
data5.lm1
data5.lm2 = step(ress,direction = 'backward')
data5.lm2
data5.lm3 = step(ress,direction = 'forward')
data5.lm3#############第六题############
# 使用while循环求1+2+3+…+100的和
# for循环
s = 0
for (i in 1:100) {s = s + i
}
print(s)
# while循环
i = 1
s = 0
while(i<=100){s = s + ii = i + 1
}
print(s)#############第七题############
fun7 = function(score){if(score>=90){print('等级A')} else if(score>=80){print('等级B')} else if(score>=70){print('等级C')} else if(score>=60){print('等级D')} else{print('等级E')}
}
fun7(23)######第八题######
fun8 <- function(x,u,sigma) {y = 1/(sqrt(2*pi)*sigma)*exp(-(x-u)^2/(2*sigma^2))return(y)
}fun8(1,0,1)
# 校验
dnorm(1,0,1)#####第九题######
# 用内置函数计算得到
fun99 = function(x){sin(x)
}
# 计算2-3的积分
integrate(fun99,2,3)# 先定义一个函数用于判断随机抽取的点在积分区域内或者外
fun9 = function(a){# a是一个矩阵a[2]-sin(a[1]) # y-y_count 第一个y是随机产生的y，y_count是x计算得到的
}
n = 100000 # 抽取10万次
x = runif(n,2,3)
y = runif(n,0,sin(2))
A = cbind(x,y)
b = apply(A, 1,fun9)
mean(b<=0)*(3-2)*(sin(2)-0)# 通过比例计算积分区域面积# 第九题补充 计算误差
i = 1
list_a = c()
while (i<=1000) {fun9 = function(a){# a是一个矩阵a[2]-sin(a[1]) # y-y_count 第一个y是随机产生的y，y_count是x计算得到的}n = 1000 # 抽取10万次x = runif(n,2,3)y = runif(n,0,sin(2))A = cbind(x,y)b = apply(A, 1,fun9)squre = mean(b<=0)*(3-2)*(sin(2)-0)# 通过比例计算积分区域面积error_9 = 0.5738426-squrelist_a = c(list_a,error_9)i = i + 1# print(mean(list_a))
}
aa = c(1:1000)
plot(aa,list_a)# 法2
n = 100000
x = runif(n,2,3)
mean(sin(x))*(3-2)#######第十题######
data=c(0.497,0.506,0.524,0.518,0.498,0.511,0.520,0.515,0.512)
fun6 = function(data,alph,u){x_bar = mean(data)t_value = (x_bar-u)/(sd(data)/sqrt(length(data)))print(t_value)if(abs(t_value)>= qt(1-alph/2,length(data)-1)){print('拒绝原假设')}else{print('接受原假设')}# dnorm为概率密度函数，返回给定x下的概率值。dnorm(z)等价于f(x=z)# pnorm为概率分布函数。pnorm(z)等价于P[x<=z]# qnorm为pnorm的逆运算，返回给定P下的z。# rnorm生成符合正太分布的随机数。
}fun6(data,alph = 0.05,u = 0.5)

1、编写函数计算 h(x,n)=1+x+x∧2+⋯⋯+x∧nh(x, n)=1+x+x^{\wedge} 2+\cdots \cdots+x^{\wedge} nh(x,n)=1+x+x∧2+⋯⋯+x∧n

#####第一题############
fun1 = function(x,n){s = 0 for (i in 0:n) {s = s + x^i}message('前n项目和为：',s)
}
fun1(2,4) # 调用函数

2、编写函数计算15，21，39，45，17的乘积，分别用for循环和while循环

######第二题###############
data1 = c(15,21,39,45,17)
## for 循环
fun2 = function(data1){m = 1 # 这里设为1来装连乘for (i in 1:length(data1)) {m =  m*data1[i]}message('最后乘积为：',m)
}
fun2(data1)###########分割线################################# while 仅此条件不同而已
fun3 = function(data1){m = 1i = 1while (i<=length(data1)) {m = m * data1[i]i =i + 1 # 需要更新变量}message('最后乘积为：',m)
}

3、用循环函数求出1-20之间能被4整除的数的和，并输出结果。

n = 20
s = 0
a = c() # 用来装能被4整除的数
for (i in 1:n) {if(i %% 4 == 0)a = c(a,i)
}
print(sum(a))

4、x向量的值为100个5以内（包括5）的随机整数，产生一个长度与x相等的y向量，对于x中值为3的元素，y的对应元素赋值为0，其余赋值为1。计算y中有多少1。


###########第四题###########
# 考察which函数使用
# 第一步产生x向量的值为100个5以内（包括5）的随机整数
x = round(runif(100,0,5.1))
y = c() # 产生空向量
x[which(x!=3)]=1 # x其余赋值为1
x[which(x==3)]=0# x中值为3的元素，y的对应元素赋值为0
x
# 注意这里的顺序一定不能乱
message("有",sum(x))

5、

（1）读入文件5.txt中的数据，建立y关于x1 x2 x3 x4的线性回归方程。

（2）采用逐步回归法建立y关于x1 x2 x3 x4的线性回归方程，并对方程和回归系数进行显著性检验。

############第五题############
## （1）读入文件5.txt中的数据，建立y关于x1 x2 x3 x4的线性回归方程。
# 读取数据
data5 = read.csv('./5.txt',sep = '\t',header = TRUE)
# 建立模型
ress = lm(y~x1+x2+x3+x4,data = data5)
summary(ress)
## （2）采用逐步回归法建立y关于x1 x2 x3 x4的线性回归方程，并对方程和回归系数进行显著性检验。
data5.lm1 = step(ress,direction = 'both')
data5.lm1
data5.lm2 = step(ress,direction = 'backward')
data5.lm2
data5.lm3 = step(ress,direction = 'forward')
data5.lm3

6、使用while循环求1+2+3+…+100的和。

i = 1
s = 0
while (i<=100) {s = i + s i = i+1 # 更新变量
}
print(s)

7、编写函数，对任一个百分制的成绩M，将其转换成对应的等级，具体转换规则如下：

90-100为A; 80-89为B; 70-79为C; 60-69为D; 0-59为E;

#########第七题###############
#########第七题###############
# print('请输入一个0-100的数，enter两次结束')
# A = scan() # enter两次结束
fun4 = function(A){if(90<=A) {print("等级A")} else if (80<=A) {print("等级B")}else if (70<=A) {print("等级C")}else if (60<=A) {print("等级D")}else {print("等级E")}}
fun4(60)

8、编写函数出任意给定x的值，计算正态分布密度函数N（）的值。

f(x)=12πσe−(x−μ)22σ2(−∞<x<+∞)f(x)=\frac{1}{\sqrt{2 \pi} \sigma} e^{-\frac{(x-\mu)^{2}}{2 \sigma^{2}}} \quad(-\infty<x<+\infty) f(x)=2πσ1e−2σ2(x−μ)2(−∞<x<+∞)

##########第八题##################
fun5 = function(x,u,sigma){y = 1/(sqrt(2*pi)*sigma)*(exp(-(x-u)^2/(2*sigma^2)))return(y)
}
fun5(1,0,1)
# 校验
dnorm(1,0,1)

9、用蒙特卡洛法计算的值。

∫23sin⁡xdx\int_{2}^{3} \sin x d x ∫23sinxdx

###########第九题################
# 基于蒙特卡洛求定积分的两种方法
# 法1
decisionCondition = function(a){a[2]-sin(a[1])
}
# 用于判断在面积内的函数n = 100000
x = runif(n,2,3)
y = runif(n,0,sin(2)) # 看图像改数据
A = cbind(x,y)
A
b = apply(A, 1, decisionCondition)
# length(which(b<0))
mean(b<=0)*1*sin(2) # 注意这里，将算出来的面积乘上总面积得到积分区域面积# 法2
n=10000
x = runif(n,2,3)
mean(sin(x))*1 # 乘以（b-a）

10、给定一组数，检验其均值等于某一常数（t检验）

μ0 = 0.5
H0:μ=μ0vs H1:μ≠μ0H_{0}: \mu=\mu_{0} \text { vs } H_{1}: \mu \neq \mu_{0} H0:μ=μ0 vs H1:μ=μ0
检验统计量为：
t=xˉ−μ0s/nt=\frac{\bar{x}-\mu_{0}}{s / \sqrt{n}} t=s/nxˉ−μ0
拒绝域为 ∣t∣>tc(n−1)|t|>t_{c}(n-1)∣t∣>tc(n−1) 。
X=c(0.497,0.506,0.524,0.518,0.498,0.511,0.520,0.515,0.512)\mathrm{X}=\mathrm{c}(0.497,0.506,0.524,0.518,0.498,0.511,0.520,0.515,0.512)X=c(0.497,0.506,0.524,0.518,0.498,0.511,0.520,0.515,0.512)
μ0=0.5\mu_{0}=0.5μ0=0.5

##########第十题########
data=c(0.497,0.506,0.524,0.518,0.498,0.511,0.520,0.515,0.512)
fun6 = function(data,alph,u){x_bar = mean(data)t_value = (x_bar-u)/(sd(data)/sqrt(length(data)))if(abs(t_value)>= qt(1-alph/2,length(data)-1)){print('拒绝原假设')}else{print('接受原假设')}# dnorm为概率密度函数，返回给定x下的概率值。dnorm(z)等价于f(x=z)# pnorm为概率分布函数。pnorm(z)等价于P[x<=z]# qnorm为pnorm的逆运算，返回给定P下的z。# rnorm生成符合正太分布的随机数。
}fun6(data,0.975,0.5)

R语言期末试题-重庆工商大学-统计学课程相关推荐

广西大学c语言期末试题,2006广西大学c课程考试试卷_答案.pdf
2006广西大学c课程考试试卷_答案 C语言程序设计考试试卷 2006年6月11日 × ×大学课程考试试卷考试用 (2005-- 2006学年度第 2 学期) 课程名称:C语言程序设计试卷类型:( ...
R语言期末复习资料----助力高绩点
多元统计分析及R语言建模(第四版) R语言期末复习资料第一章多元统计分析概述 1.列出常用的统计软件,说明其使用范围和各自的优缺点解: (1)SAS:组合软件系统,入门比较困难 (2)SPSS: ...
c语言描述考试答案,c语言期末试题
篇一:c语言期末试题试卷编号:10183 所属语言:Turbo C 试卷方案:fny-13 试卷总分:110分共有题型:7种 ━━━━━━━━━━━━━━ 一.填空共5题 (共计15分) ━━━ ...
R语言计算回归模型学生化残差（Studentized Residuals）实战：如果样本学生化残差（Studentized Residuals）绝对值大于3则是离群值
R语言计算回归模型学生化残差(Studentized Residuals)实战:如果样本学生化残差(Studentized Residuals)绝对值大于3则是离群值目录
c语言表达式106的结果是,云南师范大学C语言期末试题
<云南师范大学C语言期末试题>由会员分享,可在线阅读,更多相关<云南师范大学C语言期末试题(10页珍藏版)>请在人人文库网上搜索. 1.云南师范大学2009-2010学年下学期 ...
史上R语言最强--资源（免费课程、书籍、教程和各种高级图表）
史上R语言最强–资源(免费课程.书籍.教程和各种高级图表) 下面是一些免费R的书籍.教程.软件包.图表2单和其他材料,可以在网上学习编程和改进工作流程.有单独的概述 Python resources, ...
湖南科技大学c语言试题,湖南科技大学C语言期末试题库.doc
湖南科技大学C语言期末试题库选择题(共题),成绩() 1.C语言所提供的基本数据类型包括字符型.整型.双精度型.单精度型和____: A) 指针型 B) 结构型 C) 数组型 D) 枚举类型正确 ...
湖南科技大学c语言试题,湖南科技大学c语言期末试题库
湖南科技大学c语言期末试题库志伟分享左凯小编 1 选择题(共题),成绩() 1. C 语言所提供的基本数据类型包括字符型.整型. 双精度型.单精度型和____: A) 指针型 B) 结构型 C) ...
R语言自定义编写函数生成学生化残差的直方图（dist of studentized residuals），并叠加标准正太曲线、核密度曲线、轴须图rug曲线、检验模型是否满足正态性（normality)
R语言自定义编写函数生成学生化残差的直方图(Distribution of studentized residuals),并叠加标准正太曲线.核密度曲线.轴须图rug曲线.检验模型是否满足正态性(no ...
c语言电子期末考试题,山东师范大学11-12电子C语言期末试题
山东师范大学2011-2012学年第1学期期末考试试题 (时间:120分钟共100分) 课程编号:080920201.080940201 课程名称:C语言程序设计试题类别: A考试类型:闭卷 Pa ...

R语言期末试题-重庆工商大学-统计学课程

总结和更正

1、编写函数计算 h(x,n)=1+x+x∧2+⋯⋯+x∧nh(x, n)=1+x+x^{\wedge} 2+\cdots \cdots+x^{\wedge} nh(x,n)=1+x+x∧2+⋯⋯+x∧n

2、编写函数计算15，21，39，45，17的乘积，分别用for循环和while循环

3、用循环函数求出1-20之间能被4整除的数的和，并输出结果。

4、x向量的值为100个5以内（包括5）的随机整数，产生一个长度与x相等的y向量，对于x中值为3的元素，y的对应元素赋值为0，其余赋值为1。计算y中有多少1。

5、

（1）读入文件5.txt中的数据，建立y关于x1 x2 x3 x4的线性回归方程。

（2）采用逐步回归法建立y关于x1 x2 x3 x4的线性回归方程，并对方程和回归系数进行显著性检验。

6、使用while循环求1+2+3+…+100的和。

7、编写函数，对任一个百分制的成绩M，将其转换成对应的等级，具体转换规则如下：

8、编写函数出任意给定x的值，计算正态分布密度函数N（）的值。

9、用蒙特卡洛法计算的值。

10、给定一组数，检验其均值等于某一常数（t检验）

R语言期末试题-重庆工商大学-统计学课程相关推荐

最新文章

热门文章

R语言期末试题-重庆工商大学-统计学课程

总结和更正

1、编写函数计算 h(x,n)=1+x+x∧2+⋯⋯+x∧nh(x, n)=1+x+x^{\wedge} 2+\cdots \cdots+x^{\wedge} nh(x,n)=1+x+x∧2+⋯⋯+x∧n

2、编写函数计算15，21，39，45，17的乘积，分别用for循环和while循环

3、用循环函数求出1-20之间能被4整除的数的和，并输出结果。

4、x向量的值为100个5以内（包括5）的随机整数，产生一个长度与x相等的y向量，对于x中值为3的元素，y的对应元素赋值为0，其余赋值为1。计算y中有多少1。

5、

（1）读入文件5.txt中的数据，建立y关于x1 x2 x3 x4的线性回归方程。

（2）采用逐步回归法建立y关于x1 x2 x3 x4的线性回归方程，并对方程和回归系数进行显著性检验。

6、使用while循环求1+2+3+…+100的和。

7、编写函数，对任一个百分制的成绩M，将其转换成对应的等级，具体转换规则如下：

8、编写函数出任意给定x的值，计算正态分布密度函数N（）的值。

9、用蒙特卡洛法计算 的值。

10、给定一组数，检验其均值等于某一常数（t检验）

R语言期末试题-重庆工商大学-统计学课程相关推荐

最新文章

热门文章

9、用蒙特卡洛法计算的值。