Outubro 2012 - R-br

Modelos segmentados de regressão não linear...
by andrebvs 25 Out '12

25 Out '12

3 11

simulacao gama
by Natalia Martins 25 Out '12

25 Out '12

Prezados redistas, boa tarde estou tentando simular uma gama, usei o comando abaixo no R, no entanto eu gostaria de por um valor maximo e um valor minimo (truncá-la), alguem pode me ajudar? Muito obrigada. radiacao <- rgamma(n, shape=3, rate=2) -- Natália da Silva Martins Bacharel em Estatística - Universidade Estadual de Maringá/ UEM Mestranda em Estatística e Experimentação Agronômica - ESALQ/ USP Contato: (19) 8306-4743

4 3

Re: [R-br] Transformação de dados no dataframe
by fabrizzioconde＠gmail.com 25 Out '12

25 Out '12

Caro Augusto, Muito obrigado pelas dicas de codificação, resolveram muito bem o problema, mas ainda estou com problema na montagem da tabela final deste e-mail, pois o comando table que sugeriu conta o número de casos de AA, BB e AB na tabela final e o que necessito é posicionar de maneira diferente o dado e não contar o número de vezes que o mesmo aparece nos fatores. AA11AB2 AA21AA0 AA31BB1 AA41BB1 AA12AA0 AA22AB2 AA32AB2 AA42AA0 AA13BB1 AA23BB1 AA33AA0 AA43AA0 AA14BB1 AA24AB2 AA34AA1 AA44AB2 Desta forma, note que o elemento (1,5) na tabela acima foi colocado na posição (1,1) da tabela abaixo. O elemento (5,5) foi colocado na posição (2 1). O elemento (9,5) foi inserido na posição (3,1) e o elemento (13,5) está na posição (4,1). Assim, consigo capturar a variável AA1 de cada indivíduo (1, 2, 3 e 4). Da mesma forma, o elemento (2,5) foi colocado na posição (1 2), o elemento (6,5) foi inserido na posição (2,2) e assim sucessivamente. AA1AA2AA3AA4 12011 20220 31100 41212 Não consigo fazer esse redimensionamento da tabela antiga para a nova. Desde já agradeço imensamente sua ajuda. Abraços, Fabrízzio -------Original Message------- From: r-br-request(a)listas.c3sl.ufpr.br Date: 24/10/2012 12:00:14 To: r-br(a)listas.c3sl.ufpr.br Subject: Digest R-br, volume 20, assunto 30 Enviar submissões para a lista de discussão R-br para r-br(a)listas.c3sl.ufpr.br Para se cadastrar ou descadastrar via WWW, visite o endereço https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br ou, via email, envie uma mensagem com a palavra 'help' no assunto ou corpo da mensagem para r-br-request(a)listas.c3sl.ufpr.br Você poderá entrar em contato com a pessoa que gerencia a lista pelo endereço r-br-owner(a)listas.c3sl.ufpr.br Quando responder, por favor edite sua linha Assunto assim ela será mais específica que "Re: Contents of R-br digest..." Tópicos de Hoje: 1. Re: Transformação de dados no dataframe (Augusto Ribas) ---------------------------------------------------------------------- Message: 1 Date: Wed, 24 Oct 2012 09:50:12 -0300 From: Augusto Ribas <ribas.aca(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br Subject: Re: [R-br] Transformação de dados no dataframe Message-ID: <CACMkfRyM6o+=mvtJ7GuQXn2xMiNYgGdo1Za0m0erZ6aw6Xc09g(a)mail.gmail.com> Content-Type: text/plain; charset="windows-1252" Vamos supor um exemplo de dados: #dados de exemplo dados<-data.frame(col1=round(runif(20,1,4)),col2=sample(c("A","B"),20 replace=T),col3=sample(c("A","B"),20,replace=T)) dados #primeiro com o comando paste, vamos juntar a coluna 2 com a 3, criando os AA AB e BB dados$col4<-factor(paste(dados$col2,dados$col3,sep="")) #Depois como os walmes nos ensinou aqui na lista, a gente usa o match() #para recodificar nos valores que queremos, no caso como vc falou levels(dados$col4) recod<-c(0,1,1,2) dados$col5<-recod[match(dados$col4,levels(dados$col4))] dados #agora que vc te todos os ingredientes, é so usar o comando table() para fazer essa tabela #eu fiz um exemplo de tabela com a coluna 5 e a coluna 5, mas pelo que vi vc quer a coluna #5 e os nomes das linhas, rownames(sua planilha), entao use os nomes da linha na tabela table(dados$col5,dados$col4) Espero que ajude. Abraços Em 24 de outubro de 2012 09:16, fabrizzioconde(a)gmail.com < fabrizzioconde(a)gmail.com> escreveu: > Bom dia, > > Existe algum script que faça o seguinte: > > 1) Cria uma quinta coluna na matriz abaixo baseada na terceira e na quarta > coluna DA seguinte forma: AA = 0, BB = 1 e AB = 2. > > AA1 1 A B AA2 1 A A AA3 1 B B AA4 1 B B AA1 2 A A AA2 2 A B AA3 2 A B > AA4 2 A A AA1 3 B B AA2 3 B B AA3 3 A A AA4 3 A A AA1 4 B B AA2 4 A B AA3 > 4 A A AA4 4 A B > > Assim, a matriz se transformaria em: > > AA1 1 A B 2 AA2 1 A A 0 AA3 1 B B 1 AA4 1 B B 1 AA1 2 A A 0 AA2 2 A B 2 > AA3 2 A B 2 AA4 2 A A 0 AA1 3 B B 1 AA2 3 B B 1 AA3 3 A A 0 AA4 3 A A 0 > AA1 4 B B 1 AA2 4 A B 2 AA3 4 A A 1 AA4 4 A B 2 > > 2) Após esse passo, tenho que eliminar as colunas 3 e 4 DA matriz acima > (talvez em um novo dataframe) e redimensioná-la para o seguinte formato: > > AA1 AA2 AA3 AA4 1 2 0 1 1 2 0 2 2 0 3 1 1 0 0 4 1 2 1 2 > > Ao final, preciso que cada coluna DA matriz acima seja encarada pelo R > como uma variável (AA1, AA2, AA3, AA4), pois irei aplicar o pacote de > regressão logística penalizada stepPlr e considerar essas variáveis > preditoras. > > Muito obrigado. > Fabrízzio > [image: FREE Christmas Animations for your email ? by IncrediMail! > Click Here!]<http://www.incredimail.com/?id=620188&did=10500&ppd=2747 201111071523,9,[TypeID],[IM_UPN2]&rui=139771782&sd=20121024> > > _______________________________________________ > 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. > -- Grato Augusto C. A. Ribas Site Pessoal: http://augustoribas.heliohost.org Lattes: http://lattes.cnpq.br/7355685961127056 -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr br/pipermail/r-br/attachments/20121024/96f8c562/attachment-0001.html> -------------- Próxima Parte ---------- Um anexo não-texto foi limpo... Nome: xmas_girl_en.gif Tipo: image/gif Tamanho: 61041 bytes Descrição: não disponível URL: <http://listas.inf.ufpr br/pipermail/r-br/attachments/20121024/96f8c562/attachment-0001.gif> ------------------------------ _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Fim da Digest R-br, volume 20, assunto 30 *****************************************

4 4

Ajuda para criar um grafico
by Humberto 24 Out '12

24 Out '12

Ola pessoal, gostaria de saber se alguem poderia me ajudar na criacao desse grafico em anexo (link para baixar) a partir da tabela abaixo: Species lci or uci 1 Bigeye tuna 0.0952 0.1293 0.176 2 Yellowfin tuna 0.0916 0.1255 0.172 3 Albacore 0.0913 0.1562 0.267 4 Swordfish 0.1311 0.1565 0.187 5 Sailfish 0.0368 0.0923 0.231 6 White marlim 0.1300 0.2161 0.359 7 Blue marlim 0.0738 0.1471 0.293 8 Blue shark 0.1303 0.1653 0.210 9 Crocodile shark 0.1123 0.1792 0.286 10 Pelagic stingray 0.1074 0.2052 0.392 11 Turtle oliva 0.1275 0.3802 1.134 12 Oceanic whitetip shark 0.1658 0.3788 0.866 desde ja agradeco Humberto Eu inseri um link para 1 arquivo nesta mensagem: figure.docx <http://www.yousendit.com/download/WUJhQ3Q1Qk50NjlMWE1UQw>(172KB)YouSendIt <https://www.yousendit.com>http://www.yousendit.com/download/WUJhQ3Q1Qk50NjlMWE1UQw O Mozilla Thunderbird <http://www.getthunderbird.com> torna fácil compartilhar grandes arquivos pelo e-mail.

1 0

Transformação de dados no dataframe
by fabrizzioconde＠gmail.com 24 Out '12

24 Out '12

Bom dia, Existe algum script que faça o seguinte: 1) Cria uma quinta coluna na matriz abaixo baseada na terceira e na quarta coluna DA seguinte forma: AA = 0, BB = 1 e AB = 2. AA11AB AA21AA AA31BB AA41BB AA12AA AA22AB AA32AB AA42AA AA13BB AA23BB AA33AA AA43AA AA14BB AA24AB AA34AA AA44AB Assim, a matriz se transformaria em: AA11AB2 AA21AA0 AA31BB1 AA41BB1 AA12AA0 AA22AB2 AA32AB2 AA42AA0 AA13BB1 AA23BB1 AA33AA0 AA43AA0 AA14BB1 AA24AB2 AA34AA1 AA44AB2 2) Após esse passo, tenho que eliminar as colunas 3 e 4 DA matriz acima (talvez em um novo dataframe) e redimensioná-la para o seguinte formato: AA1AA2AA3AA4 12011 20220 31100 41212 Ao final, preciso que cada coluna DA matriz acima seja encarada pelo R como uma variável (AA1, AA2, AA3, AA4), pois irei aplicar o pacote de regressão logística penalizada stepPlr e considerar essas variáveis preditoras. Muito obrigado. Fabrízzio

