Hello,
I'm using the function nlminb of the package stats inside a loop and when the
number of trials grows, R crashes and says "R
GUI front-end has stopped working". Could you help me with this problem? I have
try in versions 2.15.1,2.15.2 and 3.0.0.
> sessionInfo()
R version 2.15.2 (2012-10-26)
Platform: i386-w64-mingw32/i386 (32-bit)
locale:
[1] LC_COLLATE=Galician_Spain.1252 LC_CTYPE=Galician_Spain.1252
[3] LC_MONETARY=Galician_Spain.1252 LC_NUMERIC=C
[5] LC_TIME=Galician_Spain.1252
attached base packages:
[1] stats graphics grDevices utils datasets methods base
This is my code:
u<-runif(3000)
k=13
n=length(u)
prob=vector(length=n)
ro=vector(length=n)
beta1=vector(length=n)
g1=vector(length=n)
sb1=vector(length=n)
low=vector(length=n)
tlow=vector(length=n)
effects=vector(length=n)
vb1=vector(length=n)
AA=list()
BB=list()
nnI=list()
library(stats)
for (i in 1:n){
v=as.numeric(u<=u[i])
p=mean(v)
A=vector(length=k)
n1=n%/%k
A<-sapply(1:k,function(m) A<-sum(v[((m-1)*n1+1):(m*n1)]))
A[k]=sum(v[((k-1)*n1+1):n])
AA[[i]]<-A
B=c(rep(n1,k-1),length(v[((k-1)*n1+1):n]))
BB[[i]]<-B
l1=vector(length=(length(A)))
L1<-function(pe,rho){
for(j in 1:length(A)){
if(A[j]==0){l1[j]=0} else {l1[j]<-sum(log(pe+(-rho/(rho-1))*(0:(A[j]-1))))}
}
return(sum(l1))}
l2=vector(length=(length(B-A)))
L2<-function(pe,rho){
for(j in 1:length(B-A)){
if((B[j]-A[j])==0){l2[j]=0} else
{l2[j]<-sum(log(1-pe+(-rho/(rho-1))*(0:(B[j]-A[j]-1))))}
}
return(sum(l2))}
l3=vector(length=(length(B)))
L3<-function(pe,rho){
for(j in 1:length(A)){
if(B[j]==0){l1[j]=0} else {l3[j]<-sum(log(1+(-rho/(rho-1))*(0:(B[j]-1))))}
}
return(sum(l3))}
L<-function(pe,rho){
L<-L1(pe,rho)+L2(pe,rho)-L3(pe,rho)
return(L)
}
Max<- function(x){ -L(x[1], x[2])}
opt<-nlminb(c(0.01,0.01), Max,lower = rep(0.001,2), upper =
rep(0.999,2),control=list(rel.tol=1e-6))
prob[i]=opt$par[1]
ro[i]=opt$par[2]
}
Thanks,
Irene Castro Conde.
[[alternative HTML version deleted]]
______________________________________________
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.
