Novembro 2012 - R-br

Re: [R-br] Fwd: Re: Fwd: Re: [Dúvida] Criando pacotes para o R.
by Pedro Rafael 06 Jul '13

06 Jul '13

Obrigado pela opinião Leonard Mendonça. [ ], Pedro Rafael Diniz Marinho. Em 17 de novembro de 2012 14:45, Pedro Rafael < pedro.rafael.marinho(a)gmail.com> escreveu: > Obrigado pela opinião Leonard Mendonça. > > [ ], > Pedro Rafael Diniz Marinho. > > > Em 17 de novembro de 2012 12:13, Leonard Mendonça de Assis [via R-br] < > ml-node+s2285057n4657102h39(a)n4.nabble.com> escreveu: > > No meu entendimento é uma nova versão >> >> []s >> Leonard de Assishttp://about.me/ldeassis >> >> Em 17/11/2012 11:32, Pedro Rafael escreveu: >> >> Benilton voce poderia me tirar uma dúvida? Já li a documentacao do R mas >> nao ficou claro para mim. Caso eu queira fazer uma nova versao do pacote >> pelo que eu entendi eu devo fazer upload do novo arquivo tar.gz informando >> a nova versao. Contudo, caso eu precise fazer pequenos ajustes, por >> exemplo, ajustes na documentacao. Como fazer isso sem mudar a versao do >> pacote? Isso pode ser feito? >> ---------- Mensagem encaminhada ---------- >> >> De: "beniltoncarvalho [via R-br]" <[hidden email]<http://user/SendEmail.jtp?type=node&node=4657102&i=0> >> > >> Data: 16/11/2012 17:50 >> Assunto: Re: [R-br] Fwd: Re: [Dúvida] Criando pacotes para o R. >> Para: "Pedro Rafael Diniz Marinho" <[hidden email]<http://user/SendEmail.jtp?type=node&node=4657102&i=1> >> > >> >> >> Aguarde 24-48h, é o tempo do ciclo dos "builds". >> On 16 Nov 2012 20:47, "Pedro Rafael" <[hidden email]<http://user/SendEmail.jtp?type=node&node=4657087&i=0>> >> wrote: >> >>> Pessoal os binarios para windows e Mac sao gerados pelo CRAN? Vi em >>> alguns pacotes que entraram nos espelhos hoje nao tem esses binarios. >>> Existe uma mensagem do lado dizendo que os pacotes sao novos. Os binarios >>> nao sao gerados logo quando o pacofe entra no CRAN? >>> ---------- Mensagem encaminhada ---------- >>> De: "beniltoncarvalho [via R-br]" <[hidden email]<http://user/SendEmail.jtp?type=node&node=4657087&i=1>> >>> >>> >>> Data: 16/11/2012 17:01 >>> Assunto: Re: [R-br] [Dúvida] Criando pacotes para o R. >>> Para: "Pedro Rafael Diniz Marinho" <[hidden email]<http://user/SendEmail.jtp?type=node&node=4657087&i=2>> >>> >>> >>> >>> Se vc rodar >>> >>> R CMD check pacote >>> >>> ou seu equivalente para o windows... o que sao os erros e o conteudo >>> dos logs? >>> >>> 2012/11/16 tiago souza marçal <[hidden email]<http://user/SendEmail.jtp?type=node&node=4657081&i=0> >>> > >>> >>> O meu path esta configurado da seguinte forma: >>> >>> C:\Perl\site\bin;C:\Perl\bin;C:\Rtools\bin;C:\MiKTeX2.9\miktex\bin;C:\R\R-2.15.2\bin\i386;C:\MinGW\bin;C:\HTMLHelpWorkshop;C:\Program >>> Files\Common Files\Microsoft >>> Shared\WindowsLive;%LINDOAPI_HOME%bin\win32;C:\Windows\system32;C:\Windows;C:\Windows\System32\Wbem;C:\Windows\System32\WindowsPowerShell\v1.0\;C:\Program >>> Files\Windows Live\Shared;C:\Program Files\Internet Explorer >>> >>> Você vê algum problema? >>> >>> Att. >>> >>> Tiago. >>> >>> ------------------------------ >>> Date: Fri, 16 Nov 2012 13:58:58 -0200 >>> >>> From: [hidden email]<http://user/SendEmail.jtp?type=node&node=4657081&i=1> >>> To: [hidden email]<http://user/SendEmail.jtp?type=node&node=4657081&i=2> >>> Subject: Re: [R-br] [Dúvida] Criando pacotes para o R. >>> >>> na minha vm windows, eu instalei primeiro o miktex, depois instalei o >>> rtools. Deve ser alguma loucura do windows (Path faltando, por exemplo) >>> >>> []s >>> Leonard de Assishttp://about.me/ldeassis >>> >>> Em 16/11/2012 13:11, tiago souza marçal escreveu: >>> >>> O Rtools eu tenho o problema e quando eu peço para gerar o PDF ( R CMD >>> Rd2pdf ) ele me retorna uma mensagem de erro falando do MikTex 1.9, >>> entretanto, eu estou utilizando o MikTex 2.9. >>> >>> Att. >>> >>> Tiago. >>> >>> ------------------------------ >>> Date: Thu, 15 Nov 2012 20:33:21 -0200 >>> From: [hidden email]<http://user/SendEmail.jtp?type=node&node=4657081&i=3> >>> To: [hidden email]<http://user/SendEmail.jtp?type=node&node=4657081&i=4> >>> Subject: Re: [R-br] [Dúvida] Criando pacotes para o R. >>> >>> RTools, errei o nome >>> >>> voce acha no CRAN >>> >>> []s >>> Leonard de Assishttp://about.me/ldeassis >>> >>> Em 15/11/2012 09:42, tiago souza marçal escreveu: >>> >>> Leonard onde eu consigo baixar esse rutil? >>> >>> Att. >>> >>> Tiago. >>> >>> ------------------------------ >>> Date: Thu, 15 Nov 2012 09:28:46 -0200 >>> From: [hidden email]<http://user/SendEmail.jtp?type=node&node=4657081&i=5> >>> To: [hidden email]<http://user/SendEmail.jtp?type=node&node=4657081&i=6> >>> Subject: Re: [R-br] [Dúvida] Criando pacotes para o R. >>> >>> No Windows tem que instalar o rutil recomendado para a versão do R que >>> estamos usando. >>> >>> Tenho uma VM windows 64 bits comigo e todos os problemas que encontrei >>> quando fuçei com pacotes foram referentes ao R não localizar todos os >>> 'bichinhos' do Rutil os quais precisava. >>> >>> []s >>> Leonard de Assishttp://about.me/ldeassis >>> >>> Em 15/11/2012 08:06, tiago souza marçal escreveu: >>> >>> Estou com um problema parecido com o seu, porém, no Windows . No meu >>> caso o check do PDF está ok mas quando eu peço para ele ser gerado (R >>> CMD Rd2pdf) uma mensagem de erro e retornada informando uma versão >>> antiga do MikTex que eu não estou usando. >>> >>> Mas no seu caso você observou mensagem de erro no processo de >>> instalação do pacote (R CMD INSTALL) ? Caso isso tenha ocorrido este >>> pode ter sido o problema, em função de algum erro na edição do arquivo.Rd. >>> Observe que foi retornado pra você um WARNING no check dos arquivos.Rd >>> com algumas mensagens de erro. >>> >>> Espero ter te ajudado. >>> >>> Att. >>> >>> Tiago. >>> >>> ------------------------------ >>> Date: Wed, 14 Nov 2012 23:27:37 -0300 >>> From: [hidden email]<http://user/SendEmail.jtp?type=node&node=4657081&i=7> >>> To: [hidden email]<http://user/SendEmail.jtp?type=node&node=4657081&i=8> >>> Subject: Re: [R-br] [Dúvida] Criando pacotes para o R. >>> >>> Pessoal estou estudando como criar pacotes no R e tive o seguinte erro >>> na hora de dar o check. Consigo gerar o arquivo tar.gz mas o pdf da >>> documentação não é criado. Alguém saberia me informar onde estaria o erro? >>> >>> OBS: Estou em uma máquina Linux. >>> >>> * using log directory >>> ‘/home/pedro/Documentos/Gama_Log_Logistica/FUNCAO_CODIGO/PRDM.Rcheck’ >>> * using R version 2.15.1 (2012-06-22) >>> * using platform: x86_64-pc-linux-gnu (64-bit) >>> * using session charset: UTF-8 >>> * checking for file ‘PRDM/DESCRIPTION’ ... OK >>> * checking extension type ... Package >>> * this is package ‘PRDM’ version ‘1.0’ >>> * checking package namespace information ... OK >>> * checking package dependencies ... OK >>> * checking if this is a source package ... OK >>> * checking if there is a namespace ... OK >>> * checking for executable files ... OK >>> * checking whether package ‘PRDM’ can be installed ... OK >>> * checking installed package size ... OK >>> * checking package directory ... OK >>> * checking for portable file names ... OK >>> * checking for sufficient/correct file permissions ... OK >>> * checking DESCRIPTION meta-information ... OK >>> * checking top-level files ... OK >>> * checking index information ... OK >>> * checking package subdirectories ... OK >>> * checking R files for non-ASCII characters ... OK >>> * checking R files for syntax errors ... OK >>> * checking whether the package can be loaded ... OK >>> * checking whether the package can be loaded with stated dependencies >>> ... OK >>> * checking whether the package can be unloaded cleanly ... OK >>> * checking whether the namespace can be loaded with stated dependencies >>> ... OK >>> * checking whether the namespace can be unloaded cleanly ... OK >>> * checking for unstated dependencies in R code ... OK >>> * checking S3 generic/method consistency ... OK >>> * checking replacement functions ... OK >>> * checking foreign function calls ... OK >>> * checking R code for possible problems ... OK >>> * checking Rd files ... WARNING >>> prepare_Rd: PRDM-package.Rd:30: All text must be in a section >>> prepare_Rd: PRDM-package.Rd:31: All text must be in a section >>> prepare_Rd: prdm.Rd:36-42: Dropping empty section \value >>> * checking Rd metadata ... OK >>> * checking Rd cross-references ... OK >>> * checking for missing documentation entries ... OK >>> * checking for code/documentation mismatches ... OK >>> * checking Rd \usage sections ... OK >>> * checking Rd contents ... OK >>> * checking for unstated dependencies in examples ... OK >>> * checking examples ... OK >>> * checking PDF version of manual ... ERROR >>> Rd conversion errors: >>> Warning: unknown option '-' >>> Warning in file(con, "r") : >>> cannot open file '/home/pedro/Documentos/MESTRADO': No such file or >>> directory >>> Error in file(con, "r") : cannot open the connection >>> >>> >>> [ ], >>> Pedro Rafael Diniz Marinho. >>> >>> >>> Em 14 de novembro de 2012 20:43, Pedro Rafael <[hidden email]<http://user/SendEmail.jtp?type=node&node=4657081&i=9> >>> > escreveu: >>> >>> Grande Paulo Justiniano, como vai? Muito obrigado pelas suas >>> informações. Uso o Linux faz uns 6 anos mas nunca contribui ou publiquei >>> meus códigos. Estava vendo que a linguagem Julia vem com a licença MIT que >>> eu não acho muito interessante porque quem quiser pode vender o código. >>> Minha dúvida sobre a licença GPL ou GPL2 é se uma pessoa poderá desmanchar >>> algum código que fiz ou mudar completamente o foco do pacote, retirar nome >>> dos autores, etc. Existe um pacote hoje muito usado no R que foi feito por >>> uma pessoa e depois de um tempo parou de ser mantido. Chegou um outro cara >>> de fora e reimplementou o pacote e tirou o nome dos autores antigos e >>> depois muita luta o nome dessas pessoas permaneceu como "autores antigos". >>> Isso ocorre com certa frequência? Existe alguma periodicidade de manutenção >>> do pacote que deve ser respeitada? >>> >>> [ ], >>> Pedro Rafael Diniz Marinho. >>> >>> >>> Em 14 de novembro de 2012 20:28, Paulo Justiniano [via R-br] <[hidden >>> email] <http://user/SendEmail.jtp?type=node&node=4657081&i=10>>escreveu: >>> >>> >>> >>> On Wed, 14 Nov 2012, Pedro Rafael wrote: >>> >>> > Pessoal, digamos que faço um pacote para o R. Como funciona o professo >>> de submissão? >>> >>> 1. Leia atentamente o manual "Writing R extensions" que descreve isto >>> em >>> detalhes. >>> >>> 2 certifique-se que voce testou o seu pacote : >>> >>> R CMD check SEU_PACOTE >>> >>> e corrija todos os eventuais problemas >>> >>> > Digamos que eu submeta um pacote "meu_pacote" e depois queira >>> > fazer modificações no mesmo. Para acrescentar essas modificações eu >>> tenho que submeter novamente o pacote com o mesmo nome para que as >>> atualizações >>> > fiquem disponíveis? >>> >>> Sim. Para isto voce alter a numero da versao no arquivo DESCRIPTION >>> >>> > >>> > Outra dúvida que tenho é sobre a segurança que tenho que meu pacote >>> não seja modificado. Qualquer pessoa pode pegar esse pacote e alterar ele >>> > completamente sem nenhuma permissão minha ou qualquer modificação tem >>> que ser acompanhada por intermédio da pessoa que submeteu o pacote pela >>> primeira >>> > vez? >>> > >>> >>> Bem vindo ao Universo do Software livre >>> >>> Note que se voce quiser voce pode baixar os fontes do R e mudar o que >>> quiser >>> >>> Falando mais seriamente no DESCRIPTION voce define a licenca. >>> Leia sobre licxancas GPL e análogas >>> >>> Minha sugestao: internalize completamente o conceito de software livre e >>> livre-se destas preocupacoes >>> >>> > >>> > São dúvidas bestas mas que nunca entendi direito. >>> >>> Nao sao somente suas. >>> Não é por acaso que há uma diversidade de tipos de licensas, não sem >>> controvérsias! >>> >>> > >>> > Obrigado à todos. >>> > [ ], >>> > Pedro Rafael Diniz Marinho. >>> > >>> > >>> > >>> _______________________________________________ >>> R-br mailing list >>> [hidden email] <http://user/SendEmail.jtp?type=node&node=4657009&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-Criando-pacotes-para-o-R-tp46… >>> To unsubscribe from R-br, click here. >>> NAML<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_view…> >>> >>> >>> >>> >>> _______________________________________________ R-br mailing list [hidden >>> email] <http://user/SendEmail.jtp?type=node&node=4657081&i=11> >>> 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. >>> >>> >>> _______________________________________________ >>> R-br mailing list[hidden email] <http://user/SendEmail.jtp?type=node&node=4657081&i=12>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. >>> >>> >>> >>> _______________________________________________ R-br mailing list [hidden >>> email] <http://user/SendEmail.jtp?type=node&node=4657081&i=13> >>> 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. >>> >>> >>> _______________________________________________ >>> R-br mailing list[hidden email] <http://user/SendEmail.jtp?type=node&node=4657081&i=14>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. >>> >>> >>> >>> _______________________________________________ R-br mailing list [hidden >>> email] <http://user/SendEmail.jtp?type=node&node=4657081&i=15> >>> 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. >>> >>> >>> _______________________________________________ >>> R-br mailing list[hidden email] <http://user/SendEmail.jtp?type=node&node=4657081&i=16>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. >>> >>> >>> >>> _______________________________________________ R-br mailing list [hidden >>> email] <http://user/SendEmail.jtp?type=node&node=4657081&i=17> >>> 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. >>> >>> _______________________________________________ >>> R-br mailing list >>> [hidden email] <http://user/SendEmail.jtp?type=node&node=4657081&i=18> >>> 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. >>> >> >> >> _______________________________________________ >> R-br mailing list >> [hidden email] <http://user/SendEmail.jtp?type=node&node=4657081&i=19> >> 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-Criando-pacotes-para-o-R-tp46… >> To unsubscribe from R-br, click here. >> NAML<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_view…> >> >> _______________________________________________ >> R-br mailing list >> [hidden email] <http://user/SendEmail.jtp?type=node&node=4657087&i=3> >> 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. >> >> _______________________________________________ >> R-br mailing list >> [hidden email] <http://user/SendEmail.jtp?type=node&node=4657087&i=4> >> 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-Criando-pacotes-para-o-R-tp46… >> To unsubscribe from R-br, click here. >> NAML<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_view…> >> >> >> _______________________________________________ >> R-br mailing list[hidden email] <http://user/SendEmail.jtp?type=node&node=4657102&i=2>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. >> >> >> >> _______________________________________________ >> R-br mailing list >> [hidden email] <http://user/SendEmail.jtp?type=node&node=4657102&i=3> >> 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-Criando-pacotes-para-o-R-tp46… >> To unsubscribe from R-br, click here<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=unsubscrib…> >> . >> NAML<http://r-br.2285057.n4.nabble.com/template/NamlServlet.jtp?macro=macro_view…> >> > >

