- R-br - DINF UFPR

trabalho em R test wilcox ou regressão
by alanarocha＠sapo.pt 18 Jun '13

18 Jun '13

Bom dia tenho a seguinte questão a resolver, alguém me sabe dizer se devo utilizar o teste de wilcoxon emparelhado ou um modelo de regressão linear 1. Avaliar as médias para cada caso, considerando 2 períodos a. Período antes é igual aos Doset-9 ao Doset-3, que passará a ser designado como H b. Período depois é igual aos Doset+3 ao Doset+9, que passará a ser designado como C DoseT-9 DoseT-6 DoseT-3 DoseT+3 DoseT+6 DoseT+9 930,0 900,0 900,0 390,0 1.170,0 2.700,0 1.860,0 930,0 900,0 420,0 840,0 630,0 1.860,0 2.790,0 2.700,0 1.620,0 1.440,0 1.530,0 930,0 930,0 900,0 420,0 840,0 720,0 1.080,0 480,0 870,0 2.970,0 2.700,0 2.520,0 2.790,0 1.590,0 930,0 390,0 720,0 2.340,0 930,0 930,0 900,0 390,0 390,0 390,0 1.860,0 750,0 780,0 720,0 780,0 720,0 3.720,0 3.720,0 3.600,0 840,0 840,0 840,0 780,0 420,0 60,0 390,0 390,0 420,0 2.700,0 2.790,0 2.790,0 720,0 780,0 360,0 930,0 930,0 900,0 420,0 420,0 360,0 2.700,0 3.720,0 3.720,0 1.170,0 1.170,0 540,0 1.860,0 930,0 900,0 420,0 420,0 720,0 ..... no total são 93 dados para cada coluna. a minha ideia é fazer o calculo das médias entre doset-9 e doset+9 entre -6 e +6 e entre -3 e +3. Ana

3 2

erro bar
by geovane barbosa 17 Jun '13

17 Jun '13

OLá Daniel tudo bem eu tenho seguinte conjunto de dados, porém tenho mais de 10.000 dados, como poderia construir o erro bar ao longo dessas idades, pegando a média de peso por idade? idade peso 35 59 35 64 35 87 36 52 36 90 36 45 36 87 37 65 37 49 40 26 40 58 41 63 #################################### #################################### Prof. Geovane Carlos Barbosa UCL - Faculdade do Centro Leste ################################### ################################### ________________________________ De: Daniel Marcelino <dmarcelino(a)live.com> Para: R-br Lista <r-br(a)listas.c3sl.ufpr.br>; geovane barbosa <geovanecb(a)yahoo.com.br> Enviadas: Segunda-feira, 17 de Junho de 2013 15:20 Assunto: Re: [R-br] erro bar Mande um exemplo minimamente reproduzível dos dados. Daniel 2013/6/17 geovane barbosa <geovanecb(a)yahoo.com.br> > > > >Olá pessoal tudo bem > > >como faço para construir um erro bar por cada idade tendo uma variável y sendo por exemplo peso. > > >abraços > > > > > >________________________________ > De: Marcos Silva <marcosfs2006(a)gmail.com> >Para: r-br(a)listas.c3sl.ufpr.br >Enviadas: Segunda-feira, 17 de Junho de 2013 10:39 >Assunto: Re: [R-br] Dúvida - gdata (notação científica) > > > >Veja se algo como options(scipen=30) funciona. > >Abs. > > >Em 17 de junho de 2013 09:15, Sérgio Henrique almeida da silva ju <sergio.edfisica(a)gmail.com> escreveu: > >Amigos >> >> >>Estou lendo um banco do excel usando o gdata >> >> >>require(gdata) >>dados <- read.xls(...) >> >> >>Tenho uma variável com dados numéricos de tamanho maior que 30. Ex: >>00000000000000000000000000000083990244081296600001 >> >> >>Quando leio o banco ela fica em notação científica (8.40e+19), porém preciso trabalhar com ela no formato original, qual comando eu posso usar? >> >> >> >> >>-- >> >>Sérgio Henrique Almeida da Silva Junior >>Doutorando em Epidemiologia em Saúde Pública >>Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ >>http://lattes.cnpq.br/1611345552843383 >>Tel: (21) 68463637 >> >>http://www.linkedin.com/profile/view?id=250437145&trk=tab_pro >> >>_______________________________________________ >>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. >> > > >-- >Marcos F. Silva >http://sites.google.com/site/marcosfs2006 >_______________________________________________ >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. > -- \begin{signature} Daniel Marcelino ☁ dm.silva(a)umontreal.ca ☎ (514) 343 6111 #3799 Skype: d.marcelino ✎ 3200 Jean Brillant, Office C5071 Montreal, QC; H3T 1N8 Canada \end{signature} "Small steps toward a much better world"

