- 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

RCurl
by Raphael Saldanha 14 Jun '13

14 Jun '13

Olá pessoal, Estou tentando usar o pacote RCurl para logar em um site (plataforma moodle) e pegar alguns relatórios de acesso, mas não estou tendo sucesso. Alguém tem experiência com este pacote? --- Atenciosamente,Raphael Saldanha rfsaldanha(a)outlook.com

2 4

arquivo netcdf
by Heloíse Pavanato 13 Jun '13

13 Jun '13

Olá, Tenho tentado abrir um arquivo netcdf mas sem êxito. Tentei com o pacote ncdf para abrir um arquivo através de um link OPeNDAP dataset <- open.ncdf(' http://dap.marinexplore.org/dap/u/911215431c07/02841d66-d38e-11e2-8d83-2200… ') e através do RNetCDF dataset <- open.nc(' http://dap.marinexplore.org/dap/u/911215431c07/02841d66-d38e-11e2-8d83-2200… ') e sempre dá o Error in open.ncdf trying to open file..... Estou trabalhando em ambiente Windows. Isso seria um problema? Alguém tem alguma dica para me dar? Grata, Heloise

5 8

Como ordenar os elementos das linhas de uma matriz de dados?
by Gilenio Borges Fernandes 12 Jun '13

12 Jun '13

Prezados(as) Como ordenar os elementos das linhas de uma matriz de dados? Exemplo, 10 linhas e 7 colunas: C1 C2 C3 C4 C5 C6 C7 L1 28 33 36 69 35 45 62 L2 21 28 33 36 16 22 59 L3 23 37 40 61 22 26 30 L4 20 61 74 78 15 22 62 L5 13 51 61 63 7 41 45 L6 13 20 21 42 4 41 62 L7 13 21 28 43 6 8 12 L8 32 42 43 44 45 48 50 L9 5 20 69 74 1 50 52 L10 20 36 51 78 12 49 79 Grato Gilenio -- 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

4 3

Atualizar R 2.14 para 2.15
by Alexandre Santos 12 Jun '13

12 Jun '13

Boa tarde pessoal, Consegui instalar rgdal no R 3.0.0 por sugestões de colegas da lista, mas começaram a aparecer outros pacotes que estão dando problema com o R 3.0.0, como não pretendo entrar em contato com cada um dos desenvolvedores de pacotes que uso, voltei para o R 2.14.1, mas para funcionar alguns pacotes preciso do R 2.15.1. Então pergunto se tem alguma forma de atualizar do R 2.14 para R 2.15 no Ubuntu 12.04, sem passar pelo: sudo apt-get install r-base Que vai me instalar novamente ao R 3.0.0. Tentei baixar o r-base-core_2.15.1-1precise2_amd64.deb e instalar via central de programas do Ubuntu, mas deu conflito com r-base-core, teria alguma solução para o que eu estou tentando fazer? Obrigado, Alexandre

3 4

Web Data mining com R
by Fabrício Barth 12 Jun '13

12 Jun '13

Pessoal, Vai acontecer neste sábado à tarde um tutorial sobre Web Data Mining com R em São Paulo (capital). Este tutorial faz parte de um evento da Mozilla http://www.bandtec.com.br/index.php/vida-no-campus/mozilla-tech-day/ O tutorial é gratuíto. Abraços, -- Fabrício J. Barth http://fbarth.net.br http://blog.fbarth.net.br http://twitter.com/fbarth http://www.flickr.com/fbarth

1 0

média entre pontos de uma mesma matrix
by Tito Conte 12 Jun '13

12 Jun '13

Bom dia, Preciso calcular a média de pontos entre matrizes que possuem NA i gnorando-os exemplo: média entre estas matrizes [,1] [,2] [,3] [1,] 3 6 9 [2,] 4 7 10 [3,] 5 8 11 [,1] [,2] [,3] [1,] 1.5 3.0 4.5 [2,] 2.0 3.5 5.0 [3,] 2.5 4.0 5.5 [,1] [,2] [,3] [1,] 2 2 NA [2,] 2 2 2 [3,] 2 2 NA e obter um resultado que fosse uma matriz do mesmo tamnho com a média de cada ponto ignoradno NAs alguém pode me ajudar com essa? Tito Conte

3 11

filtrando dados
by Hélio Gallo Rocha 12 Jun '13

12 Jun '13

Bom dia a todos, Estou com uma dúvida simples, mas não achei na lista, se puderem ajudar, agradeço desde já: ponto=rep(1:10,3) quadra=c(rep(1,10),rep(2,10),rep(3,10)) y=rnorm(1:30) dados=cbind(quadra,ponto,y);dados Dúvida: Teria como selecionar os dados de "y" retirando os dados da quadra 1? -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

3 7

Selecionar valores de um objeto
by Thiago de Andrade 12 Jun '13

12 Jun '13

Boa noite, amigos. Estou trabalhando com o pacote ltm com modelagens de TRI. Tenho a seguinte saída > EscoresAP<-factor.scores(fit_AfetosPositivosR, resp.patterns=DadosAfetosPositivosReduzido) > EscoresAP Call: grm(data = DadosAfetosPositivosReduzido, Hessian = T) Scoring Method: Empirical Bayes Factor-Scores for specified response patterns: Item4 Item7 Item17 Item18 Item19 Obs Exp z1 se.z1 1 1 1 1 1 1 101 89.462 -1.663 0.408 2 6 6 4 2 2 1 0.002 0.181 0.338 3 1 2 2 2 2 4 2.297 -0.809 0.210 4 1 1 1 1 1 101 89.462 -1.663 0.408 5 1 1 1 1 1 101 89.462 -1.663 0.408 6 2 4 2 3 1 1 0.350 -0.552 0.231 Preciso manusear os valores da coluna "z1" (escores dos sujeitos 1 a 6 no teste), realizando uma série de análises com esses valores. Gostaria de saber se há alguma função por meio da qual eu possa extrair ou selecionar os dados dessa coluna, já que eu tenho mais de 1300 casos e não é muito prático copiar um a um. Agradeço desde já. -- Thiago de Andrade Mestrando em Psicologia Social, PPGP-UFPB Núcleo de Pesquisas em Desenvolvimento Sócio-Moral (NPDSM)