5 18

Imputação múltipla: como fazer operações com os dados imputados.
by Lucas Petri Damiani 12 Dez '12

12 Dez '12

Boa noite, sou leitor assíduo da lista, embora nunca tenha participado das discussões. Porém, hoje vou me arriscar e passar um problema pra vcs. Espero que possam me ajudar. Estou realizando uma análise num banco de dados de cerca de 10 mil pacientes de Trauma. O objetivo da pesquisadora é testar a validade de um escore de gravidade (preditor de sobrevida) em pacientes que entram no pronto socorro e eventualmente promover alterações nesse escore. Porém, algumas medidas que ela deseja testar possuem cerca de 20% de missing. E a exclusão desses pacientes faz com que o novo escore predito fique um bocado diferente do original, o que me fez ir atrás de alternativas para tratar esses dados faltantes e assim me arrisquei nesse pacote *mice. *Em testes mais simples que fiz, fiquei bastante empolgado com o método. Enfim, tentei reproduzir abaixo o problema que encontrei. ##### Banco de dados de exemplo: dados <- matrix(c( 4, 6, 5, 34, 120, 5, 1, 0, NA, 1, 1, 3, 80, 6, 0, 1, 2, 4, 2, 24, 90, 5, 1, 1, 3, 5, 5, 5, 98, 5, 1, 1, NA, 6, 5, 40, 94, 4, 0, 0 , NA, 6, 5, 12, 87, 3, 0, 1, 1, 5, NA, 32, 70, 2, 0, 1, 1, 5, 1, 7, 60, 3, 1, 0, 2, NA, 2, 25, 80, 5, 1, 0, 4, 6, 5, 24, 101, 7, 1, 0, 1, 5, 1, 1, 120, 8, 1, 0, 3, 5, 3, 12, 200, 6, 0, 1, 3, 3, 1, 3, 40, 5, 1, 0, NA, 5, 1, 35, 60, 4, 0, 1, 1, 4, 1, 32, 87, 3, 1, 0, 4, 5, 1, 23, 30, 2, 1, 1, 2, 1, 1, 5, 46, 1, 0,1), byrow=T, ncol=8) dados <- as.data.frame(dados) colnames(dados) <- c("ECGL_AO", "ECGL_MRM", "ECGL_MRV", "FR", "PAS", "GRAU", "IOT", "SOBREVIDA") dados$IOT <- as.factor(dados$IOT) dados$ECGL <- apply(dados[, c("ECGL_AO", "ECGL_MRV", "ECGL_MRM")],1,sum) dados$ECGL_cut <- cut(dados$ECGL, breaks=c(3,4,6,9,13,16), right=FALSE) # Notem que a medida RTS é um escore derivado de outras medidas, que pretendo imputar. dados$RTS <- (0.9368*(dados$ECGL_cut=="[4,6)") + 2*0.9368*(dados$ECGL_cut=="[6,9)") + 3*0.9368*(dados$ECGL_cut=="[9,13)") + 4*0.9368*(dados$ECGL_cut=="[13,16)") + 0.7326*(dados$PAS>=1 & dados$PAS<=49) + 2*0.7326*(dados$PAS>=50 & dados$PAS<=75) + 3*0.7326*(dados$PAS>=76 & dados$PAS<=89) + 4*0.7326*(dados$PAS>=90) + 0.2908*(dados$FR>=1 & dados$FR<=5) + 2*0.2908*(dados$FR>=6 & dados$FR<=9) + 4*0.2908*(dados$FR>=10 & dados$FR<=29) + 3*0.2908*(dados$FR>=30)) #Imputando dados nas variáveis "ECLG": require(mice) PRED <- (1 - diag(1, ncol(dados[,1:8]))) PRED[8,1:8] <- PRED[1:8,8] <- 0 teste <- mice(dados[,1:8], method=c("pmm","pmm", "pmm", "", "", "", "", ""), seed=15, predictorMatrix = PRED) #O meu problema começa agora. Eu não sei como calcular medidas derivadas no objeto teste. E a fórmula abaixo retorna erro, já conferi várias vezes muito dos parênteses e etc, mas creio que o problema seja que a fórmula é muito grande entro da função I(x), pois quando eu reduzo o modelo, eu consigo resultados satisfatórios. Me corrijam se estiver errado, por favor. fit_error <- with(teste, glm(SOBREVIDA ~ I(0.9368*((ECGL_AO + ECGL_MRV + ECGL_MRM) >= 4) & ((ECGL_AO + ECGL_MRV + ECGL_MRM) <= 5) + 2*0.9368*(((ECGL_AO + ECGL_MRV + ECGL_MRM) >= 6) & ((ECGL_AO + ECGL_MRV + ECGL_MRM) <= 8)) + 3*0.9368*I(((ECGL_AO + ECGL_MRV + ECGL_MRM) >= 9) & ((ECGL_AO + ECGL_MRV + ECGL_MRM) <= 12)) + 4*0.9368*I(((ECGL_AO + ECGL_MRV + ECGL_MRM) >= 13) & ((ECGL_AO + ECGL_MRV + ECGL_MRM) <= 15)) + 0.7326*I(PAS>=1 & PAS<=49) + 2*0.7326*I(PAS>=50 & PAS<=75) + 3*0.7326*I(PAS>=76 & PAS<=89) + 4*0.7326*I(PAS>=90) + 0.2908*I(FR>=1 & FR<=5) + 2*0.2908*I(FR>=6 & FR<=9) + 4*0.2908*I(FR>=10 & FR<=29) + 3*0.2908*I(FR>=30)) + IOT, family="binomial")) #Enfim, gostaria de rodar algo como isso escrito abaixo, mas não sei como definir essa variável RTS dentro do arquivo teste. fit <- with(teste, glm(SOBREVIDA ~ RTS + IOT, family="binomial")) fit2 <- glm(SOBREVIDA ~ RTS + IOT, family="binomial", data=dados) pool(fit) # Uma saída que imagino é gerar os cinco bancos imputados através da função "complete(teste, m)" e ajustar a regressão. Mas dessa forma não posso utilizar a função *pool. *Teria que reprogramá-la. De qualquer forma, estou perdido. Qualquer ideia é bem vinda. Abraços, Lucas Petri Damiani

