====== Mônica Antunes Ulysséa ====== Mestranda do Programa de Pós-graduação em Zoologia da Universidade Estadual de Feira de Santana, UEFS. Laboratório de Hymenoptera do MZUSP. Tese: //Estrutura de guildas em Formigas associadas à serrapilheira de Caatingas arbóreas//. Orientador: Carlos Roberto Ferreira Brandão. ====== Exercícios ====== {{:bie5782:01_curso_atual:alunos:trabalho_final:monica:exercicio1.r|exercicio1F}} {{:bie5782:01_curso_atual:alunos:trabalho_final:monica:exercicio2final.r|exercicio2F}} {{:bie5782:01_curso_atual:alunos:trabalho_final:monica:exercico3f.r|exercicio3F}} {{:bie5782:01_curso_atual:alunos:trabalho_final:monica:exercico4.r|}} {{:bie5782:01_curso_atual:alunos:trabalho_final:monica:exercicio5f.r|}} {{:bie5782:01_curso_atual:alunos:trabalho_final:monica:exercicio6.r|exercicio6P}} {{:bie5782:01_curso_atual:alunos:trabalho_final:monica:exercicio7f.r|exercicio7F}} {{:bie5782:01_curso_atual:alunos:trabalho_final:monica:exercicio8f.r|exercicio8F}} {{:bie5782:01_curso_atual:alunos:trabalho_final:monica:exercicio9f.r|exercicio9F}} ====== Propostas de trabalho final ====== == A == **Alterada pós-comentário**\\ Criar uma função que:\\ - faça a análise de regressão linear simples entre 12 medidas morfológicas tomadas de X espécimes de formigas (nº que varia de 1 a 10 de acordo com a espécie) e a medida de tamanho total do corpo da espécie. - retorne um objeto com o sumário (para sabermos se a regressão é significativa) de todas as análises de regressão. == B == Criar uma função que calcule a média de cada uma da 12 medidas morfológicas tomadas de X espécimes (nº que varia de 1 a 10 de acordo com a espécie) para cada uma das 101 espécies de formigas coletadas e colocar o resultado em uma tabela. == Comentários == A proposta A é boa, e pode servir para qualquer conjunto de dados de medidas. Tente generalizar a função para que mais pessoas tenham interesse nela. Defina o seu output (o tipo de gráfico), e os passos para chegar nele. Sua proposta B está muito simples. ** Gabriel ** ====== Trabalho final ====== == Código da função == analise <- function(x) { if(is.matrix(x)==FALSE) { xmatrix=as.matrix(x) } else(is.matrix(x)==TRUE) { xmatrix=x } par(mfrow=c(3, 4)) for(i in 1:12) { plot(xmatrix[,i]~xmatrix[,13], ylab=colnames(dados[i]), xlab=colnames(dados[13]), main="RLS", cex.main=1) abline(lm(xmatrix[,i]~xmatrix[,13]), col="red") } sumario=vector("list", length=12) for(i in 1:12) { sumario[[i]]=summary(lm(xmatrix[,i]~xmatrix[,13])) } par(mfrow=c(1, 1)) return(sumario) } == Código comentado da função == {{:bie5782:01_curso_atual:alunos:trabalho_final:monica:trabalhofinalcodigocomentado.r|}} == Arquivo da função == {{:bie5782:01_curso_atual:alunos:trabalho_final:monica:trabalhofinalscript.r|}} == Página de ajuda == analise package:nenhum R Documentation Analise de regressão linear simples entre 12 variaveis e uma variável fixa. Description: Retorna os gráficos das relações entre variáveis com o plote da linha da análise de regressão linear simples e também o sumário das análises de regressões feitas para podermos observar o p valor. Usage: analise <- function(x) Arguments: x: objeto da classe data.frame ou matrix com valores de distintas variáveis. Details: Seu objeto deve ser construído a partir de uma tabela de dados que apresente as variáveis em colunas e as observações de cada espécime nas linhas. São usadas 13 variáveis e a última delas é a variável fixa. Ou seja, cada uma das 12 variáveis será avaliada com relaçao a variável fixa por meio da análise de regressão linear simples. Value: São gerados um conjunto de 12 gráficos para a análise exploratória de regressão linear simples entre as 12 variaveis e a variável fixa, e um objeto da classe lista que contém o sumario das regressões, onde podemos observar a significância das mesmas. Author(s): Mônica Antunes Ulysséa moniqueta_poliqueta@yahoo.com.br References: Gotelli, N.J. and Ellison, A.M. 2004. A Primer of Ecological Statistics. Sinauer Associates, Inc., Sunderland, MA, 510p. Fox, J. 2002. An R and S-Plus Companion to Applied Regression. Sage Publications, Thousand Oaks, CA, USA, xvi+312p. See Also: lm(), plot(), par(), is(), as() Examples: dados <- read.table("Dados.txt", header=TRUE, as.is=TRUE, sep="\t")\\ dados\\ analise(dados)\\ == Dados para rodar o exemplo == {{:bie5782:01_curso_atual:alunos:trabalho_final:monica:dados.txt|}}