楼主: conan092331
1803 2

[问答] 自编的函数,引用函数得到的结果是NULL,直接复制函数内容出现正常结果,为什么啊? [推广有奖]

  • 0关注
  • 1粉丝

高中生

65%

还不是VIP/贵宾

-

威望
0
论坛币
11 个
通用积分
0
学术水平
0 点
热心指数
0 点
信用等级
0 点
经验
499 点
帖子
49
精华
0
在线时间
4 小时
注册时间
2015-10-19
最后登录
2015-11-10

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群
赵安豆老师微信:zhaoandou666

经管之家联合CDA

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

感谢您参与论坛问题回答

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

+2 论坛币
自己编写了一个函数叫DCSIS,生成了随机数进行模拟,直接调用函数时得到的结果w是NULL,但是如果直接把DCSIS函数的具体内容用到生成的随机数的话得到的w就是正常的结果,请问这是为什么?是函数调用出现了问题还是函数编写有问题?还有就是同样的函数和编写方法,matlab里运行比R快很多,这是为什么啊?


编写的函数如下:

DCSIS<-function(x,y)

{

x=t(t(x));

y=t(t(y));

n<-nrow(x);

p<-ncol(x);

m<-nrow(y);

if (m!=n)  

print("Error:matrix of x and y do not match")

source("dcov.r")      这是自己编写的另一个函数,运行正常,没有问题

dcov3=dcov(y,y);

w=rep(0,p);

for (k in 1:p)

{ dcov1=dcov(x[,k],y);

   dcov2=dcov(x[,k],x[,k]);

   w[k]=dcov1/sqrt(dcov2)/sqrt(dcov3);}

}

直接调用函数得到w是NULL

n = 100;

p = 2000;

library("mvtnorm")

x<-rmvnorm(n,rep(0,p),diag(p));   

beta<-c(1,1,1,rep(0,p-3));

y=x%*%beta+rnorm(n);

source("DCSIS.r")

w=DCSIS(x,y);



直接具体求解,得到的w是正常的结果

n = 100;

p = 2000;

library("mvtnorm")

x<-rmvnorm(n,rep(0,p),diag(p));   

beta<-c(1,1,1,rep(0,p-3));

y=x%*%beta+rnorm(n);

source("dcov.r")      

dcov3=dcov(y,y);

w=rep(0,p);

for (k in 1:p)

{ dcov1=dcov(x[,k],y);

   dcov2=dcov(x[,k],x[,k]);

   w[k]=dcov1/sqrt(dcov2)/sqrt(dcov3);}





二维码

扫码加我 拉你入群

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

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

关键词:null MATLAB atlab matla CSIS matlab

沙发
jinkelazzz 发表于 2015-10-19 11:00:38 |只看作者 |坛友微信交流群
最后少了1个w和一个大括号?

使用道具

藤椅
conan092331 发表于 2015-10-19 18:12:24 |只看作者 |坛友微信交流群
jinkelazzz 发表于 2015-10-19 11:00
最后少了1个w和一个大括号?
是说自编函数部分吗?以前编的最后也不会在写w,都能正常用。而且加上w也试过了,出来的还是NULL

使用道具

您需要登录后才可以回帖 登录 | 我要注册

本版微信群
加好友,备注cda
拉您进交流群

京ICP备16021002-2号 京B2-20170662号 京公网安备 11010802022788号 论坛法律顾问:王进律师 知识产权保护声明   免责及隐私声明

GMT+8, 2024-4-28 12:01