====== Camilla Pagotto ====== {{:bie5782:01_curso2009:alunos:trabalho_final:fulano.jpg?200|Fulano de Tal é um cara legal}} Mestranda em Ecologia, Instituto Biociências, USP. O título de minha tese é: "Pequenos mamíferos arborícolas em um mosaicode vegetações em diferentes estádios sucessionais: influência da estrutura da floresta e disponibilidade de recursos sobre a distribuição das espécies.", orientado pela Profª Drª Renata Pardini. ====== Trabalho Final ====== No meu projeto investigo como a regeneração florestal influência a guilda dos pequenos mamíferos arborícolas, para isso (além de outras variáveis) medidas sobre a estrutura da vegetação foram coletadas. ===== Plano A ===== Estratificação vertical da folhagem Para a quantificação da estratificação vertical da folhagem utilizei de uma vara de 3 m como mira para o estabelecimento de uma coluna imaginária de 15 cm de diâmetro, onde o limite inferior e superior das faixas de folhagem foram anotados. Essas medidas foram coletadas em 24 pontos por sítio de amostragem, um no início e outro no fim de 12 transectos de 15m. A metragem ocupada por folhagem em diferentes estratos da floresta (0-1 m, 1-5 m, 5-10 m, 10-15 m, 15-20 m, 20-25 m, e assim por diante) será calculada para cada um dos 24 pontos de amostragem de cada sítio. Quero criar uma função onde a partir desses dados de quantificação da estratificação da vegetação, para cada um dos 24 pontos, eu posso extrair o grau de desconexão da vegetação ou seja, uma medida dos intervalos sem preenchimento por folhagem. ==== Comentários PI ==== Promissora, mas não consigo avaliar completamente pq algumas coisas não ficaram claras. Pense nelas e vamos conversando: * O primeiro passo para pensar na função é definir claramente //input// e //outputs//. Não entendi muito bem o que seria esta desconexão: um valor, um vetor com vários valores? * Da mesma forma, não ficou muito claro para mim o formato de entrada de dados: uma planilha (dataframe)? O que vai ans linhas e nas colunas? * Esta análise seria feita para cada ponto, certo? Então que tipo de gráfico ou resumo exploratório poderia ser produzido para mostrar a variação pelos pontos? ==== Resposta - Tentando fazer a função ==== Paulo, a partir dos meus dados brutos de estratificação, como estes de uma vareta: 2,4-2,5 / 9,9-10 / 11-11,2 / 12-17, onde de 2,4 a 2,5 tem preenchimento de vegetação e de 2,5 a 9,9 não tem preenchimento, e assim por diante. Eu preciso do valor dos espaços sem preenchimento, então pensando em como começar a minha função, percebi que não tenho como "puxar" esses dados de uma planilha, já que tenho mais valores do que o número de linhas no excell, sendo assim, cada um dos dados de cada vareta terá que ser inserido na "unha". Como no exemplo abaixo, com intervalos de 0,01cm: a<-c(0,seq(2.4,2.5,0.01),seq(9.9,10,0.01),seq(11,11.2,0.01),seq(12,17,0.01)) b<-diff(a) # diferença entre medidas de uma única vareta. c<-b[which(b[] > 0.02)] # para mim não importa valores menores que 0,02, pq o meu menor intervalo é de 0,05. d<-sum( c ) e<-d/max(a) # preciso da soma dos valores dos espaços sem preenchimento dividido pelo maior valor (que é a altura do dossel em cada vareta) ffff<-mean(e1,e2,e3....e12) # média de todas as varetas, um valor por sítio. Com esses passos eu consigo ter uma medida de espaços sem preenchimento divididos pela altura do dossel, o que seria uma medida de desconexão, o contrário dessa medida é a minha medida de interesse, que é a conexão de vegetação por sítio. ==== Plano B ==== Para a quantificação das variáveis que possivelmente afetam a guilda dos pequenos mamíferos arborícolas (riqueza e abundância), foram amostrados 28 sítios em 3 estádios de regeneração florestal (médio, médio-avançado e avançado), onde foram coletados dados sobre disponibilidade de recursos alimentares (frutos e artrópodes de sub-bosque) e recursos espaciais (estrutura da vegetação). Para uma análise exploratória desses dados, vou criar uma função onde serão estabelecidos previamente quem são as variáveis dependentes e independentes, e dessa forma sejam feitas correlações entre elas, além dos gráficos correspondentes. === Comentários PI=== Outra boa idéia. novamente, pense com cuidado os dados de entrada e o que vc quer dar ao usuário como saída. No caso há várias maneiras de explorar correlações, e também de calculá-las. Uma cosia útil e com menos pressupostos são aquelas linhas de ajuste não paramétrico (//splines//) que mostrei na aula. Mas há várias outras maneiras de explorar relações entre variáveis. Um bom começo é o capítulo de análise exploratória do livro do Verzani. ===Página de Ajuda=== EVF package:nenhum R Documentation Função para calcular intervalos sem preenchimento de folhagem nos estratos verticais da vegetação através de um método cuja denominação, nesta função é "varetas", tendo como "input" uma lista. Description: Calcula a média entre intervalos sem preenchimento por folhagem ao longo dos estratos verticais de vegetação (modificação do método descrito por Malcolm (1995)). Calculando a média por sítio e retornando como resultado um valor de conexão da vegetação. Usage: EVF(x) Arguments: x: Vetor numérico. Valores de uma lista. Details: Através de uma lista informada pelo usuário, os vetores serão os valores das varetas, onde os intervalos com preenchimento de vegetação estarão dentro de parênteses, e dentro deste um intervalo de sequencia (0.01) informado. A função roda mediante uma lista com os vetores concatenados. Os valores de dentro dessa lista (x) terão as médias das diferenças extraídas, uma média geral e um valor inverso(multiplicado por -1), que é o resultado final. Value: A média das diferenças entre os valores de uma lista é extraída. O resultado final é a média dos vetores, multiplicado por -1. Author(s): Camilla Presente Pagotto capagotto@usp.br References: MALCOLM, J. R. Forest structure and the abundance and diversity of neotropical small mammals. In: Lowman, M. D.; Nadtorni, N. M. Forest Canopies. Academic Press. San Diego. pp 179-197. 1995. (para a descrição do método modificado de quantificação da estratificação vertical da folhagem). CRAWLEY, M. J. The R Book. Wiley, New York, 2007. See Also: 'lapply', para criar uma lista e extrair a diferença dentro dos vetores. Examples: a<-c(0,seq(2.4,2.5,0.01),seq(9.9,10,0.01),seq(11,11.2,0.01),seq(12,17,0.01)) b<-c(0,seq(1,3.75,0.01),seq(10,13,0.01),seq(13.1,15.2,0.01),seq(17,22,0.01),seq(22.5,23,0.01)) c<-c(seq(0,2.5,0.01),seq(3,4,0.01),seq(5.5,11,0.01)) dados<-list(a,b,c) dados EVF <-function(x) ===Código da Função=== EVF <-function(x) { media.vareta<-rep(1,length(x)) for(i in 1:length(x)) { difer<-diff(x[[i]]) inter<-difer[which(difer[] > 0.02)] media.vareta[i]<-mean(inter) } media.total<-mean(media.vareta) conexao<- -1*media.total return(conexao) } EVF(dados)