====== Maíra Tir ====== {{:bie5782:01_curso_atual:alunos:trabalho_final:maira:dsc06483.jpg?200|}} Mestranda em Ecologia no Instituto de Biociências da USP com a tese "Polimorfismo genético e taxonomia convencional de cianobactérias da represa Billings", orientada pelo Prof. Dr. Marcelo Luis Martins Pompêo. Estatística nunca foi o meu forte, mas nunca desanimei com os meus fracassos. ^nome^projeto^ email^telefone^ |Maíra Paula Tir Serico|Polimorfismo genético e taxonomia convencional de cianobactérias da represa Billings| mairatir@usp.br|82216056| ===== Meus Exercícios ===== [[.:exec]] Exercício 1: {{:bie5782:01_curso_atual:alunos:trabalho_final:maira:exercício1.r|F}} Exercício 2: {{:bie5782:01_curso_atual:alunos:trabalho_final:maira:exercício2.r|F}} Exercício 3: {{:bie5782:01_curso_atual:alunos:trabalho_final:maira:exercício3.r|F}} Exercício 4: {{:bie5782:01_curso_atual:alunos:trabalho_final:maira:exercício4.r|F}} Exercício 5: {{:bie5782:01_curso_atual:alunos:trabalho_final:maira:exercício5.r|F}} Exercício 6: {{:bie5782:01_curso_atual:alunos:trabalho_final:maira:exercício6.r|F}} Exercício 7: {{:bie5782:01_curso_atual:alunos:trabalho_final:maira:exercício7.r|P}} Exercício 8: Exercício 9: {{:bie5782:01_curso_atual:alunos:trabalho_final:maira:exercício9.r|P}} ====== Trabalho Final====== Fiz a proposta 2 ====== HELP ====== biovolume package:- R Documentation Calcula o biovolume. Descrição: Realiza o cálculo do biovolume de células com diferentes formatos. Foi criado inicicialmente para células de organismos fitoplanctônicos, mas pode ser usado para outras situações. Utiliza médias de cada medida necesária para o cálculo do formato gheométrico. Uso: biovolume(formato, comprimento, largura, altura, diametro, n=1) Argumentos: formato: Corresponde ao formato da célulae, portanto, à equação que será utilizada. As possibilidades de formato são: "sphere", "prolate spheroid", "ellipsoid", "cylinder", "cylinder + 2 half spheres", "cylinder + 2 half cones", "cone", "double cone", "cone + half sphere", "box", "prism on elliptic-base", "prism on parallelogram-base", "sickle-shaped prism", "sickle-shaped cylinder", "pyramid". Demais formatos podem ser feitos através de combinações dos contemplados pela função. comprimento: Valores correspondentes ao comprimento da célula. largura: Valores correspondentes à largura da célula. altura: Valores correspondentes à altura da célula. diametro: Valores correspondentes ao diâmetro da célula. n: Valor correspondente ao número de células do indivíduo medido ou da colônia a qual ele pertence. Detalhes: O cálculo do biovolume depende da especificação do formato (forma geométrica) da célula, não sendo possível a sua execução sem essa definição. Para cada formato há uma equação adequada que pode utilizar diferentes medidas. As medidas necessárias para cada formato podem ser consultadas na própria função. Valor: A função retorna primeiro exemplos de organismos fitoplanctônicos correspondentes à forma geométrica utilizada para os cálculos, em seguida aponta quais medidas são necessárias para que o cálculo possa ser feito e, por fim, mostra o resultado do cálculo do biovolume caso as medidads tenham sido inseridas. Cuidado: É preciso estar atento à correta escolha da forma geométrica que será utilizada. A função biovolume calcula o biovolume de apenas um indivíduo (consederando que os indivíduos são filamentos, colônias e organismos unicelulares), portanto, para obter o biovolume total da amostra deve-se saber a densidade do indivíduo na amostra e multiplicá-la pelo biovolume. Autor: Maíra Paula Tir Serico Referências: Sun, J. & Liu, D. (2003) Geometric models for calculating cell biovolume and surface area for phytoplankton. Journal of Plankton Research. Exemplos: c <- c(1,2,3,4,5) #comprimento de 5 indivíduos l <- c(6,7,8,9,10) #largura de 5 indivíduos a <- c(11,12,13,14,15) #alturade 5 indivíduos d <- c(16,17,18,19,20) #diâmetro de 5 indivíduos biovolume(formato="sphere", diametro=d, n=5) #cálculo do biovolume para 5 indivíduos "sphere" biovolume(formato="sickle-shaped prism", comprimento=c, largura=l, altura=a, n =5) #cálculo do biovolume para 5 indivíduos com formato prisma em forma de foice biovolume(formato="sphere") #Consulta sobre o que deve ser medido para conseguir calcular o biovolume para o formato "sphere" ====== Código ====== biovolume <- function(formato="", comprimento="", largura="", altura="", diametro="", n=1) { #Condição para cada formato if(formato=="sphere") { volume <- (((pi*mean(diametro)^3)/6)*mean(n)) algas <- data.frame(Cyanobacteria=c("Anabaena", "Chroococcus", "Microcystis"), Chlorophyceae=c("Coelastrum", "Chlorella", "Dictiosphaerium")) medidas <- c("Medir o diâmetro da célula e conte o número (n) de células") } if(formato=="prolate spheroid") { volume <- ((pi*mean(diametro)^2*mean(altura))/6)*mean(n) algas <- data.frame(Cyanobacteria=c("Aphanothece", "Xenococcus"), Chlorophyceae=c("Oocystis","Scenedesmus"), Cryptophyceae=c("Cryptomonas","Hillea"), Chrysophyceae=c("Mallonmonas","Dinobryon"), Dinophyceae=c("Balechina","Ptychodiscus")) medidas <- c("Medir o diâmetro e a altura e conte o número (n) de células") } if(formato=="ellipsoid") { volume <- ((pi*mean(comprimento)*mean(largura)*mean(altura))/6)*mean(n) algas <- data.frame(Dinophyceae=c("Peridinium", "Amphidinium", "Gyrodinium"), Euglenophyceae=c("Trachelomonas")) medidas <- c("Medir comprimento, largura e altura da célula e conte o número (n) de células") } if(formato=="cylinder") { volume <- ((pi*mean(diametro)^2*mean(comprimento))/4)*mean(n) algas <- data.frame(Cyanobacteria=c("Cylindrospermopsis", "Oscillatoria"), Bacillariophyceae=c("Aulacoseira", "Cyclotella")) medidas <- c("Medir o diâmetro e o comprimento da célula e conte o número (n) de células") } if(formato=="cylinder + 2 half spheres") { volume <- (pi*mean(diametro)^2*(mean(altura)/4-mean(diametro)/12))*mean(n) algas <- data.frame(Bacillariophyceae=c("Chrysanthemodiscus", "Skeletonema")) medidas <- c("Medir diâmetro e altura da célula e conte o número (n) de células") } if(formato=="cylinder + 2 half cones") { volume <- (pi/4*mean(diametro)^2*(mean(altura)-mean(diametro)/3))*mean(n) algas <- data.frame(Chlorophyceae=c("Actinastrum", "Ankistrodesmus", "Closteriopsis")) medidas <- c("Medir diâmetro e altura da célula e conte o número (n) de células") } if(formato=="cone") { volume <- ((pi*mean(altura)*mean(diametro)^2)/12)*mean(n) algas <- data.frame(Chrysophyceae=c("Pyramidochrysis", "Sphaleromantis")) medidas <- c("Medir diâmetro e altura da célula e conte o número (n) de células") } if(formato=="double cone") { volume <- ((pi*mean(altura)*mean(diametro)^2)/12)*mean(n) algas <- data.frame(Dinophyceae=c("Schuettiella"), Chlorophyceae=c("Brachiomonas")) medidas <- c("Medir diâmetro e altura da célula e conte o número (n) de células") } if(formato=="cone + half sphere") { volume <- ((pi*mean(altura)*mean(diametro)^2)/4)*mean(n) algas <- data.frame(Chrysophyceae=c("Ochromonas"), Cryptophyceae=c("Chroomonas")) medidas <- c("Medir diâmetro e altura da célula e conte o número (n) de células") } if(formato=="box") { volume <- (mean(comprimento)*mean(largura)*mean(altura))*mean(n) algas <- data.frame(Cyanobacteria=c("Merismopedia"), Bacillariophyceae=c("Asterionella", "Bacillaria")) medidas <- c("Medir o comprimento, largura e altura da célula e conte o número (n) de células") } if(formato=="elliptic-base prism") { volume <- ((pi*mean(comprimento)*mean(largura)*mean(altura))/4)*mean(n) algas <- data.frame(Bacillariophyceae=c("Achnanthidium", "Fragilaria"), Chlorophyceae=c("Pediastrum"), Euglenophyceae=c("Phacus")) medidas <- c("Medir comprimento, largura e altura e conte o número (n) de células") } if(formato=="prism on parallelogram-base") { volume <- (1/2*mean(comprimento)*mean(largura)*mean(altura))*mean(n) algas <- data.frame(Bacillariophyceae=c("Gyrosigma", "Cymatoneis", "Nitzschia", "Nitzschiella")) medidas <- c("Medir comprimento, largura e altura da célula e conte o número (n) de células") } if(formato=="sickle-shaped cylinder") { volume <- ((pi*mean(altura)*mean(diametro)^2)/6)*mean(n) algas <- data.frame(Chlorophyceae=c("Ankistrodesmus", "Monorraphidium", "Closteriopsis", "Kirchineriella"), Bacillariophyceae=c("Licmophora")) medidas <- c("Medir o diâmetro e altura da célula e conte o número (n) de células") } if(formato=="sickle-shaped prism") { volume <- ((pi*mean(comprimento)*mean(largura)*mean(altura))/4)*mean(n) algas <- data.frame(Bacillariophyceae=c("Eunotia")) medidas <- c("Medir comprimento, largura e altura da célula e conte o número (n) de células") } if(formato=="pyramid") { volume <- (((mean(diametro)^2)*mean(altura))/6)*mean(n) algas <- data.frame(Dynophyceae=c("Tetradinium")) medidas <- c("Medir diâmetro (diagonal da base) e altura da célula e conte o número (n) de células") } resulta <- list(algas,medidas,volume) names(resulta) <- c("Gêneros","O que deve ser medido/contado?","Biovolume") return(resulta) } ====== Proposta 1 ====== A idéia é criar uma função que produza um gráfico de mais de uma variável em relação a profundidade de um corpo d'água. Seria útil para mim porque esse tipo de gráfico é muito ilustrativo para apresentações e interpretações de dados limnológicos e é muito difícil de encontrar programas que o façam (eu ainda não encontrei =/). As varáveis principais envolvidas nessas análises são oxigênio dissolvido, condutividade elétrica, ph, temperatura, nutrientes, entre outros... Se em um gráfico eu pudesse juntar algumas dessas variáveis seria muito útil e acredito que não só para mim. **Comentários** É factível, mas acho que poderia ter mais algo, como por exemplo diagnóses das associações entre as variáveis, talvez. --- //[[macfabio@gmail.com|Fabio de A. Machado]] 2011/04/06 19:16// ====== Proposta 2 ====== É uma proposta mais complicada e ainda não pensei muito bem como seria feita, mas gostaria de ter uma função que facilitasse o cálculo de biovolume de espécies de fitoplâncton. Atualmente faço esses cálculos com o Excel (UUUURRGGHHH!!!), mas não está sendo muito prático e gostaria de otimizar esse processo. **Comentários** Minha impressão geral é que se dá para fazer no Excel, dá para fazer no R de forma mil vezes mais simples. Eu diria para você ir no plano B que parece mais interessante, ainda mais que você vai ter uma função que você já tem demanda. Fora isso, o Excel erra em alguns contextos muito específicos. São dois pontos positivos para sua proposta 2. --- //[[macfabio@gmail.com|Fabio de A. Machado]] 2011/04/06 19:30//