Dezembro 2013 - R-br

Coords - GeoR
by João Sousa 05 Dez '13

05 Dez '13

Gostaria de ajuda/orientação de como obter os vetores das coordenadas dos bairros da cidade do Rio de Janeiro. O objetivo é montar um gráfico espacial para mostrar o número de furtos, roubos e homicídios dolosos na cidade. Tive um curso desse na graduação, há muuuito tempo e a coisa simplesmente caiu no esquecimento. []' João Sousa

3 8

Parcela subdividida: problema na declaração do erro associado a parcela
by Jose Claudio Faria 04 Dez '13

04 Dez '13

Pessoal, O CMR abaixo é da apostila do Walmes. #------------------------------------------------------------------------------------------ # dados ps <- expand.grid(BL=c("I", "II", "III", "IV"), ES=c("e1", "e2"), AD=c("a1", "a2", "a3")) ps$alt <- c(58,77,38,52,44,59,30,34,85,90,73,77,59,68,45,55,66,93,67,64,54,75,53,48) str(ps) # #------------------------------------------------------------------------------------------ # análise de variância (erro A = BL x AD) m0 <- aov(alt ~ BL + AD + Error(BL:AD) + ES + AD:ES, data=ps) # termo Error para declarar erro A m0 <- aov(alt ~ BL + AD * ES + Error(BL:AD), data=ps) # o mesmo com fórmula comprimida summary(m0) #------------------------------------------------------------------------------------------ Encontrei muitos outros exemplos (na web e em livros) com essa notação: Error(BL:AD) Contudo o R (3.0.2patched) não está aceitando essa notação! Fornece sempre a mensagem de erro: In aov(alt ~ BL + AD + Error(BL:AD) + ES + AD:ES, data = ps) : modelo Error() é singular Contudo o R está aceitando esse modelo: m0 <- aov(alt ~ BL + AD + Error(BL/AD) + ES + AD:ES, data=ps) Mas nesse caso o resíduo A (associado a bloco e parcela - AD) não fica correto. Alguém que anda trabalhando com esses modelos poderia ajudar? Ab, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica UESC/DCET/Brasil joseclaudio.faria at gmail.com Telefones: 55(73)3680.5545 - UESC 55(73)9100.7351 - TIM 55(73)8817.6159 - OI ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

4 7

Plotar média de dois grupos
by Andre Oliveira 04 Dez '13

04 Dez '13

Pessoal, dei uma pesquisada e não consegui plotar o que desejo. Estou com os comandos. require(gdata) dados<-read.xls("multinacional.xls", encoding = "latin1") str(dados) attach(dados) library(doBy) require(lattice) xyplot(salario~escolaridade|sexo) summaryBy(salario~escolaridade|sexo,dados,FUN=c(mean)) escolaridade sexo salario.mean 1 médio Feminino 12.80000 2 médio Masculino 10.77500 3 pós-graduação Feminino 32.03333 4 pós-graduação Masculino 33.28000 5 superior Feminino 21.67000 6 superior Masculino 25.13478 Como plotar em um mesmo gráfico os gráficos de barras das médias dos três escolaridade que tenho por sexo um com cada cor? Tem como fazer isto de forma automática ou vou ter que criar os dados em um data frame para plotar? obrigado

3 4

Re: [R-br] anova com dois fatores não paramétrica
by Antonio Moita 04 Dez '13

04 Dez '13

Shah, D. A., and Madden, L. V. 2004. Nonparametric analysis of ordinal data in designed factorial experiments. Phytopathology 94:33-43. Existe um artigo complementar que a revista chama de extra *The e-Xtra logo stands for “electronic extra” and indicates that the online article contains supplemental material not included in the print edition. The online supplement contains detailed instructions on use of statistical software and interpretation of the output, SAS programs (with data), and annotated output files. Segundo os autores foi implementada um pacote no R. Att, Moita -------------------------------------------- Em ter, 3/12/13, Luiz Henrique da Conceicao Leal <richfield(a)ig.com.br> escreveu: Assunto: [R-br] anova com dois fatores não paramétrica Para: r-br(a)listas.c3sl.ufpr.br Data: Terça-feira, 3 de Dezembro de 2013, 11:32 Alguém sabe qual é alternativa não paramétrica quando não se pode assumir normalidade no caso de análise de variância com dois fatores? Desde já agradeço Luiz -----Anexo incorporado----- _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.

2 1

Repositório do CRAN no Nordeste do Brasil
by Jose Claudio Faria 04 Dez '13

04 Dez '13

Olá para todos, O novo repositório experimental do CRAN no Nordeste do Brasil (UESC/Bahia) está no ar. O dados estão abaixo: - URL: http://nbcgib.uesc.br/mirros/cran/ - Country: Brazil - Name: Brazil (BA) - Cidade: Ilhéus - Host: NBCGIB - Núcleo de Biologia Computacional e Gestão de Informações Biotecnológicas. Para que possa ser incluido na lista dos repositórios oficiais do CRAN solicitamos testes. Ab, -- ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria Estatistica UESC/DCET/Brasil joseclaudio.faria at gmail.com Telefones: 55(73)3680.5545 - UESC 55(73)9100.7351 - TIM 55(73)8817.6159 - OI ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\

1 0

scan()
by Fernanda De Bastiani 04 Dez '13

04 Dez '13

