|
谢谢!
您好,想问下自己给出copula函数
代码如下
density.Gumbel<-function(u,v,alpha)
{
temp<-(-log(u))^alpha+(-log(v))^alpha
((-log(u))^(alpha-1)*(-log(v))^(alpha-1)*temp^(-2+(1/alpha))*temp^((1/alpha)-1+alpha))/(u*v*exp(temp^(1/alpha)))
}
likeli.Gumbel<-function(alpha)
{
f<-sum(log(eval(call(density.Gumbel,FSh,FSz.alpha)))
f
}
nlminb(objective=likeli.Gumbel,start=2) #求极值
alpha<-nlminb(objective=likeli.Gumbel,start=2)$parameters
paste('alpha=',alpha)#打印参数
这里提示call第一个参数必须是字符型,请问怎么解决?
|