Alexandre, Olhei o algoritmo que vc enviou e, em momento algum, vi vc chamando "n.parc" ou "n.parc[k]", a não ser nos limites do contador k. Pelo que entendi, vc quer atribuir o valor 2, 4, 8, etc. a uma variável que irá determinar o número de parcelas. Em que ponto você usa essa variável? Seria em "at2<-at/k"? S -----Mensagem Original----- From: ASANTOS Sent: Thursday, April 12, 2012 10:41 AM To: r-br@listas.c3sl.ufpr.br Subject: Re: [R-br] Problema em loop Obrigado pela dica Sandro, Mas não funcionou não. Alexandre Em 12/04/2012 10:31, Sandro Sperandei escreveu:
Alexandre,
Notei que seu número de parcelas é sempre uma potência de 2. Correto?
Se for esse o caso, atribua á variável "número de parcelas" simplesmente o valor de 2^k.
Funcionou?
S
-----Mensagem Original----- From: ASANTOS Sent: Thursday, April 12, 2012 10:20 AM To: r-br@listas.c3sl.ufpr.br Subject: [R-br] Problema em loop
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 _______________________________________________ R-br mailing list R-br@listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.