[R-br] Problema em loop

[ASANTOS]

Bom dia pessoal,

      Estou tentando dividir uma área em 2, 4, 8, 16, 32, 64, 128, 256, 
512, 1024, e 2048, para tira o CV entre as parcelas e armazenar o 
resultado em write.table(cv2, file ="area1.mod", row.names=F,append=T). 
Bom, estou conseguindo calcular um a um, mais eu gostaria de otimizar  e 
fazer um loop em um vetor que especifica o número de parcelas, como em:

n.parc<-c(2,4,8,16,32,64,128,256,512,1024,2048)##Número de parcelas

for(k in 1:length(n.parc))

mas não sei porque razão não estou conseguindo fazer esse loop, já 
tentei com n.parc e n.parc[k] e nada, segue a minha rotina completa:

####Divisão em 2, 4, 8, 16, 32, 64, 128, 256, 512, 1024, 2048 parcelas

n.parc<-c(2,4,8,16,32,64,128,256,512,1024,2048)##Número de parcelas

for(k in 1:length(n.parc)){

at2<-at/k
parc2<-round(sqrt(at2))
N_S=seq(p0[1], by=parc2, length.out=10)##Define o centro das parcelas
L_O=seq(p0[2],by=parc2, length.out=10)
grid0 <- expand.grid(N_S=seq(p0[1], by=parc2, length.out=2),
L_O=seq(p0[2],by=parc2, length.out=2))#transformar em área essas medidas 
do centro de parcela
grid=GridTopology(c(304404,  6468774), c(parc2,parc2), c(10,10))
grid=SpatialGrid(grid=grid)
plot(grid)
points(ninhos)
lines(limite)
##Passando para objetos da classe Spatial
#Pontos
pontos <- SpatialPoints(cbind(ninhos$x,ninhos$y))
pontos<- SpatialPoints(list(ninhos$x,ninhos$y))
pontos<- SpatialPoints(data.frame(ninhos$x,ninhos$y))

##Limite da área
borda=SpatialPolygons(list(Polygons(list(Polygon(limite)),"limite")))
##Plotando
plot(borda)#borda
points(pontos,col='grey')#Pontos
lines(grid,col='blue')#Malha
#remover o que esta fora do limite
fullgrid(grid)=FALSE
clip<-overlay(grid,borda)##Jogando o limite sobre a malha
sel.grid <-grid[!is.na(clip),]##Retirando os NA's que correspondem a 
malha fora do limite
points(sel.grid, col='blue')#ver o resultado final grafico
class(sel.grid)  ## SpatialPixels
ap3 <- as.SpatialPolygons.SpatialPixels(sel.grid)
table(over(pontos, ap3)) ##Frequencia em cada parcela
z<-table(over(pontos, ap3)) #contar quantos pontos tem em cada parcela
z<-as.data.frame(z) #formato somente com colunas validas
cc<- SpatialPoints(cbind(grid0$N_S,grid0$L_O))#
cc<- SpatialPoints(list(grid0$N_S,grid0$L_O))
cc<- SpatialPoints(data.frame(grid0$N_S,grid0$L_O))
clip2<-overlay(cc,ap3)
sel.cc<-cc[!is.na(clip2),]#seleção de pontos que nao estao dentro da area
#criar matriz nula
coords<-NULL#matriz sem valor para criar os pontos
coords<-sel.cc[1:length(sel.cc)] #numero para cada parcela, ate o final 
do comprimento do vetor
coords<-as.data.frame(coords)#transformar em data frame
Var1<-1:length(sel.cc)#criar novo objeto para comparação, referenciando 
a coordenada do centroide
coords<-cbind(coords,Var1) #Var1= numero das parcelas
myList <- list(coords, z)#tabela com qtdade dentro de coords e de z
final<-merge(coords, z, by.x='Var1', by.y='Var1', all.x=TRUE, 
all.y=TRUE)#comando para mostrar os NA para krigagem ordinal(limitação 
do software, não mostra as parcelas que tem 0 ninhos)
final<-final[,-1]#não precisa da 1a coluna
cv2<-sqrt(var(final$Freq))/mean(final$Freq)
write.table(cv2, file ="area1.mod", row.names=F,append=T)
result=read.table("area1.mod",header=F)
}
##########################################################


Alguém poderia me dar uma help,

Obrigado,

-- 
Alexandre dos Santos
Engenheiro Florestal, Dr.
Universidade Federal de Lavras
Departamento de Entomologia
Laboratório de Entomologia Florestal
Caixa Postal 3037
37200-000 - Lavras/MG
Fone: +55 (35) 9223-0304