[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