2 1

Problemas com Datas
by Rodrigo Bombeiro 17 Jun '13

17 Jun '13

Olá, estou trabalhando com datas e gostaria de entender por que a segunda data resulta em NA no código abaixo: tc <- "01/01/2012 13:42:04" tu <- "23/12/2012 21:07:30" test <- strptime(c(tc,tu), "%m/%d/%Y %H:%M:%S") []´s

3 2

Dúvida - gdata (notação científica)
by Sérgio Henrique almeida da silva ju 17 Jun '13

17 Jun '13

Amigos Estou lendo um banco do excel usando o gdata require(gdata) dados <- read.xls(...) Tenho uma variável com dados numéricos de tamanho maior que 30. Ex: 00000000000000000000000000000083990244081296600001 Quando leio o banco ela fica em notação científica (8.40e+19), porém preciso trabalhar com ela no formato original, qual comando eu posso usar? -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 68463637 http://www.linkedin.com/profile/view?id=250437145&trk=tab_pro

5 5

Re: [R-br] Dúvida - gdata (notação científi ca)
by Sérgio Henrique almeida da silva ju 17 Jun '13

17 Jun '13

O options(scipen=30) e o formatC(a, digits = 30) funcionam, porém ele ignora os 0 que estão na frente, mas eu preciso desses 0. O options(digits = 30) dá erro, pois só permite 22. Em 17 de junho de 2013 11:13, geovane barbosa <geovanecb(a)yahoo.com.br>escreveu: > > > Olá pessoal tudo bem > > como faço para construir um erro bar por cada idade tendo uma variável y > sendo por exemplo peso. > > abraços > > > #################################### > #################################### > Prof. Geovane Carlos Barbosa > UCL - Faculdade do Centro Leste > ################################### > ################################### > ------------------------------ > *De:* "andrebvs(a)bol.com.br" <andrebvs(a)bol.com.br> > *Para:* r-br(a)listas.c3sl.ufpr.br > *Enviadas:* Segunda-feira, 17 de Junho de 2013 10:45 > *Assunto:* Re: [R-br] Dúvida - gdata (notação científica) > > Veja também, se algumas dessas opções funcionam: > > options(digits = 30) > a = 8.40e+19 ; a > > formatC(a, digits = 30) > > sprintf("%0.30f", 0.8 + 0.4) > > *Att.* > *André BVS* > > > > ------------------------------ > Em 17/06/2013 10:39, *Marcos Silva < marcosfs2006(a)gmail.com >* escreveu: > Veja se algo como options(scipen=30) funciona. > > Abs. > > Em 17 de junho de 2013 09:15, Sérgio Henrique almeida da silva ju < > sergio.edfisica(a)gmail.com> escreveu: > > Amigos > > Estou lendo um banco do excel usando o gdata > > require(gdata) > dados > > Tenho uma variável com dados numéricos de tamanho maior que 30. Ex: > 00000000000000000000000000000083990244081296600001 > > Quando leio o banco ela fica em notação científica (8.40e+19), porém > preciso trabalhar com ela no formato original, qual comando eu posso usar? > > > > -- > Sérgio Henrique Almeida da Silva Junior > Doutorando em Epidemiologia em Saúde Pública > Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ > http://lattes.cnpq.br/1611345552843383 > Tel: (21) 68463637 > http://www.linkedin.com/profile/view?id=250437145&trk=tab_pro > > _______________________________________________ > 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. > > > > > -- > Marcos F. Silva > http://sites.google.com/site/marcosfs2006 > > _______________________________________________ > 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. > -- Sérgio Henrique Almeida da Silva Junior Doutorando em Epidemiologia em Saúde Pública Escola Nacional de Saúde Pública Sérgio Arouca - ENSP/FIOCRUZ http://lattes.cnpq.br/1611345552843383 Tel: (21) 68463637 http://www.linkedin.com/profile/view?id=250437145&trk=tab_pro

