- 阅读权限
- 255
- 威望
- 0 级
- 论坛币
- 3533 个
- 通用积分
- 1.8600
- 学术水平
- 5 点
- 热心指数
- 9 点
- 信用等级
- 2 点
- 经验
- 7070 点
- 帖子
- 336
- 精华
- 0
- 在线时间
- 151 小时
- 注册时间
- 2009-5-14
- 最后登录
- 2019-11-22
博士生
还不是VIP/贵宾
- 威望
- 0 级
- 论坛币
- 3533 个
- 通用积分
- 1.8600
- 学术水平
- 5 点
- 热心指数
- 9 点
- 信用等级
- 2 点
- 经验
- 7070 点
- 帖子
- 336
- 精华
- 0
- 在线时间
- 151 小时
- 注册时间
- 2009-5-14
- 最后登录
- 2019-11-22
| 开心 2018-1-7 20:39:21 |
---|
签到天数: 5 天 连续签到: 1 天 [LV.2]偶尔看看I
|
5论坛币
大侠:
我的数据框是由256列构成的(名称为bb1,bb2, bb3,..., bb256),我要对每一列分别作极大似然估计,用命令一个一个处理太费时间了,想请教一下怎样用for循环啊,我写个代码,但不知道错在哪里,
for(i in 6:256(bb)){
dd.res<-nlminb(c(4,2),LL,data=bb[i],lower=c(0.0001,-Inf,0.0001,Inf,-Inf,0.0001),upper=c(15,Inf,Inf,Inf,Inf))
+ dd.res$par
+ }
错误: 不适用于非函数
谢谢! |
最佳答案
tf2000 查看完整内容
首先要解决一个问题:在循环语句中的那个命令
nlminb(c(4,2),LL,data=bb,lower=c(0.0001,-Inf,0.0001,Inf,-Inf,0.0001),upper=c(15,Inf,Inf,Inf,Inf))
是不是正确的?
如果是正确的,那就是出现在循环语句中
请检查“bb”,这个似乎不对,应该是bb[,i](如果是用列数据的话)
另外dd.res也应该是个list吧,似乎没有定义
回帖推荐
tf2000 发表于2楼
查看完整内容
首先要解决一个问题:在循环语句中的那个命令
nlminb(c(4,2),LL,data=bb,lower=c(0.0001,-Inf,0.0001,Inf,-Inf,0.0001),upper=c(15,Inf,Inf,Inf,Inf))
是不是正确的?
如果是正确的,那就是出现在循环语句中
请检查“bb”,这个似乎不对,应该是bb[,i](如果是用列数据的话)
另外dd.res也应该是个list吧,似乎没有定义
tf2000 发表于10楼
查看完整内容
还有一个简单的方法,用apply函数
例如
apply(bb,2,FUNCTION,...)
具体用法参见help
zjj8211 发表于9楼
查看完整内容
# define likelihood function(normal distribution)
log.likelihood
|