Try
RN()
- Phil Spector
Statistical Computing Facility
Department of Statistics
UC Berkeley
spec...@stat.berkeley.edu
On Sat, 19 Mar 2011, garciap wrote:
Dear people,
I'm trying to do some analysis of a data using the models by Royle & Donazio
in their fantastic book, particular the following function:
http://www.mbr-pwrc.usgs.gov/pubanalysis/roylebook/panel4pt1.fn
that applied to my data and in the console is as follows:
`desman.y` <- structure(c(3L,4L,3L,2L,1L), .Names = c("1", "2", "3", "4",
"5"))
RN<-function(y=desman.y,J=5,nsites=39,Nmax=100){
+ ####
+ lik<-function(parms){
+ r<-expit(parms[1])
+ lambda<-exp(parms[2])
+ pvec<-1-(1-r)^(0:Nmax)
+ gN<-dpois(0:Nmax,lambda)
+ gN<-gN/sum(gN)
+ lik<-rep(NA,nsites)
+ for(i in 1:nsites){
+ lik[i]<-sum(dbinom(y[i],J,pvec)*gN)
+ }
+ -1*sum(log(lik))
+ }
+ ####
+
+ tmp<-nlm(lik,c(0,0),hessian=TRUE)
+ ests<-tmp$estimate
+ aic<-tmp$minimum*2 + 2*length(ests)
+ se<- sqrt(diag(solve(tmp$hessian)))
+ list(ests=ests,se=se,aic=aic)
+ }
but when I try to see the results (i.e., the list in the function), I've the
follow errors:
summary(RN)
Error in object[[i]] : object of type 'closure' is not subsettable
or obtains something like:
print(RN)
function(y=desman.y,J=5,nsites=39,Nmax=100){
####
lik<-function(parms){
r<-expit(parms[1])
lambda<-exp(parms[2])
pvec<-1-(1-r)^(0:Nmax)
gN<-dpois(0:Nmax,lambda)
gN<-gN/sum(gN)
lik<-rep(NA,nsites)
for(i in 1:nsites){
lik[i]<-sum(dbinom(y[i],J,pvec)*gN)
}
-1*sum(log(lik))
}
####
tmp<-nlm(lik,c(0,0),hessian=TRUE)
ests<-tmp$estimate
aic<-tmp$minimum*2 + 2*length(ests)
se<- sqrt(diag(solve(tmp$hessian)))
list(ests=ests,se=se,aic=aic)
}
this second one seems like the analyses are not done at all.
Any suggestion? Many thanks,
Pablo
--
View this message in context:
http://r.789695.n4.nabble.com/problem-running-a-function-tp3390117p3390117.html
Sent from the R help mailing list archive at Nabble.com.
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
