Script das funções utilizadas no R para construir os gráficos do cap. 1 (Janzen-Connell)
Script com a nova função (Ale) para construir os gráficos do cap. 1 antigos
Dados utilizados para construir os gráficos do cap. 1 (Janzen-Connell)
read.table(“copia_de_dadosagosto2008.txt”, header=T, sep=“\t”, dec=“,”, as.is=T)→ma
## comando para ler os dados no formato que Marcia colocou no Wiki
str(ma)
'data.frame': 54 obs. of 18 variables: $ bloco : chr "A5" "A5" "A5" "A5" ... $ grupo : chr "agrupado" "agrupado" "agrupado" "agrupado" ... $ herb : chr "Ab" "Fc" "Ab" "Fc" ... $ dist : int 3 3 10 10 20 20 3 3 10 10 ... $ cap.cam : num 25.7 38.9 29.7 28.6 26.4 ... $ IH1 : num 0.6611 0.4318 0.5707 0.0655 0.2885 ... $ IH2 : num 3.13 3.26 3.23 1.68 2.27 ... $ difIH : num 2.47 2.83 2.66 1.61 1.98 ... $ mortabr : num 0.0526 0.0800 0.0263 0.0235 0.0143 ... $ comp1 : num 15.5 16.4 15.1 14.3 17.0 ... $ comp2 : num 18.0 19.0 17.6 18.0 20.2 ... $ dif.comp : num 2.47 2.55 2.55 3.66 3.25 ... $ entr1 : num 3.73 3.55 4.06 3.93 3.77 ... $ entr2 : num 5.46 5.03 5.87 5.62 5.83 ... $ dif.entr : num 1.72 1.48 1.81 1.69 2.06 ... $ nfolh1 : num 3.80 4.18 3.71 4.57 3.38 ... $ nfolh2 : num 5.11 5.28 5.01 6.17 6.17 ... $ dif.nfolh: num 1.31 1.09 1.30 1.60 2.78 ...
ma$herb[ma$herb==“Ab”]←“aberto”
ma$herb[ma$herb==“Fc”]←“fechado”
## trocando o fator “Ab” e “Fc” por “aberto” e “fechado” respectivamente para funcionar como legenda do gáfico
source(“ale.coplot.R”)
## deve carregar a função mais recente no R, chamando pelo menu <File> , e <source R code>; se a extensão do arquivo for R ele aparecera como opção caso esteja no diretório certo.
args(ale.coplot)
## mostra os argumentos definidos na função
function (x = dist, y = difIH, dados = ma, fatores = c("herb",
"grupo"), xylabel = c("Distância (m)", "Índice de herbivoria"))
## ou seja, caso não defina nada a função roda esses argumentos
ale.coplot()
##gráficos com argumentos padrão
NOTE QUE NÃO HA NECESSIDADE DE MUDAR A FUNÇÃO OU CARREGA-LA NOVAMENTE caso já tenha feito isso no workspace. Para saber se a função está carregada é só ver se aparece no ls()
ale.coplot(y=dif.entr,xylabel=c(“Distância (m)”, “dif entr”))
## DEFINE COMO Y A COLUNA dif.entr e legenda como “dif entr”, o resto é mantido no padrão
ale.coplot(y=dif.comp,xylabel=c(“Distância (m)”, “ crescimento altura”))
## DEFINE COMO Y A COLUNA dif.comp e legenda como “crescimento altura”, o resto é mantido no padrão
AGORA É COM VC