Caros, alguém conhece alguma função no R que possa substituir a função scan()? No sentido de que a função scan() precisa de um file para fazer scan, ou então voce digitar todos os valores (segundo o que eu entendo). Por exemplo: tenho um object de classe list, e gostaria de ler este objeto, sem precisar salvá-lo em file primeiro, ou sem precisar digitá-lo. n = c(2, 3, 5) s = c("aa", "bb", "cc", "dd", "ee") b = c(TRUE, FALSE, TRUE, FALSE, FALSE) object = list(n, s, b, 3) # x contains copies of n, s, b eu gostaria poder fazer scan(object), existe algo semelhante que funcione?! Obrigada, Fernanda

5 6

Re: [R-br] Gráfico de pontos
by Alex 04 Dez '13

04 Dez '13

Tiago, obrigado por se oferecer a ajudar. Informo que encontrei a solução e compartilho com os demais integrantes da lista r-br. Ressalto que o gráfico exibe pontos com cores e tamanhos segundo o quartil do dado observado. Abraços, Alex Santos quantile(dadosT1$data) points(dadosT1,xlab="Coordenada X (coluna)",ylab="Coordenada Y (linha)", pt.divide="quartiles") legend(270,70, c("1º quartil", "2º quartil", "3º quartil", "4º quartil"), col=c('black', 'black', 'black', 'black'), pt.bg=c('blue', 'green', 'yellow', 'red'), pch=c(21,21,21,21), pt.cex=c(0.6,0.9,1.3,1.6),bty='0') Dec 04, 2013; 8:21am [image: Tiago de Souza]<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=user_nodes…> Tiago de Souza<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=user_nodes…> Nos envie um CMR para que possamos tentar ajudá-lo. Att. Tiago. ######################################################## Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas ######################################################## Em 4 de dezembro de 2013 08:12, Alex <alxcart(a)gmail.com> escreveu: > Ao grupo r-br, > por favor, eu criei um gráfico de pontos onde o tamanho e cor de cada um é > baseado no quartil ao qual ele pertence, o trecho do código está abaixo. > Entretanto estou com dificuldades de criar a legenda. Alguém poderia me > orientar? > > > *attach(dadosT1)* > *par(mfrow=c(1,1))* > *ramp<-colorRamp(c("yellow", "red"))* > *points(dadosT1,xlab="Coordenada X (coluna)",* > *ylab="Coordenada Y (linha)",* > *pt.divide="quartiles", col=rgb(ramp(seq(0,1, length=5)), max=255))* > > > > Desde já agradeço, > Alex Santos > > >

1 0

Obter médias a cada 5 colunas em um data frame
by Rodrigo Muniz 04 Dez '13

04 Dez '13

Caros membros da lista, Gostaria pedir um ajuda bastante pontual, tenho um banco de dados que possui 1850 colunas que representam o comprimento de ondas de alguns tratamento, o comprimento de ondas está apresentado de 1 a 1 nanômetro, sendo que não há diferenças se for apresentado de 5 em 5, por isso gostaria de fazer uma médias a cada 5 colunas, para reduzir o volume de dados. Já tentei transpor a matriz de dados e utilizar a função tapply, porém eu ficaria com um outro data frame com 720 colunas, o que não ajudaria muito. Alguém poderia me falar alguma função para que realize essa operação? Desde já agradeço. ATT. -- *Rodrigo A. Muniz* Eng. Agrônomo. Ms Produção Vegetal (UENF) Doutorando em Engenharia de Sistemas Agrícolas (ESALQ/USP) Cel (19) 8300-4333 (Pessoal) Tel (19) 3375-1196 (Residêncial) Tel (19) 3447-8537 (Trabalho)

3 2

Algoritmo de Fisher's Scoring para o procedimento REML
by Tiago Souza Marçal 04 Dez '13

04 Dez '13

Bom dia pessoal, Estou trabalhando com a estimação de componentes de variância via modelos mistos usando o procedimento de máxima verossimilhança restrita (REML), no entanto, em minhas pesquisas pela internet encontrei o algoritmo de Fisher's Scoring para solução do seguinte modelo: y = Xa + Zg + e em que: y é o vetor de dados a e g são vetores de efeito fixo e aleatório respectivamente. X e Z são as matrizes de incidência dos efeitos fixo e aleatório respectivamente. Gostaria de poder adicionar mais efeitos aleatórios no modelo como por exemplo y = Xa + Zg + Wp + e, e encontrar as soluções dos componentes de variância pelo algoritmo de Fisher's Scoring. Alguém me indica algum material onde possa encontrar algo a respeito? Att. Tiago. ################################################################# Tiago de Souza Marçal - Graduando em Agronomia pelo CCA-UFES Bolsista de Iniciação Científica da área de Genética e Melhoramento de Plantas #################################################################

1 0

Plotar intervalo de confiança
by Fátima Lima Paula 04 Dez '13

04 Dez '13

Prezados, tenho a taxa de mortalidade por idade e queria plotar o intervalo de confiança. Alguém pode me ajudar? pop=c(37374, 54783, 86049, 126048, 287384) obfaF=c(6,15,27,55,239) txobfaF=(obfaF*365/popfaF)*100 txobfaF idade=c("1","2","3","4","5") plot(idade, txobfaF,type="b",main="Sexo Feminino",xlab="Idade",ylab="Taxa de Mortalidade",lty=2,ylim=c(1,40),xaxt='n') Obrigada Fátima -- "Minha felicidade depende da qualidade dos meus pensamentos"

6 14