2 1

Gráficos animado no R
by Jodavid Ferreira 24 Out '12

24 Out '12

Gostaria de saber como fazer um gráfico animado no R tipo, um gráfico, que vai se transformando em outro! -- Jodavid Ferreira Responsável pela comissão de Comunicação do MJ da RCCPB - http://jovem.rccpb.com.br Membro da equipe Estadual do Ministério de Comunicação da RCCPB - http://rccpb.com.br Responsável pelo Ministério de Comunicação do Grupo de Oração Boa Esperança - http://grupoboaesperanca.com Bacharelando em Estatística - UFPB Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig Site Pessoal - http://jodavidferreira.com Twitter - @JodavidFerreira <http://twitter.com/JodavidFerreira>

3 2

Regressão Logística
by Edson Lira 23 Out '12

23 Out '12

Estou usando as rotinas abaixo, PARA LER OS DADOS h<-read.csv2("C:\\Consulta\\CIESA\\bd.csv") h[h==""]<-NA names(h) table(h$quem_elab) str(h) h$tempo<-2012-h$ano_ini h$melh_orig<-h$niv_melhora[drop=T] h$niv_melhora<-h$niv_melhora[drop=T] h$bsc<-h$bsc[drop=T] levels(h$niv_melhora) levels(h$niv_melhora)[levels(h$niv_melhora)%in% c("elevado","muito elevado")] <- "Melhora" levels(h$niv_melhora)[levels(h$niv_melhora)%in% c("nenhum","pouco", "intermediário")] <- "Sem melhora" #levels(h$niv_melhora)[levels(h$niv_melhora)%in% c("intermediário")] <-"" levels(h$quem_elab) h$part<-h$quem_elab levels(h$part)[levels(h$part)%in% c("alta direção")]<- "Direção" levels(h$part)[levels(h$part)%in% c("alta direção\ngerentes", "alta direção\ngerentes\nfuncionários","alta direção\ngerentes\noutros","gerentes", "funcionários")] <- "Gerentes" table(h$part) #PARA SELECIONAR SOMENTE AS VARIÁVEIS QUE PENSO USAR NO MODELO nb<-data.frame(melhora=h$melh_orig,bsc=h$bsc,champion=h$champion,ap_dir=h$apoio_bsc, niv_tre=h$niv_treina,cons_ext=h$cext_bsc,reu_f_est=h$reu_func_estr, com_est=h$comu_estr,ajuste=h$ajustes,ind_apr=h$ind_aprop,ind_desemp=h$ind_desd_dep, mapa_est=h$mapa_estr,tempo=h$tempo) levels(nb$melhora) levels(nb$melhora)[levels(nb$melhora)%in% c("elevado","muito elevado")] <- 1 levels(nb$melhora)[levels(nb$melhora)%in% c("nenhum","pouco", "intermediário")] <- 0 nb$desemp<-as.numeric(as.character(nb$melhora)) #PARA ELIMINAR OS NA'S cd <- subset(nb, subset=!is.na(desemp) & !is.na(bsc)& !is.na(champion)& !is.na(ap_dir)& !is.na(niv_tre)& !is.na(cons_ext) & !is.na(reu_f_est)& !is.na(com_est)& !is.na(ajuste)& !is.na(ind_apr)& !is.na(ind_desemp) & !is.na(mapa_est)& !is.na(tempo)& !is.na(desemp)) dim(cd) #MODELO DE REGRESSÃO glm3<-glm(desemp~bsc+champion+ap_dir+niv_tre+cons_ext+reu_f_est+com_est+ ajuste+ind_apr+ind_desemp+mapa_est+tempo+desemp,family=binomial,data=cd) modelstep <- step(glm3, direction = "both") logistic.display(modelstep) summary(modelstep) Estou tentando aplicar uma regressão logística, quero estimar um modelo aonde a variável resposta é o desempenho de determinada empresa, aonde tenho os níveis melhora=1, sem melhora=0 com outras variáveis. No link abaixo tem uma amostra do banco de dados. http://www.datafilehost.com/download-996309b6.html Alguém poderia dar uma avaliada no que estou fazendo? [ ]'s. Edson Lira Estatístico Manaus-Amazonas

