[R-br] Problema com representação de parcelas usando expand.grid()

ASANTOS alexandresantosbr em yahoo.com.br
Quinta Setembro 4 23:21:00 BRT 2014


Boa noite Pessoal,

         Estou tentando representar as parcelas de um experimento em 
campo e estou tendo dois problemas sendo:

1) _Representação do centro das parcelas_:

## Coordenadas do grid no espaçamento das plantas de 2x4m

## distâncias verticais e horizontais entre plantas
v <- 2; h <- 4

## número de plantas vertical e horizontal
npv <- 4; nph <- 5

## distância entre vãos vertical e horizontal (distância entre parcelas)
V <- 7; H <- 7

## número de parcelas em  linhas e colunas
ncol <- 4; nrow <- 12

grid0 <- expand.grid(L=seq(2,length.out=nrow,by=npv*v+V),
                      C=seq(2,length.out=ncol,by=nph*h+H),
                      x=seq(2,length.out=npv,by=v),
                      y=seq(2,length.out=nph,by=h))
grid1 <- with(grid0, data.frame(x=L+x, y=C+y))

plot(y~x, grid1, asp=1) ### Representação das plantas


     Até aí tudo bem, porém quando tento representar por pontos o centro 
das parcelas, começa o problema, sendo:

### O centro de cada parcela vai se localizar a cada 13 m na horizontal 
e a cada 23 na vertical

grid2 <- expand.grid(seq(7,length.out=13*4,by=13),
                      seq(12,length.out=23*12,by=23))

points(grid2[,1],grid2[,2],col="red") ## Vejo que os pontos caem errado 
e não sei por que?


2) _Criar uma coluna para identificar a parcela a qual pertence cada 
planta_:


pos<-1:48
text(grid2[,1],grid2[,2]+0.5,labels=pos, cex=1,col=2)

e depois que o centro cair no local certo, gostaria de fazer grid1$pos 
para identificar a qual centro da parcela cada ponto pertence.


        Espero que alguém possa dar um help,

Obrigado,


-- 
======================================================================
Alexandre dos Santos
Proteção Florestal
IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso
Campus Cáceres
Caixa Postal 244
Avenida dos Ramires, s/n
Bairro: Distrito Industrial
Cáceres - MT                      CEP: 78.200-000
Fone: (+55) 65 8132-8112 (TIM)   (+55) 65 9686-6970 (VIVO)
e-mails:alexandresantosbr em yahoo.com.br
         alexandre.santos em cas.ifmt.edu.br
Lattes: http://lattes.cnpq.br/1360403201088680
======================================================================

-------------- Próxima Parte ----------
Um anexo em HTML foi limpo...
URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140904/7412f21d/attachment.html>


Mais detalhes sobre a lista de discussão R-br