Agosto 2011 - R-br

Re: [R-br] Correlação Canônica!!
by Thiago De paula protásio 27 Ago '11

27 Ago '11

Obrigado Gustavo; Dei uma olhada no pacote vegan, mas acho que ele é mais apropriado para dados ecológicos.... Obrigado! Thiago de Paula Protásio Acadêmico de Engenharia Florestal Universidade Federal de Lavras Laboratório de Energia da Biomassa Florestal Laboratório de Biomateriais (035) 9183-2246

2 1

Re: [R-br] Salvar Base de Dados com +400.000 linhas
by Roney Fraga Souza 26 Ago '11

26 Ago '11

Daniel, Fiz o procedimento que vc deseja com a PNAD 2009. Não obtive problema, o arquivo .csv ficou com 42 MB e o PSPP (alternativa livre ao SPSS) importou normalmente. Meu computador é um mac, c2d (64bits) com 2gb de ram e utilizo o R no Eclipse. > load('PNAD2009.RData') load('PNAD2009.RData') > ls() ls() [1] "PNAD2009" > dim(PNAD2009) dim(PNAD2009) [1] 399387 383 > PNAD2009 <- PNAD2009[,1:40] PNAD2009 <- PNAD2009[,1:40] write.csv(PNAD2009, file='pnad2009.csv', quote=FALSE) > write.csv(PNAD2009, file='pnad2009.csv', quote=FALSE) Att Roney

1 0

Interpretação saída do pacote contrast
by ASANTOS 26 Ago '11

26 Ago '11

Bom dia pessoal, Bom já sei que esse assunto já foi motivo de tópicos passados no nosso mail-list e no internacional ,tem material na internet etc, mais essa saída do pacote contraste não entrou na minha cabeça ainda não e venho pedir ajuda do grupo. Estou comparando minha variável resposta número de óvulos de fêmeas de um parasitóide com as variáveis explicativa tratamento (com 2 níveis: 1 casal e 10 casais) e geração (com 4 níveis: 3, 6, 9 e 12 gerações). Na saída da funcao contrast(), ele solta o trata- mento10 casais no intercepto(imagino) e não faz as comparacoes dentro de cada geração, ele só desdobra para o tratamento 1casal, mesmo assim sumiu com a geracao 12, onde pergunto onde que estou errando é na funcao ou na interpretação? Sendo os argumentos que usei: > ##Contraste tratamentos x desenvolvimento oviposicao a mumia > m.comple<-glm(des.ov.mu~tratamento*geracao,data=dados6, family="quasipoisson") > anova.comple<-anova(m.comple,test="Chi") > anova.comple Analysis of Deviance Table Model: quasipoisson, link: log Response: des.ov.mu Terms added sequentially (first to last) Df Deviance Resid. Df Resid. Dev P(>|Chi|) NULL 82 9.1366 tratamento 1 1.0969 81 8.0397 2.265e-07 *** geracao 3 2.2880 78 5.7518 4.445e-12 *** tratamento:geracao 3 2.7715 75 2.9802 1.328e-14 *** --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 > #summary(m.comple) 10 casais > # > ##Fazendo os contrastes > require(contrast) > contr<-contrast(m.comple, list(geracao=levels(geracao),tratamento="1casal"),list(geracao=levels(geracao),tratamento="10casal")) > print(contr,X=TRUE) glm model parameter contrast Contrast S.E. Lower Upper t df Pr(>|t|) -0.298869262 0.03278233 -0.36417501 -0.23356352 -9.12 75 0.0000 0.006271251 0.03351287 -0.06048982 0.07303232 0.19 75 0.8521 -0.103435587 0.03532313 -0.17380286 -0.03306831 -2.93 75 0.0045 0.049914391 0.03237891 -0.01458771 0.11441650 1.54 75 0.1274 Contrast coefficients: (Intercept) tratamento1casal geracaog3 geracaog6 geracaog9 0 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 1 0 0 0 tratamento1casal:geracaog3 tratamento1casal:geracaog6 tratamento1casal:geracaog9 0 0 0 1 0 0 0 1 0 0 0 1 > ## > summary(m.comple)## Dando uma olhada no modelo Call: glm(formula = des.ov.mu ~ tratamento * geracao, family = "quasipoisson", data = dados6) Deviance Residuals: Min 1Q Median 3Q Max -0.38017 -0.09964 0.00000 0.03447 0.87147 Coefficients: Estimate Std. Error t value Pr(>|t|) (Intercept) 2.09063 0.02250 92.934< 2e-16 *** tratamento1casal -0.29887 0.03278 -9.117 8.93e-14 *** geracaog3 -0.15385 0.03229 -4.765 9.04e-06 *** geracaog6 -0.15621 0.03313 -4.715 1.09e-05 *** geracaog9 -0.05985 0.03230 -1.853 0.067820 . tratamento1casal:geracaog3 0.30514 0.04688 6.509 7.62e-09 *** tratamento1casal:geracaog6 0.19543 0.04819 4.055 0.000121 *** tratamento1casal:geracaog9 0.34878 0.04608 7.570 7.83e-11 *** --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Null deviance: 9.1366 on 82 degrees of freedom Residual deviance: 2.9802 on 75 degrees of freedom AIC: NA Number of Fisher Scoring iterations: 4 Obrigado, Alexandre