4 3

Separar uma matriz por trechos completos de uma sequência
by Cristiano Santos 12 Jun '13

12 Jun '13

Olá pessoal, Alguém conhece uma forma fácil de implementar isso: Eu tenho essa matriz: cbind(c(1,2,3,7,8,9,10,15,16,17), c(1,1,1,1,1,1,1,1,1,1)) E quero separar a matriz pela sequência da primeira coluna, ou seja, obter uma matriz para cada trecho completo da sequência: list(cbind(c(1,2,3), c(1,1,1)), cbind(c(7,8,9,10), c(1,1,1,1)), cbind(c(15,16,17), c(1,1,1))) Ats, Cristiano Santos.

1 0

RES: código R para plotar um coração com animação
by Alan R Panosso 11 Jun '13

11 Jun '13

Caro Walmes, Genial... Amanhã começo funções e o R com o pessoal da graduação, posso utilizar esse seu exemplo para motiva-los? Atenciosamente 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 ==========================================================================

1 0

[OFF-TOPIC] código R para plotar um coração com animação
by walmes . 11 Jun '13

11 Jun '13

Saudações pessoal, Alguém tem cópia daquele código R, uma vez compartilhado na lista, acho que na época do R-STAT ainda. O código que falo é aquele que tem um laço for(), as cores mudam com a execução. Dou aula sobre funções hoje e quinta e queria motivar os alunos à presentear suas respectivas com um CMR amoroso. Ou provocar logo a separação dependendo do senso de humor e diferença de nerdismo do casal. Versões simples eu encontrei em: http://stackoverflow.com/questions/8082429/plot-a-heart-in-r http://blog.revolutionanalytics.com/2013/02/make-a-valentines-heart-with-r.… http://rud.is/b/2013/02/14/happy-valentines-day-mrshrbrmstr/ 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 linux user number: 531218 ==========================================================================

3 3

Problema para instalar pacote rgdal no R 3.0.0 em Ubuntu 12.04
by ASANTOS 11 Jun '13

11 Jun '13

Boa noite pessoal, Não consigo instalar o pacote rgdal no R 3.0.0 no Ubuntu, quando utilizo: install.packages("rgdal") ocorre o segundo problema: ** preparing package for lazy loading Error : package 'lattice' was built before R 3.0.0: please re-install it ERROR: lazy loading failed for package 'sp' * removing '/home/asantos/R/x86_64-pc-linux-gnu-library/3.0/sp' Warning in install.packages : installation of package 'sp' had non-zero exit status ERROR: dependency 'sp' is not available for package 'rgdal' * removing '/home/asantos/R/x86_64-pc-linux-gnu-library/3.0/rgdal' Warning in install.packages : installation of package 'rgdal' had non-zero exit status The downloaded source packages are in '/tmp/Rtmp5gpcLF/downloaded_packages' Alguém poderia me ajudar e ver por que esta ocorrendo este erro Obrigado -- ====================================================================== Alexandre dos Santos Proteção Florestal Coordenador do curso Técnico em Florestas Vice Coordenador do curso de Engenharia Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br ======================================================================

2 1

lm e NA
by Daniela Recchia 11 Jun '13

11 Jun '13

Bom dia, estou com dificuldades em rodar um modelo no R que acredito ser devido aos NA. Apesar do software aparentemente remover os NA's na hora de rodar a funcao lm, quando ajusto meu modelo model<-lm(cpstotalnorm~indconfident+indpride+indenergy+indgrowth+indseethingsnew+indinputstrategy +indtrainingavailable+grppointsofview+grpconflictresol+grprightpeople+grpdifferentviews+grprethink +grplessonss) recebo um erro "error contrasts can be applied only to factors with 2 or more levels", pesquisei sobre o erro e vi que poderiam ser dois casos: ou as variáveis poderiam estar definidas como caracteres e nao como fatores - o que nao é o caso, ou há NA's no conjunto de dados - é o caso. Tentei retirar entao os NA's dos dados para rodr o modelo, porém usando a funcao: datanew<-na.omit(data) o R retira todas as observacoes, e nao só o NA, por exemplo: > datanew$indseethingsnew factor(0) Levels: 1 2 3 4 5 6 7 enquanto se retiro somente os NA's dessa variável, no conjunto de dados original, ele mantém todas as outras observacoes: > c<-na.omit(data$indseethingsnew) > c [1] 2 2 2 3 1 3 1 2 1 1 2 3 3 3 3 2 2 1 3 4 1 1 4 3 2 3 1 2 2 1 5 2 1 3 1 1 1 6 3 1 2 1 1 1 1 2 1 2 1 3 1 [52] 3 2 5 1 5 1 1 3 2 1 2 2 1 2 2 1 1 2 2 1 2 1 2 2 1 1 2 2 3 2 2 1 2 2 2 2 3 1 3 2 2 1 3 1 2 2 1 2 3 2 1 [103] 2 1 2 4 2 1 1 2 2 1 3 3 2 2 1 1 6 2 1 2 4 2 3 3 2 2 2 2 1 3 6 1 1 3 2 3 3 1 2 3 2 2 7 2 3 2 1 3 4 1 2 [154] 3 2 1 3 attr(,"na.action") [1] 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 [26] 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 [51] 85 86 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 [76] 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 [101] 211 212 213 214 215 216 217 218 219 attr(,"class") [1] "omit" Levels: 1 2 3 4 5 6 7 Entao meu novo conjunto de dados datanew contém apenas zeros! Alguém teria alguma sugestao? Obrigada, Daniela -- 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

5 7