4 5

renomeando níveis de fatores
by Daniela Recchia 17 Jun '13

17 Jun '13

Boa tarde, tenho uma solucao nao muito inteligente para o que venho perguntar e gostaria de obter algo mais otimizado. Tenho um data frame com 40 variaveis (fatores) e cada um com diversos níveis. Algo por exemplo: Variável satisfacao Variável aprendizado Variável outrosplanos muito satisfeito aprendo todos os dias posso mudar planos nem satisfeito nem insatisfeito aprendo 1 vez na semana preciso consultar o chefe mas posso mudar muito insatisfeito aprendo 1 vez ao mes nao posso mudar nunca aprendo raramente nao aprendo algo novo nunca Tenho por exemplo 10 variáveis do estilo da "satisfacao", 20 da "aprendizado" e 10 da "outrosplanos", agora quero ter nos näiveis números e nao frases, "1","2","3". A princ#ipio pensei em fazer manualmente para cada uma e depois unir todas as variáveis em um novo conjunto de dados: newsatisfacao<-factor(data$satisfacao, levels=c("muito satisfeito","nemsatisfeito nem insatisfeito","muito insatisfeito"),labels=c("1","2","3")) ..... newdata<-cbind(newsatisfacao, newaprendizado...) Assim eu teria o que quero, mas é muito trabalhoso e nao prático! Gostaria de poder renomear os níveis de forma mais dinamica e rápida, de repente sem ter que criar uma nova variável, só renomeá-los. Poder dizer que para todas as variáveis que tem as opcoes "muito satisfeito","nemsatisfeito nem insatisfeito","muito insatisfeito", devem ser renomeadas para "1","2","3" e assim por diante. Abracos, -- Daniela Rodrigues Recchia M.Sc. in Statistics Tel. (Mobil): +49 01578 7588382 Kißlegg - Germany “The standard of success in life is absolutely the amount of joy you feel” Esther & Jerry Hick

2 1

Dado não estruturado
by Vinicius Brito Rocha 17 Jun '13

17 Jun '13

Pessoal, estou querendo começar a analisar dados não estruturados como logs, paginas webs etc. Já dei uma boa olhada em materiais de texmining, livros e aplicações com o pacote tm do R. Alguém mais experiente no assunto recomenda algum pacote para desse tipo de dado? Abs -- *Vinicius Brito Rocha.* *Estatístico e Atuário ** M.Sc. Engenharia de Produção/PO*

3 4

Teste de hipoteses em análise de grupo de experimentos
by Gilenio Borges Fernandes 15 Jun '13

15 Jun '13

Gostaria de contar com ajuda para os seguinte problema - analise conjunta de 4 experimentos em blocos casualizados completos com 3 repetições resp=c(109.30,124.00,100.00,87.45,115.12,93.10,95.60,101.00,98.50, 68.00,47.00,78.00,118.00,141.60,117.70,112.30,91.30,125.65,115.00, 109.00,97.50,104.00,78.00,87.00,117.70,127.20,116.50,108.00,156.70, 93.10,83.00,111.50,114.00,128.00,65.00,84.00,122.50,114.60,129.00, 108.00,99.75,135.40,111.00,140.00,83.50,89.00,88.50,60.00,126.60, 112.00,123.60,103.20,113.00,158.50,130.00,119.00,87.50,74.00,82.00,67.00) trat=gl(5,12,label=paste("N",1:5,sep="")) exp<-rep(gl(4,3,label=paste("E",1:4,sep="")),5) bloco<-rep(gl(3,1,label=paste("bloco",1:3,sep="")),20) dados=data.frame(trat=trat,exp=exp,bloco=bloco,resp=resp) mod = aov(resp ~ trat * exp + exp/bloco, data = dados) summary(mod) A declaração do modelo (mod) com a função aov testa tudo com o resíduo (32 gl). Quero testar tratamento (trat) e experimento (exp) com a a interação trat*exp e esta com o resíduo (32 gl). Como declarar os efeitos no modelo? -- Gilenio Borges Fernandes Professor Associado III Universidade Federal da Bahia Instituto de Matemática Departamento de Estatística Av. Adhemar de Barros, s/n – Ondina. 40.170-110 - Salvador - BA, Brasil Tel.: (071)3283-6280/6336 Fax: (071)3283-6276 URL: http://lattes.cnpq.br/6764860618464860