2 3

Salvar Base de Dados com +400.000 linhas
by Daniel Dantas 26 Ago '11

26 Ago '11

Boa tarde a cada um de vcs! Em que formato eu salvo uma base de dados que tem mais de 400.000 linhas e no max 30-40 colunas? A ideia é salvar essa base de dados para futuramente abrir no SPSS (seria o melhor) ou qualquer outro software com uma interface mais prática que a do R, pois preciso encaminhá-la para uma pessoa que não entende nada de R. Na minha incessante procura, o que encontrei foram pacotes que somente importavam e não exportavam formato .sav. Se puderem me ajudar com esse problema, e até mesmo comentar de outros formatos e softwares interessantes que possam ser utilizados para divulgação desses resultados obtidos com o R para usuários de softwares que são comumente usados por não estatísticos-matemáticos. Obrigado, Daniel

7 11

OFF-topic[artigo]
by Pedro Emmanuel Alvarenga Americano do Brasil 26 Ago '11

26 Ago '11

Amigos de R, alguem tem acesso a esse artigo completo, que poderia pegar pra mim? http://www.jstor.org/pss/2530910 , Abraço forte e que a força esteja com você, Dr. Pedro Emmanuel A. A. do Brasil Instituto de Pesquisa Clínica Evandro Chagas Fundação Oswaldo Cruz Rio de Janeiro - Brasil Av. Brasil 4365 Tel 55 21 3865-9648 email: pedro.brasil(a)ipec.fiocruz.br email: emmanuel.brasil(a)gmail.com ---Apoio aos softwares livres www.zotero.org - gerenciamento de referências bibliográficas. www.broffice.org ou www.libreoffice.org - textos, planilhas ou apresentações. www.epidata.dk - entrada de dados. www.r-project.org - análise de dados. www.ubuntu.com - sistema operacional

2 2

Leitura de vários arquivos
by Leandro Marino 26 Ago '11

26 Ago '11

*Caros, estou há algum tempo pesquisando no histórico da lista e nao encontro. Estou com uma rotina que diariamente tenho que importar vários arquivos. Existe alguma forma de importar todas para data.frames diferentes cada um com seu nome original?! Ex.: Arq1.txt ar1312.txt lll.txt Arq1 <- read.table('Arq1.txt') ar1312 <- **read.table('ar1312.txt') . . . Como os recebo via ftp e cada dia com um nome mais criativo que o outro fica complicado. Já uso o dir() para listar e salvar os nomes. Assim vou num editor de texto e manipulo esses nomes, entretanto, gostaria de saber se existiria uma forma de importalos mais rápido... ** *Atenciosamente, Leandro Marino http://www.leandromarino.com.br (Fotógrafo) http://est.leandromarino.com.br/Blog (Estatístico) Cel.: + 55 21 9845-7707 Cel.: + 55 21 8777-7907

5 6

Reorganização de banco de dados
by Greice Laureano 26 Ago '11

26 Ago '11

Bom dia, Gostaria de ajuda no seguinte problema: Tenho a um conjunto de dados assim: ID AF1MEDIDA1 AF1MEDIDA2 AF2MEDIDA1 AF2MEDIDA2 AF3MEDIDA1 AF3MEDIDA2 1 11.7 5.0 21.7 21.0 33.0 20.3 E gostaria de transformar ele nesse tipo: ID REP AF VALOR 1 MEDIDA1 AF1 11.7 1 MEDIDA2 AF1 5 1 MEDIDA1 AF2 21 1 MEDIDA2 AF2 21 1 MEDIDA1 AF3 33 1 MEDIDA2 AF3 20.03 Obrigado, Att, Greice.