Re: [R-br] Remover NA: PCA
by Helena Turon 09 Jun '13

09 Jun '13

No meu caso era um vetor, que eu queria plotar o boxplot sem o NA. Assim poderia aparecer os pontos discrepantes. Mas, você tem razão se é uma matriz afetará o resultado. Aí deu certo. Helena

2 1

Remover NA: PCA
by Vinícius Lionel Mateus 09 Jun '13

09 Jun '13

Prezados, Estou executando o seguinte CRM: read.table(file = "clipboard", header = TRUE, sep = "\t", dec = ",")-> TSF > TSF > names(TSF) [1] "Date" "Label" "PMF" "BC" "Acet" "Form" "Cl." "NO3." "PO43." "SO42." "Na." "K." [13] "Mg2." "Ca2." "NH4." > TSFa <- TSF[,2:15] > model <- prcomp(TSFa, scale = TRUE) Erro em colMeans(x, na.rm = TRUE) : 'x' deve ser numérico > model <- prcomp(TSFa, na.rm = TRUE, scale = TRUE) Erro em colMeans(x, na.rm = TRUE) : 'x' deve ser numérico Tentei diferentes métodos para remover o "NAs", mas não obtive sucesso. Alguém tem alguma dica? -- Vinícius -- -- Atenciosamente, VINÍCIUS LIONEL MATEUS, M.Sc (http://lattes.cnpq.br/6501001637020665) Bacharel em Química - Doutorando em Química Analítica Laboratório de Química Atmosférica - Departamento de Química Pontifícia Universidade Católica - Rio de Janeiro (PUC - Rio) Rua Marquês de São Vicente, 225, Gávea - Rio de Janeiro, RJ - Brasil CEP.: 22451-900 Telefone: (+55) (21) 3527-1327 (+55) (21) 9358-8051 www.puc-rio.br

4 10

Remover NA: PCA
by Helena Turon 08 Jun '13

08 Jun '13

Depois de transformar a variável para numérica eu ordenei, sem NA, deu certo. Tsfa$x = as.numeric(Tsfa$x) Tsfa$x <- sort(Tsfa$x[!is.na(Tsfa$x)]) Boa sorte. Helena

2 1

[Dúvida] Estimação por Máxima Verossimilhança.
by Pedro Rafael 07 Jun '13

07 Jun '13

Senhores tenho uma dúvida. Na verdade não é dúvida, apenas quero sugestões. Tenho uma função densidade de probabilidade "complicada". Trate-se de uma distribuição chamada Kwmarashwamy Weibull Poisson. Tenho alguns bancos de dados e gostaria de verificar o ajustamento dessa distribuição à estes bancos de dados. Estou estimando os parâmetros pelo método de máxima verossimilhança. Essa distribuição tem suporte nos reais positivos (x>0) e todos os seus parâmetros são positivos. Optei em utilizar o método L-BFGS-G para restringir a busca nos reais positivos. Segue abaixo o comando. Nesse exemplo não houve convergência. Percebi que os chutes iniciais influenciam muito as estimativas dos parâmetros no caso em que há convergência. Usando métodos de maximização diferentes em muitos casos há grandes diferenças nas estimativas. Existe alguma forma mais tranquila e direta para encontrar as estimativas pelo método de máxima verossimilhança em R? O código segue abaixo: vero <- function(par,x){ a = par[1] b = par[2] c = par[3] lambda = par[4] beta = par[5] -sum(log((a*b*c*lambda*(beta^c)*(x^(c-1))*((1-exp(-(x*beta)^c))^(a-1)) * ((1-(1-exp(-(beta*x)^c))^a)^(b-1)) * exp(-lambda*(1-(1-(1-exp(-(beta*x)^c))^a)^b) - (beta*x)^c))/(1-exp(-lambda)))) } dados = c(17.23, 28.92, 33.00, 41.52, 42.12, 45.60, 48.80, 51.84, 51.96, 54.12, 55.56, 67.80, 68.64, 68.64,68.88, 84.12, 93.12, 98.64, 105.12, 105.84, 127.92, 128.04, 173.40) optim(par=c(1,1,1,1,1),fn=vero, method="L-BFGS-B",x=dados/1000, lower=c(0.001,0.001,0.001,0.001,0.001), upper=c(Inf,Inf,Inf,Inf,Inf)) O engraçado é que essa distribuição é encaixada com a distribuição Weibull. Quando tento ajustar a Weibull à esses dados há convergência. O que vocês acham que devo fazer? [ ], Pedro Rafael Diniz Marinho.

6 12

Inserir NA em objeto do tipo lista
by Aline Lipsky 07 Jun '13

07 Jun '13

Bom dia, Eu tenho um objeto do tipo lista com mais de 6 mil vetores com tamanhos diferentes.Como faço para acrescentar NA nos vetores de tamanho menor e transformar esse objeto em um data frame? Grata, Aline

2 1

encoding de páginas web readLines XML
by walmes . 07 Jun '13

07 Jun '13

Saudações, Estou fazendo um estudo com vídeos do youtube. A partir dos links eu retiro informações deles, como o título. No entanto, não estou conseguindo acertar o encoding, seja para ler corretamente a página ou para tranformar as strings após a leitura. Com o CMR abaixo eu leio o código de um vídeo do youtube e pego seu título, o problema é que os acentos não são corretamente representados, e eu preciso que sejam. Tentei atribuir um encoding na readLines() mas nenhum dos usados (latin1 e utf-8) funcionaram. Tentei usar a iconv() mas sem sucesso. Uso ubuntu 12.04. Sessão R está no CMR. library(XML) # lê a página do vídeo com esse endereço lines <- readLines("http://www.youtube.com/watch?v=P3m6rJR2yDI") # lines <- readLines("http://www.youtube.com/watch?v=P3m6rJR2yDI", encoding=?) # organiza h <- htmlTreeParse(lines, asText=TRUE, useInternalNodes=TRUE) summary(h) # pega o título do vídeo n <- unlist(getNodeSet(doc=h, path="//title", fun=xmlValue)) n Encoding(n) # retira os caracteres estranhos substituindo por vazio iconv(nap, from="UTF-8", to="ASCII", sub="") # iconv(n, from=?, to=?) *# Como sai # "MatemÃ¡tica - EstatÃstica: Conceito e frequÃªncias - YouTube" # Como eu gostaria # "Matemática - Estatística: Conceito e frequências - YouTube"* ## > sessionInfo() ## R version 3.0.1 (2013-05-16) ## Platform: i686-pc-linux-gnu (32-bit) ## locale: ## [1] LC_CTYPE=pt_BR.UTF-8 LC_NUMERIC=C ## [3] LC_TIME=pt_BR.UTF-8 LC_COLLATE=pt_BR.UTF-8 ## [5] LC_MONETARY=pt_BR.UTF-8 LC_MESSAGES=pt_BR.UTF-8 ## [7] LC_PAPER=C LC_NAME=C ## [9] LC_ADDRESS=C LC_TELEPHONE=C ## [11] LC_MEASUREMENT=pt_BR.UTF-8 LC_IDENTIFICATION=C ## attached base packages: ## [1] stats graphics grDevices utils datasets methods base ## other attached packages: ## [1] XML_3.96-1.1 ## loaded via a namespace (and not attached): ## [1] compiler_3.0.1 tools_3.0.1 ## > Alguém saberia solucionar esse problema? Desde já 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 linux user number: 531218 ==========================================================================

4 4

Adicionar Na em uma lista
by Aline Lipsky 06 Jun '13

06 Jun '13

Boa tarde, Eu tenho um objeto do tipo lista com mais de 6 mil vetores com tamanhos diferentes. Eu consigo acrescentar o NA em um elemento da lista: especies[[1]][[7]] <- NA Se alguém puder me ajudar, como faço para repetir esse passo em todos os elementos da minha lista? Tentei com a função sapply, mas não dá certo. Grata, Aline

5 5

MANTEL
by Cézar Borges 05 Jun '13

05 Jun '13

Olá caros,Realizei um teste de mantel para 19 fragmentos para verificar a dependência espacial da comunidade arbórea, o resultado (r=-0.01; p=0.55) me mostrou que os fragmentos são independentes espacialmente. No entanto, estou com a seguinte dúvida: posso interpretar este resultado como se todos os fragmentos fossem diferentes entre si??? ou se não, como posso estar verificando quais os fragmentos que são iguais e diferentes entre si???Abraços! Att, Cézar Fonseca Borges Mestrando em Ciências Ambientais UFPA/MPEG/EMBRAPA - Amazônia Oriental Geógrafo Email: cezar-fonseca(a)hotmail.com/ caborges(a)museu-goeldi.br Lattes: http://lattes.cnpq.br/9805721240137485

1 0

teste de Jonckheere
by Cristina 04 Jun '13

04 Jun '13

Olá pessoal! Gostaria de saber se alguém já usou, no R, o teste de Jonckheere. tentei usar o teste mas consigo encontrar qual pacote está instalado. Obrigada Teresa Cristina

3 3

Coordenadas UTM para Lat Lon
by Tito Conte 04 Jun '13

04 Jun '13

Alguém poderia me explicar como faço a converção de utm para lat lon, estou tentando pelo rgdal, mas não consegui entender direito como funciona as funções de conversão Tito Conte

4 4

erro na validação cruzada
by Hélio Gallo Rocha 04 Jun '13

04 Jun '13

Caros da Lista, Está dando este erro ao tentar rodar a validação cruzada na direção de 45 graus: xvalid: performing cross-validation at location ... 1, Erro em solve.default(Vcov, trend.d) : sistema é computacionalmente singular: condição recíproca número = 8.6266e-21 Grato a todos -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

2 2

Remover compilação
by Sérgio Henrique almeida da silva ju 03 Jun '13

03 Jun '13

Amigos Estou usando a versão mais nova do R (3.0.1) no ubuntu. Tive que compilar a versão 2.14, para fazer um trabalho, agora quero remover essa compilação, como faço? Sempre quando chamo o R no terminal aparece a versão mais antiga, como faço pra voltar a chamar o R 3.0.1? Abraços -- 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

2 2

Separação de duas distribuições de dados
by ASANTOS 03 Jun '13

03 Jun '13

Boa tarde pessoal, Simulei duas distribuições normais e representei em um histograma que obviamente assumiu a forma bimodal, sendo: ### set.seed(10) dist1<-rnorm(n=400,mean=10,sd=5) dist2<-rnorm(n=400,mean=35,sd=8) dist.c<-rbind(dist1,dist2) # hist(dist.c) # Gostaria de saber se existe algum procedimento que me permita separar dist1 e dist2 com seus respectivos intervalos de confiança, imaginando que eu não conheça os parâmetros media e desvio padrão que os criaram? Obrigado -- ====================================================================== Alexandre dos Santos Proteção Florestal Coordenador do curso Técnico em Florestas Vice Coordenador do curso de Engenharia Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br ======================================================================

2 1

Plotagem 3D
by Tiago Souza Marçal 03 Jun '13

03 Jun '13

Bom dia pessoal, estou querendo plotar um gráfico tridimensional com os escores canônicos que estão abaixo: dados<-structure(list(V1 = c(-18.7175, -20.1915, -21.1557, -17.5965, -18.142, -13.857, -18.0086, -14.0436, -12.571), V2 = c(38.3004, 36.4955, 34.4964, 39.9303, 40.4819, 38.3763, 34.9869, 35.0723, 35.9022), V3 = c(-20.9237, -20.12, -17.8586, -17.0183, -17.9961, -20.8457, -19.5551, -17.5703, -18.7393)), .Names = c("V1", "V2", "V3"), class = "data.frame", row.names = c(NA, -9L)) require(lattice) zlab <- "VC2" xlab <- "VC1" ylab <- "VC3" wireframe(V2~V1*V3, scales=list(arrows=FALSE), zlab=list(zlab, rot=90), xlab=list(xlab, rot=24), ylab=list(ylab, rot=-37),xlim=c(-12,-20),ylim=c(-18,-22),type=c("on","bottom")) No entanto, não estou conseguindo fazer com que os pontos apareçam no gráfico e gostaria de eliminar as 3 arestas frontais. Agradeço desde já pela ajuda. 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 #################################################################

2 2

gráfico
by alanarocha＠sapo.pt 03 Jun '13

03 Jun '13

Boa tarde, porque será que o gráfico não é desenhado? pode ser por as medianas serem iguais? código: setwd("G:/SIN/Users/Ana Rua/dados não Panakeia e Panakeia") dados<-read.table("Dadoshgbnaophgbp2sem2011.csv",header=TRUE,sep=";", stringsAsFactors=TRUE, dec=",") dados #objectivo é: #verificar se há diferenças estatisticamente significativas entre os custos e os valores de hemoglobina para clinicas não panakeia no 2semestre de 2011. # Imprime as hipoteses do teste cat("\n\t Teste Wilcoxon\n") cat("Teste bicaudal:\n\tH0: há diferenças significativas entre os valores de hemoglobina para clínicas não Panakeia e Panakeia, mediana=0", "\n\tH1: Não há diferenças significativas entre os valores de hemoglobina para clínicas não Panakeia e Panakeia, mediana diferente de 0 \n") wilcox.test(dados$hgbnaopanakeia, dados$hgbpanakeia, paired=TRUE) boxplot(dados, axes=FALSE) # O parametro axes=FALSE faz com que os eixos não sejam desenhados title(col.main='blue',main= 'Valores de hemoglobina das clínicas não Panakeia e Panakeia no 2semestre de 2011') # Acrescenta o título do gráfico title(col.lab='blue',xlab= 'meses do 2semestre do ano 2011') # Acrescenta o titulo do eixo X title(col.lab="blue", ylab= 'Valores') # Acrescenta o titulo do eixo Y box() # Desenha a caixa em torno do gr?fico axis(2) # Coloca o eixo Y (o 2 ali representa o eixo Y) no gráfico axis(1, 1:6, c('Jul', 'Ago', 'Set', 'Out', 'Nov', 'Dez')) # Coloca no eixo X (representado pelo 1), nos pontos 1 a 6, as etiquetas Jul a Dez. dados: X2sem2011 hgbnaopanakeia hgbpanakeia 1 JUL 75.9 75.4 2 AGO 77.5 77.1 3 SET 79.4 77.0 4 OUT 76.8 77.2 5 NOV 76.1 78.0 6 DEZ 78.3 76.7 Ana

2 1

MObilidade urbana
by Edson Lira 03 Jun '13

03 Jun '13

Caros amigos, estou tentando estabalecer uma amostra para o seguinte problema: Tenho uma origem e quatro destinos. Na origem consigo coletar 1326 placas de carros em 60 minutos, São 4 destinos e consigo identificar 620 placas desta origem, como estabelecer uma amostra nos 4 destinos para identificar o origem? [ ]'s. Edson Lira Estatístico Manaus-Amazonas

2 1

Seleção em objeto tipo list
by ASANTOS 02 Jun '13

02 Jun '13

Boa noite pessoal, Uso muito pouco objetos do tipo list e tenho uma pergunta trivial, quando faço: df <- data.frame(sample=rnorm(1,0,1),params=I(list(list(mean=0,sd=1,dist="Normal")))) df <- rbind(df,data.frame(sample=rgamma(1,5,5),params=I(list(list(shape=5,rate=5,dist="Gamma"))))) df <- rbind(df,data.frame(sample=rbinom(1,7,0.7),params=I(list(list(size=7,prob=0.7,dist="Binomial"))))) df <- rbind(df,data.frame(sample=rnorm(1,2,3),params=I(list(list(mean=2,sd=3,dist="Normal"))))) df <- rbind(df,data.frame(sample=rt(1,3),params=I(list(list(df=3,dist="Student-T"))))) # df[[2]] ## Objeto 2 novo<-df[[2]]=="Binomial" ## Não funciona novo<-df[[2=="Binomial",]] ## Tambem não funciona Quero selecionar somente "Binomial" dentro de [[2]] em novo objeto chamado novo e não estou conseguindo, como faço isto? Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal Coordenador do curso Técnico em Florestas Vice Coordenador do curso de Engenharia Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br ======================================================================

2 1

gráfico para teste de wilcoxon emparelhado
by alanarocha＠sapo.pt 02 Jun '13

02 Jun '13

Bom dia não percebo porque o objecto não é encontrado. codigo: boxplot(dados, axes=FALSE) # O parametro axes=FALSE faz com que os eixos não sejam desenhados title(col.main='blue',main= 'Custos em ? das clínicas não Panakeia e Panakeia') # Acrescenta o título do gráfico title(col.lab='blue',xlab= 'meses do Ano 2011') # Acrescenta o titulo do eixo X title(col.lab="blue", ylab= 'Valores') # Acrescenta o titulo do eixo Y box() # Desenha a caixa em torno do gr?fico axis(2) # Coloca o eixo Y (o 2 ali representa o eixo Y) no gráfico axis(1, 1:6, Jul:Dez) # Coloca no eixo X (representado pelo 1), nos pontos 1 a 6, as etiquetas Jul a Dez. Error in as.graphicsAnnot(labels) : object 'Jul' not found Cumprimentos Ana Quoting r-br-request(a)listas.c3sl.ufpr.br: > 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: DE - PARA (Tito Conte) > 2. Re: Detecçção de pontos externos (Tito Conte) > 3. Re: Data Compare (Epi info) > (Pedro Emmanuel Alvarenga Americano do Brasil) > 4. Previsão ARIMA (Paola Tame) > 5. gráfico para teste wilcoxon (alanarocha(a)sapo.pt) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Wed, 29 May 2013 14:42:40 -0300 > From: Tito Conte <tito.conte(a)gmail.com> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br>, Edmar > Caldas <edmar_caldas(a)yahoo.com.br> > Subject: Re: [R-br] DE - PARA > Message-ID: > <CACqq46xog5kXTFOg3W_ba84OAD-y3Zu=s=my3WyDLXUmaQbmcw(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Não entendi o que você quer edmar > > Tito Conte > > > > 2013/5/28 Edmar Caldas <edmar_caldas(a)yahoo.com.br> > >> Bom Dia! >> >> Alguém poderia me dizer se dá pra fazer um de - para. >> >> Ex: >> >> Tenho uma tabela: >> >> Cod_item Produto Jan Fev mAR ABRIL >> maio junho jlho agosto set out >> nov dez >> 001 AAA 10 15 10 8 >> 9 0 0 0 00 >> 0 0 0 >> 002 BBB 0 0 0 >> 0 0 5 4 >> 5 4 4 4 4 >> >> agora quero fazer que o cod_item 002 que tem valor diferente de ZERO >> paraa para a serie do produto AAA fica assim então. >> >> Cod_item Produto Jan Fev >> mAR ABRIL maio junho >> jlho aosto set out nov >> dez >> 001 AAA 10 15 >> 10 8 9 >> 5 4 5 4 4 >> 4 4 >> >> >> >> >> Edmar >> >> _______________________________________________ >> 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. >> > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: > <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130529/b1e88b2c/atta…> > > ------------------------------ > > Message: 2 > Date: Wed, 29 May 2013 14:45:28 -0300 > From: Tito Conte <tito.conte(a)gmail.com> > To: "r-br(a)listas.c3sl.ufpr.br" <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Detecçção de pontos externos > Message-ID: > <CACqq46z8k11XMrNvzGKG1=Z1ENG4CD3DSgxtoshnu+b3a3SLYA(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Então na verdade o problema com estes dados foram levantados por terceiros > e repassados. > O problema é que estou questionando os resultados, e quera eu mesmo > interpola-los, mas ao que parece esses já são os interpolados. Pra esse > caso solicitei os brutos. Estou aguardando. > > Agradeço a sua disposição e discussão, me ajudaram a criticar um pouco o > que tenho recebido desse pessoal. > > Tito Conte > > > > 2013/5/27 Hélio Gallo Rocha <heliogallorocha(a)gmail.com> > >> Tito, >> >> Andei olhando, parece que o comando over extrai dados como uma seleção, >> >> Talvez tenha me expresso mal, o que tava pensando era detectar os pontos >> mais externos da área, para fazer a borda na krigagem. >> >> Como estão os estudos da krigagem? >> >> Consegui ver seus dados usando o pacote RgoogleMaps, é em são luis do >> maranhão? >> >> mas pra te dizer a verdade, não tem o que krigar ali, já tem informação >> demais ali. >> >> >> >> Hélio >> >> >> Em 25 de maio de 2013 11:54, TitoConte [via R-br] < >> ml-node+s2285057n4659465h83(a)n4.nabble.com> escreveu: >> >>> Funcao over do pacote sp ( desculpe sofri correção automática) >>> >>> >>> On Saturday, May 25, 2013, Tito Conte wrote: >>> >>>> Overdose do pacote sp >>>> >>>> On Saturday, May 25, 2013, Hélio Gallo Rocha wrote: >>>> >>>>> Bom dia a todos da lista, >>>>> >>>>> Pegando o exemplo do Tito onde seus dados possuem 25000 coordenadas, >>>>> teria uma função onde se possa separar apenas os pontos externos >>>>> para fazer >>>>> a borda? >>>>> >>>>> Grato >>>>> >>>>> -- >>>>> Hélio Gallo Rocha >>>>> IFSULDEMINAS - Câmpus Muzambinho >>>>> >>>> >>>> >>>> -- >>>> Tito Conte >>>> >>>> >>>> >>> >>> -- >>> Tito Conte >>> >>> >>> >>> _______________________________________________ >>> R-br mailing list >>> [hidden email] <http://user/SendEmail.jtp?type=node&node=4659465&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-Detecccao-de-pontos-externos-tp46594… >>> 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…> >>> >> >> >> >> -- >> Hélio Gallo Rocha >> IFSULDEMINAS - Câmpus Muzambinho >> >> _______________________________________________ >> 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. >> > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: > <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130529/5e66d46d/atta…> > > ------------------------------ > > Message: 3 > Date: Wed, 29 May 2013 16:20:06 -0300 > From: Pedro Emmanuel Alvarenga Americano do Brasil > <emmanuel.brasil(a)gmail.com> > To: R BR lista <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] Data Compare (Epi info) > Message-ID: > <CAFfGvyJGHLq7aK8F17_s8FfgvTCQ187tEe4Rdmiy2f5DJjyHjg(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Juan, > > Interessante essa sua questão. Eu ja fiz isso no EpiData e no Stata ha > muitos anos mas nunca fiz no R. Procurei rapidamente no R-sitesearch e > apreceram algumas coisas. > > Dei uma lida rápida mas não tenho certeza se alguma coisa faz exatamente o > que voce quer. Se nenhuma desses links te ajudar eu te ajudo. > > http://finzi.psych.upenn.edu/R/library/compare/doc/compare-intro.pdf > http://finzi.psych.upenn.edu/R/library/compare/html/00Index.html > http://finzi.psych.upenn.edu/R/library/compare/html/compareEqual.html > > > Dr. Pedro Emmanuel A. A. do Brasil > http://blog.ipec.fiocruz.br/lapclin-chagas/ > Curriculum Lattes: http://lattes.cnpq.br/6597654894290806 > ResearchGate.net: https://www.researchgate.net/profile/Pedro_Brasil2/ > Instituto Nacional de Infectologia/Instituto de Pesquisa Clínica Evandro > Chagas > Fundação Oswaldo Cruz > Rio de Janeiro - Brasil > Av. Brasil 4365, > CEP 21040-360, > Tel 55 21 3865-9648 > e-mail: pedro.brasil(a)ipec.fiocruz.br > e-mail: 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 > > > Em 24 de maio de 2013 15:05, Sérgio Henrique almeida da silva ju < > sergio.edfisica(a)gmail.com> escreveu: > >> Obrigado Benilton >> >> Desculpe tê-lo feito digitar. >> >> Mas eu gostaria que a função me retornasse além do ID do erro de dgitação, >> mas também os valores discordantes, isso tudo em uma matrix ou um >> data.frama. >> >> Obrigado mais uma vez >> >> Abraços >> >> >> Em 24 de maio de 2013 15:42, Benilton Carvalho <beniltoncarvalho(a)gmail.com >> > escreveu: >> >> Poxa... um exemplo reproduzivel ajudaria tanto quem tirar o proprio >>> tempo para poder te ajudar.... (leia-se: se eu tenho q digitar >>> comandos para ter o seu conjunto de dados de exemplo, entao nao e' um >>> exemplo reproduzivel.... caso ideal para um consultor, que recebe por >>> hora trabalhada....) >>> >>> df1 = data.frame(ID=1:5, idade=c(15, 20, 18, 19, 20), sexo=c('M', 'M', >>> 'F', 'F', 'M')) >>> df2 = data.frame(ID=1:5, idade=c(15, 21, 18, 19, 20), sexo=c('M', 'M', >>> 'F', 'M', 'M')) >>> which(!apply(df1 == df2, 1, all)) >>> >>> b >>> >>> Em 24 de maio de 2013 14:36, Sérgio Henrique almeida da silva ju >>> <sergio.edfisica(a)gmail.com> escreveu: >>> > Prezados >>> > >>> > Alguém tem algum script no R que compara dois bancos de dados, como no >>> data >>> > compare do epiinfo ou validate do epidata? >>> > >>> > De qualquer forma eles comparam os dados de dois bancos de dados, >>> contendo >>> > as mesmas variáveis, mas que foram digitados por pessoas diferentes. >>> > >>> > Ex.: >>> > >>> > Digitador 1 (Banco1) >>> > >>> > ID idade sexo >>> > 1 15 M >>> > 2 20 M >>> > 3 18 F >>> > 4 19 F >>> > 5 20 M >>> > >>> > Digitador 2 (Banco2) >>> > >>> > ID idade sexo >>> > 1 15 M >>> > 2 21 M >>> > 3 18 F >>> > 4 19 M >>> > 5 20 M >>> > Nesse caso houve discordância na digitação no ID 2 e no 4. >>> > >>> > Gostaria de fazer uma função que me aponte essas discordância, alguém >>> tem >>> > alguma ideia? >>> > >>> > Abraços >>> > >>> > >>> > -- >>> > 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. >>> _______________________________________________ >>> 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 >> >> _______________________________________________ >> 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. >> > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: > <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130529/25b57a27/atta…> > > ------------------------------ > > Message: 4 > Date: Wed, 29 May 2013 19:12:43 -0300 > From: Paola Tame <paola.tame(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Previsão ARIMA > Message-ID: > <CAGT6hUtY3k_EUnX4MiiO+VD+i-d-hu4HkwKWNNUCTGPXQDGxpg(a)mail.gmail.com> > Content-Type: text/plain; charset="iso-8859-1" > > Boa noite pessoal, > > estou com dúvida conceitual (e um tanto quanto básica) com relação ao > cálculo do resíduo de um modelo ARIMA(0,0,1). Gerei o modelo da seguinte > forma: fit <- arima(dados, order=c(0,0,1)). Para gerar os resíduos utilizei > o fit$res. Foi exibido o resíduo para todas as observações do banco de > dados. > > Minha dúvida é: como é calculado o valor previsto para o primeiro item, já > que ele depende do erro do modelo no período anterior? > > Att., > -- > Paola Tame > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: > <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20130529/a6ce859e/atta…> > > ------------------------------ > > Message: 5 > Date: Thu, 30 May 2013 14:03:39 +0100 > From: alanarocha(a)sapo.pt > To: R-br(a)listas.c3sl.ufpr.br > Subject: [R-br] gráfico para teste wilcoxon > Message-ID: > <20130530140339.Horde.km58b5IzaPJRp04rLNO3MJA(a)mail.sapo.pt> > Content-Type: text/plain; charset=UTF-8; format=flowed; DelSp=Yes > > olá, > queria perguntar porque é que o grafico não é desenhado? > boxplot(dados], axes=FALSE) # O parametro axes=FALSE faz com que os > eixos n?o sejam desenhados > title(col.main='blue',main= 'custos clínicas não panakeia e panakeia') > # Acrescenta o t?tulo do gr?fico > title(col.lab='blue',xlab= 'Ano') # Acrescenta o titulo do eixo X > title(col.lab="blue", ylab= 'custos em euros') # Acrescenta o titulo > do eixo Y > box() # Desenha a caixa em torno do gr?fico > axis(2) # Coloca o eixo Y (o 2 ali representa o eixo Y) no gr?fico > axis(1, 1:3, 2012) # Coloca no eixo X (representado pelo 1), nos > pontos 1 a 5, as 2012. > > > cumprimentos Ana > > > ------------------------------ > > Subject: Legenda do Digest > > _______________________________________________ > 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 29, assunto 29 > *****************************************

3 2

Data Compare (Epi info)
by Sérgio Henrique almeida da silva ju 01 Jun '13

01 Jun '13

Prezados Alguém tem algum script no R que compara dois bancos de dados, como no data compare do epiinfo ou validate do epidata? De qualquer forma eles comparam os dados de dois bancos de dados, contendo as mesmas variáveis, mas que foram digitados por pessoas diferentes. Ex.: Digitador 1 (Banco1) ID idade sexo 1 15 M 2 20 M 3 18 F 4 19 F 5 20 M Digitador 2 (Banco2) ID idade sexo 1 15 M 2 21 M 3 18 F 4 19 M 5 20 M Nesse caso houve discordância na digitação no ID 2 e no 4. Gostaria de fazer uma função que me aponte essas discordância, alguém tem alguma ideia? Abraços -- 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

Fwd: ks.test e dados efeitos alelopáticos
by Polliana Zocche 30 Mai '13

30 Mai '13

Olá, Eu estou procurando que teste seria apropriado para dados que contém zeros em seus vetores, vetores com valores repetidos e vetores com valores crescentes até 100. Eu usei o ks.test() para verificar as diferentes entre número de germinações/dia entre diferentes tratamentos avaliando alelopatia, mas o teste retorna um warning por causa dos valores repetidos (ties): "Mensagens de aviso perdidas: In ks.test(alelo$control, alelo$ph5.7, alelo$v100, alelo$v200, alelo$v200ph, : cannot compute exact p-values with ties" Alguém tem uma sugestão de que teste seria apropriado para fazer comparações com dados dessa natureza? Muito obrigada! -- Polliana Zocche de Souza Bióloga/Mestre em Ecologia Doutoranda em Ecologia Departamento de Biologia Vegetal, IB, UNICAMP

2 5

gráfico para teste wilcoxon
by alanarocha＠sapo.pt 30 Mai '13

30 Mai '13

olá, queria perguntar porque é que o grafico não é desenhado? boxplot(dados], axes=FALSE) # O parametro axes=FALSE faz com que os eixos n?o sejam desenhados title(col.main='blue',main= 'custos clínicas não panakeia e panakeia') # Acrescenta o t?tulo do gr?fico title(col.lab='blue',xlab= 'Ano') # Acrescenta o titulo do eixo X title(col.lab="blue", ylab= 'custos em euros') # Acrescenta o titulo do eixo Y box() # Desenha a caixa em torno do gr?fico axis(2) # Coloca o eixo Y (o 2 ali representa o eixo Y) no gr?fico axis(1, 1:3, 2012) # Coloca no eixo X (representado pelo 1), nos pontos 1 a 5, as 2012. cumprimentos Ana

