R软件 请大神指教><

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

1、生成一个10行10列的矩阵，其元素为1:100. 对所生成的矩阵分别按三种循环方式求所有元素的和。

2、基于二分法，用while循环的方式编写函数求方程                              在区间[1,2]内的根，精度要求为10-6。

   3、编写两样本T统计量的函数

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：r软件 while T统计量 Hil 统计量 二分法 软件 统计 样本 元素

参考薛毅老师的教材
第二个问题，在函数调用时f改成具体函数形式，a=1,b=2即可

1. fzero <- function(f, a, b, eps=1e-6){
   
2. if (f(a)*f(b)>0)
   
3. list(fail="finding root is fail!")
   
4. else{
   
5. repeat {
   
6. if (abs(b-a)<eps) break
   
7. x <- (a+b)/2
   
8. if (f(a)*f(x)<0) b<-x else a<-x
   
9. }
   
10. list(root=(a+b)/2, fun=f(x))
    
11. }
    
12. }

复制代码

计算两样本的t统计量

1. twosam <- function(y1, y2) {
   
2. n1 <- length(y1); n2 <- length(y2)
   
3. yb1 <- mean(y1); yb2 <- mean(y2)
   
4. s1 <- var(y1); s2 <- var(y2)
   
5. s <- ((n1-1)*s1 + (n2-1)*s2)/(n1+n2-2)
   
6. (yb1 - yb2)/sqrt(s*(1/n1 + 1/n2))
   
7. }

复制代码

第一个问题，第一种方法for循环语句
> abc<-matrix(c(1:100),10,10)
> tem<-0
> for (i  in 1:dim(abc)[1]){for(j in 1:dim(abc)[2]){tem<-tem+abc[i,j]}}

第二种方法while语句
> tem=0
> for (i  in 1:dim(abc)[2]){while(j <11){tem<-tem+abc[i,j];j=j+1;};j=1;}
第三种方法repeat语句
> tem=0
> for (i  in 1:dim(abc)[2]){repeat {tem<-tem+abc[i,j];j=j+1;if(j>10)break};j=1;}
> tem

victorchan0633 发表于 2015-12-4 16:43
参考薛毅老师的教材
第二个问题，在函数调用时f改成具体函数形式，a=1,b=2即可

人好好 感激不尽><