4 9

Re: [R-br] [Dúvida] Invocar funções C++ no R
by FHRB Toledo 26 Ago '11

26 Ago '11

Me perdoem por enviar arquivos... Espero compreensão de todos! att, FH 2011/8/26 FHRB Toledo <fernandohtoledo(a)gmail.com> > Siga esse tutorial p/ invocar funções em C no R p/ windows... tudo > funcionando! > > att, > FH > > 2011/8/26 Pedro Rafael <pedro.rafael.marinho(a)gmail.com> > >> Já vi o erro. Fiz uma besteira de compilar o código em 32bits e tava >> tentando rodar no R 64bits. Na verdade eu acho que um código 32bits era para >> rodar sem problemas em no 64bits. >> >> Em 26 de agosto de 2011 09:40, Pedro Rafael < >> pedro.rafael.marinho(a)gmail.com> escreveu: >> >> Benilton realmente no Windows o arquivo gerado e uma dll como você tinha >>> falado. Consegui gerar um programa chamado hello.c. O conteúdo do programa >>> apresenta-se logo abaixo: >>> >>> #include <R.h> >>> void hello(int *n) >>> { >>> int i; >>> for(i=0; i < *n; i++) { >>> Rprintf("Hello, world!\n"); >>> } >>> } >>> >>> Fui no terminal do DOS e digitei R CMD SHLIB hello.c e o arquivo >>> hello.dll foi gerado. Setei o R no diretório da dll e fiz: >>> >>> > dyn.load("hello.dll") >>> Erro em inDL(x, as.logical(local), as.logical(now), ...) : >>> impossível carregar objeto compartilhado 'C:/PEDRO/PACOTES_R/hello.dll': >>> LoadLibrary failure: %1 não é um aplicativo Win32 válido. >>> specificado. >>> >>> Como você pode ver deu um erro. O que pode ter sido? >>> >>> Em 26 de agosto de 2011 09:19, beniltoncarvalho [via R-br] < >>> ml-node+3770707-1188459073-223914(a)n4.nabble.com> escreveu: >>> >>> Confirme que sua variavel PATH aponta para a localizacao correta do >>>> RTools e MinGW. Por exemplo, se o RTools foi instalado em "c:\RTools", >>>> entao confirme que a variavel PATH tenha os seguintes valores >>>> inclusos: >>>> >>>> c:\Rtools\bin;c:\Rtools\MinGW\bin; >>>> >>>> Eu so' nao tenho certeza se, depois de mudar os valores de PATH, vc >>>> vai precisar reiniciar o computador... >>>> >>>> Uma vez que voce tenha certeza de que a variavel PATH tem os valores >>>> corretos, voce abre um novo terminal DOS e tenta o grep.... >>>> >>>> b >>>> _______________________________________________ >>>> R-br mailing list >>>> [hidden email] <http://user/SendEmail.jtp?type=node&node=3770707&i=0> >>>> 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. >>>> >>>> >>>> ------------------------------ >>>> If you reply to this email, your message will be added to the >>>> discussion below: >>>> >>>> http://r-br.2285057.n4.nabble.com/R-br-Duvida-Invocar-funcoes-C-no-R-tp3768… >>>> To unsubscribe from R-br, click here<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscrib…>. >>>> >>>> >>> >>> >>> >>> -- >>> Saudações, >>> Pedro Rafael Diniz Marinho. >>> >>> >>> >> >> >> -- >> Saudações, >> Pedro Rafael Diniz Marinho. >> >> >> >> _______________________________________________ >> 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. >> > >

1 0

[Dúvida] Invocar funções C++ no R
by Pedro Rafael 26 Ago '11

26 Ago '11

Como se procede para importar uma função C++ ou até mesmo C para o R? -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6

2 15

Identação do Código R
by Pedro Rafael 26 Ago '11

26 Ago '11

Pessoal, uma dúvida besta. Existe um padrão para identação de um código R? Várias linguagens usam um padrão específico como Python por exemplo. -- Saudações, Pedro Rafael Diniz Marinho - Estatístico SES-PB Currículo Lattes: http://buscatextual.cnpq.br/buscatextual/visualizacv.jsp?id=K4250792T6

6 8