2 1

problema com o pacote coxphw
by Paula Marques Sientchkovski 15 Jun '13

15 Jun '13

Estou encontrando dificuldades em rodar o pacote coxphw ele está apresentando a seguinte mensagem de erro:Houve 16 avisos (use warnings() para vê-los) e os avisos são: 1: In Surv(obj$resp[, 1], obj$resp[, 2], obj$resp[, 3]) : Stop time must be > start time, NA created 2: In min(obskm$surv[obskm$time < obj$resp[i, 2]]) : nenhum argumento não faltante para min; retornando Inf 3: In min(obskm$surv[obskm$time < obj$resp[i, 2]]) : nenhum argumento não faltante para min; retornando Inf 4: In min(obskm$surv[obskm$time < obj$resp[i, 2]]) : nenhum argumento não faltante para min; retornando Inf 5: In min(obskm$surv[obskm$time < obj$resp[i, 2]]) : nenhum argumento não faltante para min; retornando Inf 6: In min(obskm$surv[obskm$time < obj$resp[i, 2]]) : nenhum argumento não faltante para min; retornando Inf 7: In min(obskm$surv[obskm$time < obj$resp[i, 2]]) : nenhum argumento não faltante para min; retornando Inf 8: In min(obskm$surv[obskm$time < obj$resp[i, 2]]) : nenhum argumento não faltante para min; retornando Inf 9: In Surv(obj$resp[, 1], obj$resp[, 2], obj$resp[, 3]) : Stop time must be > start time, NA created 10: In min(obskm$surv[obskm$time < obj$resp[i, 2]]) : nenhum argumento não faltante para min; retornando Inf 11: In min(obskm$surv[obskm$time < obj$resp[i, 2]]) : nenhum argumento não faltante para min; retornando Inf 12: In min(obskm$surv[obskm$time < obj$resp[i, 2]]) : nenhum argumento não faltante para min; retornando Inf 13: In min(obskm$surv[obskm$time < obj$resp[i, 2]]) : nenhum argumento não faltante para min; retornando Inf 14: In min(obskm$surv[obskm$time < obj$resp[i, 2]]) : nenhum argumento não faltante para min; retornando Inf 15: In min(obskm$surv[obskm$time < obj$resp[i, 2]]) : nenhum argumento não faltante para min; retornando Inf 16: In min(obskm$surv[obskm$time < obj$resp[i, 2]]) : nenhum argumento não faltante para min; retornando Inf Já tentei ajustar a sintaxe ,mas não obtive êxito. Obrigada. Paula;

1 0

regressão logística
by Alan R Panosso 15 Jun '13

15 Jun '13

Caros amigos, Venho por meio deste pedir boas referências para estudos sobre regressão logística. Artigos, aplicação e livros. Grato a todos Alan De: R-br [mailto:r-br-bounces@listas.c3sl.ufpr.br] Em nome de walmes . Enviada em: terça-feira, 11 de junho de 2013 19:50 Para: r-br(a)listas.c3sl.ufpr.br Assunto: Re: [R-br] [OFF-TOPIC] código R para plotar um coração com animação Ok pessoal, valeu aí por disponibilizar. Ainda bem que a resposta veio a tempo de eu terminar a aula (17:30-19:30) então deu para mostrar a função. Os alunos gostaram. Vai fazer sucesso. Segue o script da referida aula. http://www.leg.ufpr.br/~walmes/ensino/ce083-2013-01/ce083-2013-01-aula15i.R Grato. 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 <http://www.leg.ufpr.br/%7Ewalmes> linux user number: 531218 ==========================================================================

2 2