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)
Você deve se certificar de que os pacotes que possuem as funções que você usa foram todos carregados. Para ajudar, dê um RSiteSearch("constrOptim") RSiteSearch("grad") como ponto de partida para a sua investigação. Nestes caso e super importante fornecer o sessionInfo(), pois pode ser problema de versão coisa e tal. Envite o envio de arquivos em anexo pois isso não faz o seu código ser reproduzível com um copia e cola. Consulte o link no rodapé da mensagem para conhecer alternativas de envio de dados para a lista. À disposição. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes@ufpr.br twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================
participantes (2)
-
Thiago Touya -
Walmes Zeviani