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