2 1

Previsão ARIMA
by Paola Tame 29 Mai '13

29 Mai '13

Boa noite pessoal, estou com dúvida conceitual (e um tanto quanto básica) com relação ao cálculo do resíduo de um modelo ARIMA(0,0,1). Gerei o modelo da seguinte forma: fit <- arima(dados, order=c(0,0,1)). Para gerar os resíduos utilizei o fit$res. Foi exibido o resíduo para todas as observações do banco de dados. Minha dúvida é: como é calculado o valor previsto para o primeiro item, já que ele depende do erro do modelo no período anterior? Att., -- Paola Tame

1 0

Detecçção de pontos externos
by Hélio Gallo Rocha 29 Mai '13

29 Mai '13

Bom dia a todos da lista, Pegando o exemplo do Tito onde seus dados possuem 25000 coordenadas, teria uma função onde se possa separar apenas os pontos externos para fazer a borda? Grato -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

2 4

DE - PARA
by Edmar Caldas 29 Mai '13

29 Mai '13

Bom Dia! Alguém poderia me dizer se dá pra fazer um de - para. Ex: Tenho uma tabela: Cod_item Produto Jan Fev mAR ABRIL maio junho jlho agosto set out nov dez 001 AAA 10 15 10 8 9 0 0 0 00 0 0 0 002 BBB 0 0 0 0 0 5 4 5 4 4 4 4 agora quero fazer que o cod_item 002 que tem valor diferente de ZERO paraa para a serie do produto AAA fica assim então. Cod_item Produto Jan Fev mAR ABRIL maio junho jlho aosto set out nov dez 001 AAA 10 15 10 8 9 5 4 5 4 4 4 4 Edmar

