[R-br] Problema com o constrOptim
Thiago Touya
thiagotouya em gmail.com
Quarta Setembro 21 23:34:00 BRT 2011
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)
-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110921/2636f37b/attachment.html>
-------------- Próxima Parte ----------
Um anexo não-texto foi limpo...
Nome: datasets-cleveland-clinic-avc.csv
Tipo: text/csv
Tamanho: 18976 bytes
Descrição: não disponível
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20110921/2636f37b/attachment.csv>
Mais detalhes sobre a lista de discussão R-br