2 2

Distribuição Gamma Incompleta
by Leonardo Monteiro 05 Dez '12

05 Dez '12

Boa tarde, to trabalhando com séries de dados diários de precipitação da ANA (Hidroweb) e gostaria de saber se alguem sabe um pct ou script para distr. gamma incompleta, com o objetivo de preenchimento dessas falhas? obrigado, -- Leonardo Amaral Monteiro - Eng° Agrônomo Mestre em "Física do Ambiente Agrícola" - Agrometeorologia Doutorando do PPG em "Engenharia de Sistemas Agrícolas" Escola Superior de Agricultura "Luiz de Queiroz" - Universidade de São Paulo Piracicaba,SP. Fone: (19)34294123 - Ramal:236

3 2

Áreas correlatas - Hospitais são acusados de substituir médicos por enfermeiros em cirurgias
by Richard Santos 03 Dez '12

03 Dez '12

Se profissionais de áreas correlatas podem exercer as profissões correlatas, esse caso não precisaria ter punição, certo? Porque ciências humanas inclusive é mais intuitivo do que ciências exatas. Quanto à quantidade de conhecimento estatística é tão robusta quanto medicina. E ambas causam mortes quando feito errado. É essa opinião que eu não entendo, para mim, profissão é coisa muito séria. Abraços, Richard http://www.correiobraziliense.com.br/app/noticia/cidades/2012/11/27/interna…