2 1

adicionando pontos com mouse
by Hélio Gallo Rocha 28 Mai '13

28 Mai '13

Caros , Seria possível adicionar pontos em um mapa, usando o mouse com posterior extração das coordenadas? -- Hélio Gallo Rocha IFSULDEMINAS - Câmpus Muzambinho

2 2

Uso de vetores grandes
by Itamar Santos 28 Mai '13

28 Mai '13

Bom dia a todos da lista! Eu uso o R 64 bits em uma máquina com 24Gb de RAM, mas, mesmo assim, por vezes o R não aloca vetores grandes. Na última vez que ocorreu isso, era um vetor de 3.7Gb logo no início da sessão. Consegui resolver aquele problema mudando a estratégia de programação, mas a princípio julgava que isso não aconteceria, já que o R é capaz de reconhecer bem mais memória RAM que isso. Provavelmente isto está ligado a algum setup. Tentei procurar a resposta, sem sucesso. Alguém sabe como fazer para que o R aceite esses vetores grandes? Obrigado e abraços, Itamar.

2 5

teste de wilcoxon
by alanarocha＠sapo.pt 28 Mai '13

28 Mai '13

para o teste dos dados o melhor gráfico será um plot ou um boxplot? X2sem2011 custosnaopanakeia hgbnaopanakeia 1 JUL 29.41 75.9 2 AGO 32.28 77.5 3 SET 26.88 79.4 4 OUT 25.94 76.8 5 NOV 35.27 76.1 6 DEZ 24.50 78.3 Wilcoxon signed rank test data: as.numeric(dados$custosnaopanakeia) and as.numeric(dados$hgbnaopanakeia) V = 0, p-value = 0.03125 alternative hypothesis: true location shift is not equal to 0 Cumprimentos Ana Rocha

1 0