2 1

Fwd: Re: Teste no R
by Leonard de Assis 22 Out '12

22 Out '12

Esqueci de postar na lista > median.test<-function(y1,y2, na.rm = TRUE){ + z <- c(y1,y2) + g <- rep(1:2, c(length(y1), length(y2))) + m <- median(z, na.rm) + fisher.test(z<m,g)$p.value + } > turma1=c(72,79,93,91,70,95,82,80,74,86, NA) > turma2=c(78,66,65,84,69,73,71,75,68,90,76) > median.test(turma1,turma2) [1] 0.1983806 []s Leonard de Assis assis <dot> leonard <at> gmail <dot> com Em 22/10/2012 17:40, Zé Augusto escreveu: > Luiz Henrique > > Seu script retorna o erro > "Erro em median.default(dados, na.rm = TRUE) : need numeric data" > no R2.15.0 64bits. > > ##################################################### > > Leonardo > > Desculpe a minha ignorância, mas seu script, para > y1 = turma1 > y2 = turma2 > não retorna nada. > > Grato por qualquer informação. > José Augusto > > > Em 22 de outubro de 2012 17:05, Leonard de Assis > <assis.leonard(a)gmail.com> escreveu: > median.test<-function(y1,y2){ > z<-c(y1,y2) > g <- rep(1:2, c(length(y1),length(y2))) > m<-median(z) > fisher.test(z<m,g)$p.value > } >> Fonte: Lista R >> >> []s >> Leonard de Assis >> assis <dot> leonard <at> gmail <dot> com >> >> Em 22/10/2012 16:36, Luiz Henrique da Conceicao Leal escreveu: >> >> # Dados >> turma1=c(72,79,93,91,70,95,82,80,74,86,NA) >> turma2=c(78,66,65,84,69,73,71,75,68,90,76) >> Notas=data.frame(turma1,turma2) >> >> >> # Teste da mediana >> teste.da.mediana=function(dados,alfa=0.05){ >> M=median(dados,na.rm=TRUE) >> Tabela=table(dados[,1]<=M) >> for(j in 2:ncol(dados)) Tabela=rbind(Tabela,table(dados[,j]<=M)) >> colnames(Tabela)=c("Acima da mediana","Abaixo da mediana") >> rownames(Tabela)=colnames(dados) >> if(length(table(Tabela<=5))==1) {Teste=chisq.test(Tabela,correct = FALSE)} >> else Teste=fisher.test(Tabela) >> Conclusão=if(Teste$p.value<alfa) {"Rejeita-se a hipótese de igualdade das >> medianas"} else "Não se pode rejeitar a hipótese de igualdade das medianas" >> Medianas=apply(Notas,2,median,na.rm=TRUE) >> list(Tabela=Tabela,Teste=Teste,Conclusão=Conclusão,Medianas=Medianas) >> } >> >> teste.da.mediana(Notas,alfa=0.05) >> >> >> Em 22 de outubro de 2012 15:54, Márcia Varela da Silva >> <marcinhavarela(a)hotmail.com> escreveu: >>> Como fazer o "Teste da Mediana" no R? >>> >>> >>> >>> _______________________________________________ >>> 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. >> >> >> >> _______________________________________________ >> 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. >> >> >> >> _______________________________________________ >> 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

Teste no R
by Márcia Varela da Silva 22 Out '12

22 Out '12

Como fazer o "Teste da Mediana" no R?

5 4

[off-topic] Concurso público 03 vagas para Estatística UTFPR
by Walmes Zeviani 22 Out '12

22 Out '12

Saudações, Encaminho link do edital da UTFPR sobre divulgação de abertura de 03 vagas para Professor de Estatística dedicação exclusiva. O concurso tem vagas para outras áreas do conhecimento também. http://www.utfpr.edu.br/concursos/campi/ct/arquivos-gerais/Edital_068_2012_… Aos interessados, boa sorte. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 VoIP: (3361 3600) 1053 1173 e-mail: walmes(a)ufpr.br skype: walmeszeviani twitter: @walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

1 0