help: R GUI front-end has stopped working

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]]