Boa noite,
estou com o seguinte problema
"Erro em dR(theta, theta.old, ...) :
não foi possível encontrar a função "grad" "
quando rodo o programa abaixo.
require(survival)
dados<-read.table('datasets-cleveland-clinic-avc.csv',header=T,sep=";")
dados
values<-dados$OPMOS
value<-as.numeric(values)
value
i<-c(1:310)
i
ajustando<-survreg(Surv(value)~1, dist='lognormal')
ajustando
verossi<-function(param){
mi<-param[1]
sigma<-param[2]
lv<-log((1/sqrt(2*3.1415)*value*sigma)*exp(-0.5*((log(value)-mi)/sigma)^2))
sum(-lv)
}
ValoresIniciais<-c(1,5)
est<-constrOptim(ValoresIniciais,verossi,NULL,method="BFGS",
ui=rbind(c(1,0),
c(0,1)),
ci=c(0,0.0001),
hessian=FALSE)
est
Alguém poderia me ajudar, pois no help é dito que se nao for usar o gradiente deve-se colocar "NULL" na posição correspondente.
O banco de dados vai em anexo.
--
Att,
Thiago Morais de Carvalho (61 - 8169-5094)