- R-br - DINF UFPR

Selecionar linhas e colunas
by Edmar Caldas 24 Mar '17

24 Mar '17

Bom Dia! Preciso de mais uma ajuda com seleção de casos. Estou tentando selecionar na coluna Fratura todo mundo que tem a categoria ausente, Porém só quero que essa seleção seja valido para o ID que seja >= 50. raiox[c(raiox$Fratura=="Ausente","id >= 50",] Edmar

3 2

Ajuda no calculo entre duas datas
by Edmar Caldas 24 Mar '17

24 Mar '17

Pessoal, Preciso de uma ajuda em um calculo entre 02 datas. Preciso que o resultado seja numero inteiro que nesse caso é 38 anos. > data_nascimento <- as.Date("1978-06-14") > data_nascimento [1] "1978-06-14" > data_hoje <-as.Date("2017-03-23") > data_hoje [1] "2017-03-23" > idade <-difftime(data_hoje, data_nascimento)/365 > idade Time difference of 38.8 days Edmar

4 3

Salvar e ler de bando de dados
by Adriele Giaretta Biase 23 Mar '17

23 Mar '17

Olá pessoal, eu estou trabalhando no R com conjunto de dados muito grande (numa fase de teste) e a tendência é que essa banco de dados aumente ainda mais, numa próxima etapa que não será teste. E em algumas de minhas programações, o R retornou mensagem de aviso dizendo que não podia alocar o vetor de tamanho de 3.7 GB. Dessa forma, ao invés de concatenar uma matriz em baixo da outra num objeto criado (onde ele me retorna que não foi capaz de armazenar tudo), eu gostaria de ir salvando cada matriz gerada de resultado no diretório local. Minha dúvida, é como criar esse algoritmo de forma que eu consiga salvar dentro de um loop, data.frames com nomes diferentes, sem que ocorra substituição de um data.frame pelo outro. Eu criei um código mínimo (de forma que vocês possam me dar alguma ajuda), e eu possa adaptar para a situação que estou trabalhando. Esse código está substituindo o data.frame a cada rodada dentro do loop. Eu queria que o nome do data.frame alterasse para que guardasse todos os resultados com nomes diferentes, acompanhando a variação do loop (1,2, ..., N): library(MASS) salva_BD_i <- function(N){ for (i in 1:N){ p <- 30 # numero de variáveis a serem geradas (p) ME <- rep(1, p) rho <- 0.5 # correlacao sigma2 <- 1 sigma <- sigma2 * ((1-rho)*diag(p)+rho*matrix(1, p, p)) n <- 10000 Y <- mvrnorm(n, ME, sigma) write.table(Y, "i.txt") } } N=30000 salva = salva_BD_i(N) Numa segunda etapa, eu precisaria de ler (puxar todos os data.frame para dentro de um objeto do R, num loop. Também criei um código mínimo executável, mas acontece o mesmo, o índice do loop não é processado dentro do read.table: > setwd('C:\\Users\\Adriele\\ teste') > dir() [1] "1.txt" "2.txt" "3.txt" "4.txt" "5.txt" "6.txt" "7.txt" "8.txt" "9.txt" "10.txt" "Adriele.R" Leia_BD_i <- function(N){ aux= c() for (i in 1:N){ BD_i <- read.table("i.txt", header = TRUE) media = apply(BD_i, 2, mean) aux = rbind(aux, media) } return(medias = aux) } N=10 Leitura = Leia_BD_i(N) Error in file(file, "rt") : cannot open the connection In addition: Warning message: In file(file, "rt") : cannot open file 'i.txt': No such file or directory Agradeço se alguém puder contribuir, -- Adriele Giaretta Biase. Mestre em Estatística e Experimentação Agropecuária - UFLA. Doutora em Estatística e Experimentação Agronômica - ESALQ/ USP Contato: (19) 98861-0619.

2 1

Verificar multicolinearidade
by Elias Carvalho 23 Mar '17

23 Mar '17

Boa noite pessoal Alguém poderia me indicar os melhores métodos via R para identificar multicolinearidade de dados ? Pelo que pesquisei no forum a maioria das informações parece ser antigas (2011 a 2015) por isso gostaria de algo mais recente. -- *In Jesu et Maria* *Obrigado* *Prof. Elias Carvalho* *"Felix, qui potuit rerum cognoscere causas" (Virgil 29 BC)"Blessed is he who has been able to understand the cause of things"*

2 1

Hel in R: Erro (Permission denied)
by Elvira De Bastiani 23 Mar '17

23 Mar '17

Bom dia pessoal! Sou iniciante no programa R e estou utilizando o pacote Bipartite para minhas análises. Eu estou tentando rodar um modelo com 1.000 amostras de redes de interação que são criadas aleatoriamente a partir de um banco de dados. E agora durante a análise, começou a aparecer o seguinte erro: Error in file(file, ifelse(append, "a", "w")) : cannot open the connection In addition: Warning message: In file(file, ifelse(append, "a", "w")) : cannot open file 'web-1-1.pairs': Permission denied Aparentemente, durante a análise é gerado um arquivo que neste caso é o 'web-1-1.pairs', e o mesmo tem acesso negado, não finalizando a análise para as 1.000 redes de interação. Pelo que andei pesquisando, este erro pode acontecer quando: Tem um arquivo aberto em outro programa, como o Excel, enquanto tenta gravá-lo..ou se o diretório de trabalho não tenha acesso para gravação a esse diretório.. mas acredito que isso não seja meu caso. Alguém conhece ou já teve esse erro, e poderia compartilhar como solucioná-lo ? Muito obrigada Elvira D'Bastiani Mestranda em Ecologia e Conservação Laboratório de Ecologia e Evolução de Interações (Leei) Universidade Federal do Paraná - UFPR Curitiba, Paraná, Brasil

1 0

duvida sobre experimento fatorial
by fúlvia cristina oliveira 21 Mar '17

21 Mar '17

Boa tarde Eu estou tentando rodar os dados de um experimento com dois fatores, os dados estão desbalanceadosEu preciso fazer o desdobramentos desses dados?Outra duvida é sobre os teste de comparação de médias para DIC com dados desbalanceados, no teste de Tukey não houve diferenças para p<0.05, mas no teste de Duncan eu obtive diferença significativa. Qual teste eu devo levar em conta? Posso ser questionada se usar o Duncan? Fúlvia Cristina Oliveira Bióloga - UFMS Mestranda no Programa de Pós-Graduação em ZootecniaProdução animal no Cerrado-Pantanal Universidade Estadual de Mato Grosso do Sul - UEMS Unidade Universitária de Aquidauana - UUA(67) 99610-2812(11) 99538-4088fulcris(a)yahoo.com.br

2 1

Dúvida de como usar o split com o Anova do pacote Car
by Olympio Neto 21 Mar '17

21 Mar '17

Prezados, Quero fazer um desdobramento interação em experimento fatorial muito semelhante ao que foi apresentado na página R-íduculas (Desdobramento de interação em experimento fatorial - http://www.leg.ufpr.br/doku.php/ridiculas) no entanto, precisaria usar a função Anova do pacote CAR que não aceita o argumento "split" para, justamente, calcular os desdobramentos dados os contrastes. Existe alguma forma de usar algo parecido com o argumento split, mas com a função "Anova" do CAR? Muito obrigado Olympio

3 3

Re: [R-br] Box's M Test (Pollyanna Alves) (Pollyanna Alves)
by Maurício Sangiogo 20 Mar '17

20 Mar '17

Boa tarde Pollyana, Tentei rodar aqui com todas suas variáveis e realmente deu erro.... tentei encontrar um motivo lógico e não encontrei, por fim fui por tentativa e erro eliminando variáveis uma a uma e percebi que o problema provavelmente seja o número de variáveis em teste, ao que entendi essa função só consegue testar até 20 variáveis resposta por vez. Perceba que se tu delimitar tuas variáveis em 20 dará certo: require(biotools) dados<-read.xlsx(file.choose(), sheet=1) boxM(dado[,(2:22)],dado[,1]) Nesse comando eu usei as 20 primeiras variáveis (deixei a última de fora no caso), e deu certo. H0 rejeitada... Se precisar de mais algo só falar.. Att .................. Eng. Agr. Maurício Sangiogo TEL:(55) 96822030 ________________________________ De: R-br <r-br-bounces(a)listas.c3sl.ufpr.br> em nome de r-br-request(a)listas.c3sl.ufpr.br <r-br-request(a)listas.c3sl.ufpr.br> Enviado: segunda-feira, 20 de março de 2017 12:00 Para: r-br(a)listas.c3sl.ufpr.br Assunto: Digest R-br, volume 75, assunto 14 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 R-br Página de Detalhes<https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br> listas.inf.ufpr.br A [R-br] é a lista Brasileira oficial de discussão do programa R e tem o propósito de permitir a troca de informações entre os usuários de R (em português). 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: Box's M Test (Pollyanna Alves) (Pollyanna Alves) ---------------------------------------------------------------------- Message: 1 Date: Mon, 20 Mar 2017 09:16:22 -0300 From: Pollyanna Alves <pollyalvesb(a)gmail.com> To: Maurício Sangiogo <ms_sangiogo(a)hotmail.com>, a lista Brasileira oficial de discussão do programa R. <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] Box's M Test (Pollyanna Alves) Message-ID: <CAF=U39x5Xx_AYT51B7aR-4UipY_wDf2=ZoOdgeDHMFgDD3W8AA(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" Olá Maurício, obrigada pela dica, mas quando o apliquei foi gerado o seguinte erro: > boxM(castoria[,-1],castoria[,1]) Error in boxM(castoria[, -1], castoria[, 1]) : there are one or more levels with less observations than variables! Seguem meus dados em anexo, caso possa me ajudar. Obrigada, Pollyanna Alves de Barros Mestranda em Biologia Animal Graduada em Ciências Biológicas Graduação Sanduíche na "University of Western Australia" Museu de Zoologia João Moojen Universidade Federal de Viçosa Em 18 de março de 2017 13:48, Maurício Sangiogo via R-br < r-br(a)listas.c3sl.ufpr.br> escreveu: > Olá Poliana > > > O "grouping" seria teu vetor de tratamentos, no seu caso (se entendi > corretamente), seria a coluna com suas 189 espécies. > > > Se tu for explorar o exemplo dessa função tu verá que nesse caso o vetor > de tratamentos é conseguido por "iris[,5]" (seria a coluna das "species" > dentro do data frame "iris". Ainda sobre o exemplo, "iris[, -5]" seria os > dados com todas as variáveis resposta a serem analisadas (o comando está > indicando tudo "menos" a coluna dos "tratamentos", dentro do data frame > "iris"). > > > > Supondo que em teu caso em especial, o nome do teu data frame com os dados > seja "dados" e a primeira coluna seja das tuas espécies e as demais com as > variáveis. Nesse caso o comando para teus dados seria mais ou menos assim: > > > boxM(dados[, -1], dados[, 1]) > > > Se não conseguir mesmo assim passe parte dos seus dados para que possamos > ajudar melhor. > > > Att, > > *Maurício Sangiogo * > > > > ------------------------------ > *De:* R-br <r-br-bounces(a)listas.c3sl.ufpr.br> em nome de > r-br-request(a)listas.c3sl.ufpr.br <r-br-request(a)listas.c3sl.ufpr.br> > *Enviado:* sábado, 18 de março de 2017 12:00 > *Para:* r-br(a)listas.c3sl.ufpr.br > *Assunto:* Digest R-br, volume 75, assunto 12 > > 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 R-br Página de Detalhes<https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br> listas.inf.ufpr.br A [R-br] é a lista Brasileira oficial de discussão do programa R e tem o propósito de permitir a troca de informações entre os usuários de R (em português). > R-br Página de Detalhes > <https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br> R-br Página de Detalhes<https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br> listas.inf.ufpr.br A [R-br] é a lista Brasileira oficial de discussão do programa R e tem o propósito de permitir a troca de informações entre os usuários de R (em português). > listas.inf.ufpr.br > A [R-br] é a lista Brasileira oficial de discussão do programa R e tem o > propósito de permitir a troca de informações entre os usuários de R (em > português). > > > 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. Box's M Test (Pollyanna Alves) > 2. Re: Dúvida de como usar o split com o Anova do pacote Car > (Fernando Antonio de souza) > 3. Re: grafico de distribuicao acumulada de uma variavel > continua (Cesar Rabak) > 4. Ajuda comando restruturação (Edmar Caldas) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Fri, 17 Mar 2017 16:54:38 -0300 > From: Pollyanna Alves <pollyalvesb(a)gmail.com> > To: a lista Brasileira oficial de discussão do programa R. > <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Box's M Test > Message-ID: > <CAF=U39wgQ5m=A22L-_yfJqS4KRv+Gk=6w8SMttvHe-OxcApCiw@mail. > gmail.com> > Content-Type: text/plain; charset="utf-8" > > Oi pessoal, > > estou tentando realizar o teste Box's M que avalia a homogeneidade da > matriz var-covar dos dados, pressuposto para MANOVA, mas não estou > entendendo o que seria esse *grouping *no script: > > *Pacote *biotools > > *Usage* > boxM(data, grouping) > > *Arguments* > data a numeric data.frame or matrix containing n > observations of p variables; it is expected > that n > p. > > grouping a vector of length n containing the class of each > observation; it is usualy a factor. > > > *O exemplo fornecido foi esse:* > > data(iris) > boxM(iris[, -5], iris[, 5]) > # End (not run) > > Os meus dados possuem 23 variáveis (colunas: medidas) e 189 linhas > (espécimes) sem contar com o cabeçário. > > Alguém poderia me ajudar? > > Agradeço desde já, > > Att., > > > Pollyanna Alves de Barros > Mestranda em Biologia Animal > Graduada em Ciências Biológicas > Graduação Sanduíche na "University of Western Australia" > Museu de Zoologia João Moojen > Universidade Federal de Viçosa > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/ > 20170317/f3284740/attachment-0001.html> > > ------------------------------ > > Message: 2 > Date: Fri, 17 Mar 2017 19:07:03 -0300 > From: Fernando Antonio de souza <nandodesouza(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br, Cesar Rabak <cesar.rabak(a)gmail.com> > Subject: Re: [R-br] Dúvida de como usar o split com o Anova do pacote > Car > Message-ID: > <CAFrWFs=nNtO1xcZyi=Wb3WUsgS9_vwD1FwpogvaW29J6+5HAiQ@mail. > gmail.com> > Content-Type: text/plain; charset="utf-8" > > Veja o pacote phia ele oferece um meio mais fácil de desdobrar interações. > O caminho do "split" é embaraçoso e complicado. > > Em 13/03/2017 16:43, "Cesar Rabak via R-br" <r-br(a)listas.c3sl.ufpr.br> > escreveu: > > > Olympio, > > > > O que diz a documentação da função car::Anova()? > > > > > > 2017-03-13 15:17 GMT-03:00 Olympio Neto via R-br < > r-br(a)listas.c3sl.ufpr.br > > >: > > > >> Prezados, > >> > >> Quero fazer um desdobramento interação em experimento fatorial muito > >> semelhante ao que foi apresentado na página R-íduculas (Desdobramento de > >> interação em experimento fatorial - http://www.leg.ufpr.br/doku.ph > >> p/ridiculas), no entanto, precisaria usar a função Anova do pacote CAR > >> que não aceita o argumento "split" para, justamente, calcular os > >> desdobramentos dados os contrastes. > >> Existe alguma forma de usar algo parecido com o argumento split, mas > com > >> a função "Anova" do CAR? > >> > >> Muito obrigado > >> > >> Olympio > >> > >> > >> _______________________________________________ > >> R-br mailing list > >> R-br(a)listas.c3sl.ufpr.br > >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > R-br Página de Detalhes > <https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br> > listas.inf.ufpr.br > A [R-br] é a lista Brasileira oficial de discussão do programa R e tem o > propósito de permitir a troca de informações entre os usuários de R (em > português). > > > >> 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 > R-br Página de Detalhes > <https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br> > listas.inf.ufpr.br > A [R-br] é a lista Brasileira oficial de discussão do programa R e tem o > propósito de permitir a troca de informações entre os usuários de R (em > português). > > > > 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/ > 20170317/41a07cd9/attachment-0001.html> > > ------------------------------ > > Message: 3 > Date: Sat, 18 Mar 2017 01:51:51 -0300 > From: Cesar Rabak <cesar.rabak(a)gmail.com> > To: Nilson Guiotoku <nguiotoku2017(a)gmail.com>, a lista Brasileira > oficial de discussão do programa R. <r-br(a)listas.c3sl.ufpr.br> > Subject: Re: [R-br] grafico de distribuicao acumulada de uma variavel > continua > Message-ID: > <CAKrF98nFbZY=_EmUA7yXLX6edTc5bCbCsdywYOwyY7tPhyYT4g@mail. > gmail.com> > Content-Type: text/plain; charset="utf-8" > > http://bfy.tw/Ahi3 > <http://bfy.tw/Ahi3> > LMGTFY <http://bfy.tw/Ahi3> > bfy.tw > LMGTFY > > > > 2017-03-16 14:49 GMT-03:00 Nilson Guiotoku via R-br < > r-br(a)listas.c3sl.ufpr.br>: > > > Bom dia! > > > > > > Eu preciso construir o gráfico da densidade acumulada de uma variável > > continua (domínio dos Reais) que não tem um modelo identificado. Usei o > > gráfico de pareto construído por tabelas de frequência, porém, no eixo x, > > ele coloca as frequência em ordem decrescente. Eu gostaria de construir > > esse mesmo gráfico, todavia, em que no eixo x ficasse as classes das > > tabelas na sequencia crescente (Ex: (-17.1,32.9] (32.9,82.9] > > (82.9,133] (133,183] (183,233] (233,283] (283,333] ), não > > usando a sequencia decrescente da frequência como no pareto. Alguma ideia > > de como fazer isso? > > > > > # Grafico de pareto > > > > > MT <- Global$Z > > > > > head(MT) > > > > [1] 127.45555 23.10762 36.16764 88.45139 142.60974 151.68347 > > > > > > > > > intervalo <- cut(MT, breaks= seq(min(MT), max(MT), 50 )) > > > > > > > > > x = table(intervalo) > > > > > print(x) > > > > intervalo > > > > (-17.1,32.9] (32.9,82.9] (82.9,133] (133,183] (183,233] > > (233,283] (283,333] > > > > 5 6 9 > > 7 4 5 6 > > > > > > > > > require(qcc) > > > > > pareto.chart(x,ylab = "Frequencia", main = "classes do lucro MT", > > > > + col=rainbow(length(x))) > > > > > > > > Pareto chart analysis for x > > > > Frequency Cum.Freq. Percentage Cum.Percent. > > > > (82.9,133] 9 9 21.42857 21.42857 > > > > (133,183] 7 16 16.66667 38.09524 > > > > (32.9,82.9] 6 22 14.28571 52.38095 > > > > (283,333] 6 28 14.28571 66.66667 > > > > (-17.1,32.9] 5 33 11.90476 78.57143 > > > > (233,283] 5 38 11.90476 90.47619 > > > > (183,233] 4 42 9.52381 100.00000 > > > > > > > > > > > > > > > > Ver Figura 1 do pareto em anexo. > > > > No eixo x (aparece as classes fora de ordem), mas na ordem decrescente da > > frequência. Gostaria que as classes ficassem na ordem no eixo x > > ((-17.1,32.9] (32.9,82.9] (82.9,133] (133,183] (183,233] > > (233,283] (283,333] ) > > > > > > > > O ideal é que apareça na ordem crescente no eixo x, como na Figura em > > anexo (Figura 2). > > > > > > > > Nilson Silva. > > > > _______________________________________________ > > R-br mailing list > > R-br(a)listas.c3sl.ufpr.br > > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > R-br Página de Detalhes > <https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br> > listas.inf.ufpr.br > A [R-br] é a lista Brasileira oficial de discussão do programa R e tem o > propósito de permitir a troca de informações entre os usuários de R (em > português). > > > > 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/ > 20170318/ae93d34b/attachment-0001.html> > > ------------------------------ > > Message: 4 > Date: Sat, 18 Mar 2017 13:44:23 +0000 (UTC) > From: Edmar Caldas <edmar_caldas(a)yahoo.com.br> > To: R-br Lista <r-br(a)listas.c3sl.ufpr.br> > Subject: [R-br] Ajuda comando restruturação > Message-ID: <189887133.2734358.1489844663180(a)mail.yahoo.com> > Content-Type: text/plain; charset="utf-8" > > > Bom Dia! > Preciso de uma ajuda e não sei como começar. > tenho uma arquivo que está na coluna e quero que ele fique na linha. segue > exemplo. Estou procurando um comando de restruturação e não encontro. > alguém pode me direcionar. > > > > > Edmar > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/ > 20170318/b251304d/attachment-0001.html> > -------------- Próxima Parte ---------- > Um anexo não-texto foi limpo... > Nome: datauri-file.png > Tipo: image/png > Tamanho: 20423 bytes > Descrição: não disponível > URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/ > 20170318/b251304d/attachment-0001.png> > > ------------------------------ > > 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 > R-br Página de Detalhes > <https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br> > listas.inf.ufpr.br > A [R-br] é a lista Brasileira oficial de discussão do programa R e tem o > propósito de permitir a troca de informações entre os usuários de R (em > português). > > > > > ------------------------------ > > Fim da Digest R-br, volume 75, assunto 12 > ***************************************** > > > _______________________________________________ > 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/20170320/6c08195e/atta…> -------------- Próxima Parte ---------- Um anexo não-texto foi limpo... Nome: Castoria Dados Estimados 170317.csv Tipo: text/csv Tamanho: 24765 bytes Descrição: não disponível URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20170320/6c08195e/atta…> ------------------------------ 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 75, assunto 14 *****************************************

1 0

Re: [R-br] Box's M Test (Pollyanna Alves)
by Maurício Sangiogo 20 Mar '17

20 Mar '17

Olá Poliana O "grouping" seria teu vetor de tratamentos, no seu caso (se entendi corretamente), seria a coluna com suas 189 espécies. Se tu for explorar o exemplo dessa função tu verá que nesse caso o vetor de tratamentos é conseguido por "iris[,5]" (seria a coluna das "species" dentro do data frame "iris". Ainda sobre o exemplo, "iris[, -5]" seria os dados com todas as variáveis resposta a serem analisadas (o comando está indicando tudo "menos" a coluna dos "tratamentos", dentro do data frame "iris"). Supondo que em teu caso em especial, o nome do teu data frame com os dados seja "dados" e a primeira coluna seja das tuas espécies e as demais com as variáveis. Nesse caso o comando para teus dados seria mais ou menos assim: boxM(dados[, -1], dados[, 1]) Se não conseguir mesmo assim passe parte dos seus dados para que possamos ajudar melhor. Att, Maurício Sangiogo ________________________________ De: R-br <r-br-bounces(a)listas.c3sl.ufpr.br> em nome de r-br-request(a)listas.c3sl.ufpr.br <r-br-request(a)listas.c3sl.ufpr.br> Enviado: sábado, 18 de março de 2017 12:00 Para: r-br(a)listas.c3sl.ufpr.br Assunto: Digest R-br, volume 75, assunto 12 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 R-br Página de Detalhes<https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br> listas.inf.ufpr.br A [R-br] é a lista Brasileira oficial de discussão do programa R e tem o propósito de permitir a troca de informações entre os usuários de R (em português). 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. Box's M Test (Pollyanna Alves) 2. Re: Dúvida de como usar o split com o Anova do pacote Car (Fernando Antonio de souza) 3. Re: grafico de distribuicao acumulada de uma variavel continua (Cesar Rabak) 4. Ajuda comando restruturação (Edmar Caldas) ---------------------------------------------------------------------- Message: 1 Date: Fri, 17 Mar 2017 16:54:38 -0300 From: Pollyanna Alves <pollyalvesb(a)gmail.com> To: a lista Brasileira oficial de discussão do programa R. <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] Box's M Test Message-ID: <CAF=U39wgQ5m=A22L-_yfJqS4KRv+Gk=6w8SMttvHe-OxcApCiw(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" Oi pessoal, estou tentando realizar o teste Box's M que avalia a homogeneidade da matriz var-covar dos dados, pressuposto para MANOVA, mas não estou entendendo o que seria esse *grouping *no script: *Pacote *biotools *Usage* boxM(data, grouping) *Arguments* data a numeric data.frame or matrix containing n observations of p variables; it is expected that n > p. grouping a vector of length n containing the class of each observation; it is usualy a factor. *O exemplo fornecido foi esse:* data(iris) boxM(iris[, -5], iris[, 5]) # End (not run) Os meus dados possuem 23 variáveis (colunas: medidas) e 189 linhas (espécimes) sem contar com o cabeçário. Alguém poderia me ajudar? Agradeço desde já, Att., Pollyanna Alves de Barros Mestranda em Biologia Animal Graduada em Ciências Biológicas Graduação Sanduíche na "University of Western Australia" Museu de Zoologia João Moojen Universidade Federal de Viçosa -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20170317/f3284740/atta…> ------------------------------ Message: 2 Date: Fri, 17 Mar 2017 19:07:03 -0300 From: Fernando Antonio de souza <nandodesouza(a)gmail.com> To: r-br(a)listas.c3sl.ufpr.br, Cesar Rabak <cesar.rabak(a)gmail.com> Subject: Re: [R-br] Dúvida de como usar o split com o Anova do pacote Car Message-ID: <CAFrWFs=nNtO1xcZyi=Wb3WUsgS9_vwD1FwpogvaW29J6+5HAiQ(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" Veja o pacote phia ele oferece um meio mais fácil de desdobrar interações. O caminho do "split" é embaraçoso e complicado. Em 13/03/2017 16:43, "Cesar Rabak via R-br" <r-br(a)listas.c3sl.ufpr.br> escreveu: > Olympio, > > O que diz a documentação da função car::Anova()? > > > 2017-03-13 15:17 GMT-03:00 Olympio Neto via R-br <r-br(a)listas.c3sl.ufpr.br > >: > >> Prezados, >> >> Quero fazer um desdobramento interação em experimento fatorial muito >> semelhante ao que foi apresentado na página R-íduculas (Desdobramento de >> interação em experimento fatorial - http://www.leg.ufpr.br/doku.ph >> p/ridiculas), no entanto, precisaria usar a função Anova do pacote CAR >> que não aceita o argumento "split" para, justamente, calcular os >> desdobramentos dados os contrastes. >> Existe alguma forma de usar algo parecido com o argumento split, mas com >> a função "Anova" do CAR? >> >> Muito obrigado >> >> Olympio >> >> >> _______________________________________________ >> R-br mailing list >> R-br(a)listas.c3sl.ufpr.br >> https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br R-br Página de Detalhes<https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br> listas.inf.ufpr.br A [R-br] é a lista Brasileira oficial de discussão do programa R e tem o propósito de permitir a troca de informações entre os usuários de R (em português). >> 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 R-br Página de Detalhes<https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br> listas.inf.ufpr.br A [R-br] é a lista Brasileira oficial de discussão do programa R e tem o propósito de permitir a troca de informações entre os usuários de R (em português). > 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/20170317/41a07cd9/atta…> ------------------------------ Message: 3 Date: Sat, 18 Mar 2017 01:51:51 -0300 From: Cesar Rabak <cesar.rabak(a)gmail.com> To: Nilson Guiotoku <nguiotoku2017(a)gmail.com>, a lista Brasileira oficial de discussão do programa R. <r-br(a)listas.c3sl.ufpr.br> Subject: Re: [R-br] grafico de distribuicao acumulada de uma variavel continua Message-ID: <CAKrF98nFbZY=_EmUA7yXLX6edTc5bCbCsdywYOwyY7tPhyYT4g(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" http://bfy.tw/Ahi3 [http://lmgtfy.com/assets/sticker-b222a421fb6cf257985abfab188be7d6746866850e…]<http://bfy.tw/Ahi3> LMGTFY<http://bfy.tw/Ahi3> bfy.tw LMGTFY 2017-03-16 14:49 GMT-03:00 Nilson Guiotoku via R-br < r-br(a)listas.c3sl.ufpr.br>: > Bom dia! > > > Eu preciso construir o gráfico da densidade acumulada de uma variável > continua (domínio dos Reais) que não tem um modelo identificado. Usei o > gráfico de pareto construído por tabelas de frequência, porém, no eixo x, > ele coloca as frequência em ordem decrescente. Eu gostaria de construir > esse mesmo gráfico, todavia, em que no eixo x ficasse as classes das > tabelas na sequencia crescente (Ex: (-17.1,32.9] (32.9,82.9] > (82.9,133] (133,183] (183,233] (233,283] (283,333] ), não > usando a sequencia decrescente da frequência como no pareto. Alguma ideia > de como fazer isso? > > > # Grafico de pareto > > > MT <- Global$Z > > > head(MT) > > [1] 127.45555 23.10762 36.16764 88.45139 142.60974 151.68347 > > > > > intervalo <- cut(MT, breaks= seq(min(MT), max(MT), 50 )) > > > > > x = table(intervalo) > > > print(x) > > intervalo > > (-17.1,32.9] (32.9,82.9] (82.9,133] (133,183] (183,233] > (233,283] (283,333] > > 5 6 9 > 7 4 5 6 > > > > > require(qcc) > > > pareto.chart(x,ylab = "Frequencia", main = "classes do lucro MT", > > + col=rainbow(length(x))) > > > > Pareto chart analysis for x > > Frequency Cum.Freq. Percentage Cum.Percent. > > (82.9,133] 9 9 21.42857 21.42857 > > (133,183] 7 16 16.66667 38.09524 > > (32.9,82.9] 6 22 14.28571 52.38095 > > (283,333] 6 28 14.28571 66.66667 > > (-17.1,32.9] 5 33 11.90476 78.57143 > > (233,283] 5 38 11.90476 90.47619 > > (183,233] 4 42 9.52381 100.00000 > > > > > > > > Ver Figura 1 do pareto em anexo. > > No eixo x (aparece as classes fora de ordem), mas na ordem decrescente da > frequência. Gostaria que as classes ficassem na ordem no eixo x > ((-17.1,32.9] (32.9,82.9] (82.9,133] (133,183] (183,233] > (233,283] (283,333] ) > > > > O ideal é que apareça na ordem crescente no eixo x, como na Figura em > anexo (Figura 2). > > > > Nilson Silva. > > _______________________________________________ > R-br mailing list > R-br(a)listas.c3sl.ufpr.br > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br R-br Página de Detalhes<https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br> listas.inf.ufpr.br A [R-br] é a lista Brasileira oficial de discussão do programa R e tem o propósito de permitir a troca de informações entre os usuários de R (em português). > 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/20170318/ae93d34b/atta…> ------------------------------ Message: 4 Date: Sat, 18 Mar 2017 13:44:23 +0000 (UTC) From: Edmar Caldas <edmar_caldas(a)yahoo.com.br> To: R-br Lista <r-br(a)listas.c3sl.ufpr.br> Subject: [R-br] Ajuda comando restruturação Message-ID: <189887133.2734358.1489844663180(a)mail.yahoo.com> Content-Type: text/plain; charset="utf-8" Bom Dia! Preciso de uma ajuda e não sei como começar. tenho uma arquivo que está na coluna e quero que ele fique na linha. segue exemplo. Estou procurando um comando de restruturação e não encontro. alguém pode me direcionar. Edmar -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20170318/b251304d/atta…> -------------- Próxima Parte ---------- Um anexo não-texto foi limpo... Nome: datauri-file.png Tipo: image/png Tamanho: 20423 bytes Descrição: não disponível URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20170318/b251304d/atta…> ------------------------------ 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 R-br Página de Detalhes<https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br> listas.inf.ufpr.br A [R-br] é a lista Brasileira oficial de discussão do programa R e tem o propósito de permitir a troca de informações entre os usuários de R (em português). ------------------------------ Fim da Digest R-br, volume 75, assunto 12 ***************************************** [http://listas.inf.ufpr.br/pipermail/r-br/attachments/20170318/b251304d/atta…]

2 1

Ajuda comando restruturação
by Edmar Caldas 18 Mar '17

18 Mar '17

Bom Dia! Preciso de uma ajuda e não sei como começar. tenho uma arquivo que está na coluna e quero que ele fique na linha. segue exemplo. Estou procurando um comando de restruturação e não encontro. alguém pode me direcionar. Edmar

2 1