楼主: bbyyoo
1848 2

R 自定义似然函数 出错,帮忙看看 [推广有奖]

  • 0关注
  • 0粉丝

PHD

已卖:1289份资源

博士生

21%

还不是VIP/贵宾

-

威望
0
论坛币
3571 个
通用积分
3.0045
学术水平
5 点
热心指数
9 点
信用等级
2 点
经验
7070 点
帖子
336
精华
0
在线时间
151 小时
注册时间
2009-5-14
最后登录
2019-11-22

楼主
bbyyoo 发表于 2009-7-28 22:52:30 |AI写论文

+2 论坛币
k人 参与回答

经管之家送您一份

应届毕业生专属福利!

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

经管之家联合CDA

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

感谢您参与论坛问题回答

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

+2 论坛币
在用R做极大似然估计是出错了,哪位大侠帮忙看看错误在那里?谢谢
源代码:
>  LL<-function(params,data){
+ t1<-function(params,data){
+ t1<-params[1]*params[2](data)^(params[1]-1)*(1-data^params[1])^(params[2])
+ }
+ II<-sum(log(t1))
+  return(-II)
+ }
> dd.res<-nlminb(c(0.1,0.1),LL,data=bb1,lower=c(0.0001,-Inf,0.0001,Inf,-Inf,0.0001),upper=c(15,Inf,Inf,Inf,Inf))
错误于log(t1) : 数学函数中用了非数值变元


再往后就不能继续了。
二维码

扫码加我 拉你入群

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

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

关键词:似然函数 自定义 极大似然估计 似然估计 极大似然 函数 帮忙 定义

回帖推荐

vrooadk 发表于2楼  查看完整内容

did you miss a "*" here ? t1

本帖被以下文库推荐

沙发
vrooadk 发表于 2009-7-28 23:18:11
did you miss a  "*" here ?

t1<-params[1]*params[2]     *      (data)^(params[1]-1)*(1-data^params[1])^(params[2])
已有 1 人评分论坛币 热心指数 收起 理由
crystal8832 + 10 + 1 精彩帖子

总评分: 论坛币 + 10  热心指数 + 1   查看全部评分

藤椅
bbyyoo 发表于 2009-7-28 23:33:26
太感谢了,就是少了个乘号,这下可以了。

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

本版微信群
加好友,备注cda
拉您进交流群
GMT+8, 2025-12-28 15:46