5 8

Cotação Histórica de Ações Bovespa
by T Branquinho 01 Dez '12

01 Dez '12

Caros Colegas da Lista, Alguém já conseguiu importar as cotações históricas das Ações da Bovespa pelo R? Se sim, como o fez? http://www.bmfbovespa.com.br/shared/iframe.aspx?idioma=pt-br&url=http://www… No mais agradeço, Tales

6 6

leitura de dados
by Bernardo Rodrigues 01 Dez '12

01 Dez '12

Boa tarde. gostaria de saber como leios os dados de um documento em r. coloquei dados<-read.table("dados.r",header=F) como faço para os dados .txt e nao deu. Obrigada

5 5

Instalação do Revolution R - pede serial number
by Richard Santos 30 Nov '12

30 Nov '12

Caros, A instalação do Revolution R pede um serial number, que parece ser de 12 dígitos. Qual digito? rsrssrsr Richard

6 12

Construção de Mapa
by Edson Lira 30 Nov '12

30 Nov '12

Bom dia a todos, tenho um banco de dados de doação apta de 2000 a 2012 (dados até outubro). Abaixo os links dos arquivos que estou trabalhando. Sou novo na plotagm de mapas, estou querendo aprender. Estou usando o material do Jakson (por sinal, muito bom). http://www.datafilehost.com/download-f19b102f.html# script http://www.datafilehost.com/download-62c89b81.html #dados http://www.datafilehost.com/download-88dbbc7c.html # shape bairros zonas http://www.datafilehost.com/download-88dbbc7c.html #prj http://www.datafilehost.com/download-090c879a.html# sbn http://www.datafilehost.com/download-c7d90384.html #sbx Gostaria de fazer dois mapas com esses dados. 1. Plotar o mapa de manaus com as zonas geográficas, tendo no mapas os nomes(labels) das zonas e na legenda as doações aptas (quantidade de doações aptas). 2.Fazer a mesma rotina, só que por bairros, neste caso seriam necessários duas legendas(uma para bairro, pois acho que os nomes não caberiam, e a outra para a quantide de doação por bairro). Espero ter sido claro. [ ]'s. Edson Lira Estatístico Manaus-Amazonas

4 4

dúvida
by alanarocha＠sapo.pt 30 Nov '12

30 Nov '12

estive a ver o help(abline) mas o unico exemplo que encontro é para o modelo y~x onde posso encontrar um exemplo da função abline() para lm(pthi~dose+grupo) ?Ana

5 4

modelo e gráfico
by alanarocha＠sapo.pt 30 Nov '12

30 Nov '12

Presados, tenho este modelo como os dois níveis como sei o valor estimado do nível default criado pelo R? modCH<-lm(pthi~dose+grupo) modCH summary(modCH) outra questão é que sei que abline() é para modeos do tipo y~x eu queria um exemplo de como fazer as rectas para o meu modelo que envolve factor desde já obrigada Ana

3 2