- R-br - DINF UFPR

Maximizar plot window no RStudio
by Elias Carvalho 13 Abr '17

13 Abr '17

Bom dia pessoal Alguém sabe como posso maximizar a janela plot do RSstudio via linha de comando ? Explicando: Estou fazendo uma rotina para verificação de outliers baseado no boxplot então ao examinar os dados eu gero os gráficos abaixo, no entanto, se a janela estiver pequena demais os gráficos não aparecem e ainda surge uma mensagem Error in plot.new() : figure margins too large in R <http://stackoverflow.com/questions/12766166/error-in-plot-new-figure-margin…> : -- *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"*

3 4

Combinar multiplos datasets - Inputação Múltipla com MICE
by Elias Carvalho 11 Abr '17

11 Abr '17

*Boa noite pessoal* *Eu preciso tratar dados missing do meu dataset, li bastante e o MICE me pareceu uma boa ferramente para isso, então p**rossegui da seguinte maneira:* mice(df.factor.variable, m=5, maxit = 50, method = 'polyreg', seed = 500) *e depois para juntar os 5 datasets usei:* complete.factor.data <- complete(factor.variable.imputed, 'long') *Como tenho 1999 registros no meu dataset e foram gerados 5 conjuntos de dados imputados, o resultado foi de 9995 observações (1999 * 5).* *Alguém sabe me dizer se o procedimento está correto ou há uma maneira de unir os 5 datasets imputados em um só gerando um dataset de 1999 registros sem NAs?* *ps: vou usar esse novo dataset sem missing para aprender uma rede Bayesiana* -- *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"*

1 0

código para construir tabela com dados pareados
by Mauricio Cardeal 09 Abr '17

09 Abr '17

Boa tarde! Por favor, como faço para construir uma tabela como a do modelo do final do script, onde o conteúdo das 4 caselas são contagens do número de pares que combinam melhora e não melhora de 2 tratamentos aplicados simultaneamente no mesmo individuo (perna direita e perna esquerda), por exemplo ? # CMR: par <- c(1,1,2,2,3,3,4,4,5,5,6,6,7,7) # refere-se ao número do individuo trata <- c(1,0,1,0,1,0,1,0,1,0,1,0,1,0) # tratamento 1 (A) ou 0 (B) melhora <- c(1,0,1,0,1,0,0,1,0,1,0,0,1,1) # resultado do tratamento 1 (melhora) 0 (não melhora) tratac <- factor(trata) levels(tratac)<-list("A"=1,"B"=0 ) melhorac <- factor(melhora) levels(melhorac)<-list("+"=1,"-"=0) data.frame(par,tratac,melhorac) par tratac melhorac 1 1 A + 2 1 B - 3 2 A + 4 2 B - 5 3 A + 6 3 B - 7 4 A - 8 4 B + 9 5 A - 10 5 B + 11 6 A - 12 6 B - 13 7 A + 14 7 B + # COMO FAÇO PARA GERAR ESSA TABELA A PARTIR DOS DADOS ACIMA? # tratamento B # melhora melhora # + - #tratamento A melhora + 1 3 # melhora - 2 1 Obrigado, Maurício Cardeal UFBA

2 3

Suavização do gráfico
by Mateus Dias Nunes 08 Abr '17

08 Abr '17

Bom dia, Eu estou com um problema para suavizar meu gráfico. Eu tenho dados horários e desejo plotar individualmente cada data (no período das 06h as 18h). Entretanto o gráfico tem ficado com muito ruído devido as observações serem horárias e existirem falhas (NA) nas observações. E por isso eu gostaria de suavizar o gráfico com as médias horárias, porém não estou conseguindo implementar isso. Espero que alguém possa me ajudar a resolver. dado utilizado: https://www.dropbox.com/s/b80zzzeaesnh8qp/dado.txt?dl=0 library(openair) mydata<-read.table("dado.txt", h = T, sep = ",") mydata$direta <- as.numeric(as.character(mydata$direta)) mydata$date <- as.POSIXct(strptime(mydata$date, format = "%d/%m/%Y %H:%M:%S",tz = "Etc/GMT-1")) means <- aggregate(mydata["direta"], format(mydata["date"],"%d-%m-%Y %H:%M:%S"), mean, na.rm = TRUE) means[ ,"date"] <- seq(min(mydata[ , "date"]), max(mydata[ , "date"]), length = nrow(means)) Figura1<-plot(mydata$date[2:721], mydata$direta[2:721], type = "l", col="red",xlab = "minutes",main = "01 de Janeiro de 2007", ylab = "Direct Solar Radiation (W/m²)") Figura1 dev.copy(device=jpeg,filename="Brasilia_01dejaneiro2007.jpg", width = 800, height = 500) dev.off () desde já, agradeço. _______________________________________________________________________________________ MATEUS DIAS NUNES Postgraduate Program in Meteorology - PGMET - CPTEC / INPE Center for Weather Forecasting and Climate Studies / National Institute for Space Research Cachoeira Paulista - SP - Brasil FONE: +55 (53) 98112-5154 SKYPE: omsndomateus(a)hotmail.com ____________________________________________ ___________________________________________ <https://mailtrack.io/> Sent with Mailtrack <https://mailtrack.io/install?source=signature&lang=en&referral=nunes.mateus…> <#>

2 1

Duvida em estatística no R
by Tiago G. Santos 08 Abr '17

08 Abr '17

Boa tarde, Alguém sabe como faço para apresentar os resultados dos testes tukey ou dunnet como as letras que indicam significância estatística ..... -- Tiago Gonçalves Santos http://lattes.cnpq.br/7645251131511453 Doutorando em Fìsica pelo Instituto Tecnológico de Aeronáutica - ITA Laboratório de Caracterização de Dispositivos Semicondutores Instituto de Estudos Avançados - IEAv Trevo Coronel Aviador José Alberto Albano do Amarante, no. 1 12.228-001 São José dos Campos – SP – Brasil

2 1

Retornando de NaN para raiz cúbica de números negativos alocados em vetores
by marcos paulo 07 Abr '17

07 Abr '17

Olá prezados boa tarde. Alguém saberia me explicar porque ao concatenar números negativos e positivos em um vetor não é possível obter a raiz cúbica dos números negativos? Como procederia para contornar esse problema, mantendo os dados no vetor ? Exemplo: > vec<-c(27,-27,27,27,-27) > vec^(1/3) [1] 3 NaN 3 3 NaN > -27^(1/3) [1] -3 Situação desejada: [1] 3 -3 3 3 -3 À disposição Marcos Paulo Técnico Agrícola em Zootecnia Bacharel em Agronomia Mestrando em Produção Vegetal - Universidade Federal de Goiás Estagiário de Pós-Graduação da Embrapa Arroz e Feijão/CNPAF-Santo Antônio de Goiás-GO Contato: (62) 85075783 http://lattes.cnpq.br/4322347592884852

4 3

Ajuda: If else
by Edmar Caldas 07 Abr '17

07 Abr '17

Pessoal estou com outra dúvida, preciso criar uma coluna de inconsistência que funciona assim. Zero quando a condição está ok e 1 quando tem erro de inconsistência.minha logica é simples. se estcivil == solteiro e felicidade_casamento == NRA coloca Zero caso contrario coloque 1. erro <-1if (estcivil == "solteiro"){ + else if (Felicidade_casamento != "NRA")else erro <- 0 Edmar

3 2

APSIM.R
by Simone toni ruiz corrêa 07 Abr '17

07 Abr '17

Bom dia a todos, Gostaria de saber se alguém já trabalhou com o packege APSIM.R Trabalho com o modelo APSIM e vi que existe essa ferramenta no R, para auxilio de calibração e análise de sensibilidade. Gostaria de tirar dúvidas... Obrigada, Simone

1 0

- Pacote bipartite
by Elvira De Bastiani 06 Abr '17

06 Abr '17

Boa noite 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 o 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 (41) 9 9653-2860

2 2

Ajuda: Append
by Edmar Caldas 04 Abr '17

04 Abr '17

Pessoal estou com outra duvida. Estou tentando unir dois banco com a função append. tenho arquivo A que contém a variável id, nome, lazer e arquivo B que contém a variável id, nome, trabalho. quero unir os 02 arquivos e ficar com id, nome, lazer, trabalho mesmo que alguns casos fiquem NULL. append <- rbind(LAZER_1,LAZER_2) dá erro de numero de colunas ou argumentos não coincidem. estou tentando fazer isso: LAZER_1$lazer <-NULLLAZER_1 LAZER_2$trabalho <-NULLLAZER_2 depois append <- rbind(LAZER_1,LAZER_2) continua com o mesmo erro? Edmar

4 4

Ajuda: arrendondando uma variável
by Edmar Caldas 03 Abr '17

03 Abr '17

Boa Tarde! Tenho uma dúvida, nesse código a variável idade não é pra ter dígitos. como faço pra tirar isso? # Chamando o arquivocadastro_funcionario data_atual <- as.Date ("2017-04-03")data_atual cadastro_funcionario$data_atual <- c(data_atual)cadastro_funcionario idade <-difftime(data_atual,datanasc)/365.25round(idade, digits = 0) cadastro_funcionario$idade <-c(idade) resultado Edmar

2 1

ENC: Parcela dividida - coletas de amostras dentro da parcela e no tempo
by Fabiola Santos 03 Abr '17

03 Abr '17

Assunto: Parcela dividida - coletas de amostras dentro da parcela e no tempo Boa tarde Tenho um experimento da seguinte forma: 3 tratamentos (tipo de manejo diferente) com 4 repeticoes em bloco. Parcela de 1ha As avaliações foram feitas em subparcelas onde escolhia 5 pontos aleatorios e coletava informações dentro dessas parcelas(de 1ha) em 10 plantas avaliando o numero de lagartas e especie. Porem foram feitas semanalmente as avaliações sendo em varios estagios do milho Gostarida de saber se o modelo a seguir estaria correto: O que eu nao estou conseguindo e fazer os teste de Tukey ou SK no proprio R. mod1 <- aov(PeqS~Rep+Trat*Est + Ponto + Error(Rep:Est, Rep:Ponto), data=x) Error: Rep:Est Df Sum Sq Mean Sq F value Pr(>F) Rep 3 8.010e-27 2.669e-27 0.746 0.541 Est 5 1.815e-26 3.630e-27 1.015 0.443 Residuals 15 5.365e-26 3.577e-27 Error: Within Df Sum Sq Mean Sq F value Pr(>F) Trat 2 0 0.0 0 1 Ponto 1 648 648.0 10689 <2e-16 *** Trat:Est 10 0 0.0 0 1 Residuals 3563 216 0.1 Obrigada Abraços Ótimo final de semana Vídeo Institucional https://www.youtube.com/watch?v=u1Z3_YZ2WpY https://www.facebook.com/udipesquisaedesenvolvimento/ https://www.linkedin.com/company/10401615?trk=tyah <https://www.linkedin.com/company/10401615?trk=tyah&trkInfo=clickedVertical% 3Acompany%2CclickedEntityId%3A10401615%2Cidx%3A2-1-2%2CtarId%3A1479812006594 %2Ctas%3AUDI%20PESQUISA%20E%20DESENV> &trkInfo=clickedVertical%3Acompany%2CclickedEntityId%3A10401615%2Cidx%3A2-1- 2%2CtarId%3A1479812006594%2Ctas%3AUDI%20PESQUISA%20E%20DESENV Atenciosamente / Best Regards / Mit freundlichen Grüßen / Saludos www.udipesquisa.om <http://www.udipesquisa.om> Esta mensagem, incluindo anexos, é confidencial, privilegiada e para o único benefício da pessoa devida ou endereçada. Caso seja recebido erroneamente, favor apagá-la do sistema, não faça cópias, não exponha seu conteúdo a uma terceira parte ou utilize-a para seu benefício pessoal ou de outra pessoa. Favor avisar (fabiolasantos(a)udipesquisa.com <mailto:fabiolasantos@udipesquisa.com> ) do seu recebimento, o mais rápido possível. Agradecemos a sua colaboração. The information contained in this e-mail and in the attached files are for the exclusive use of the addressee herein nominated, and may contain trade secrets, privileged and other confidential information, protected by the applicable laws. In case you are not the right addressee, you are hereby notified that any reviewing, reading, copying and/or distributing of this e-mail's content is strictly prohibited and unauthorized. Please, delete the e-mail's content and notify the sender (fabiolasantos(a)udipesquisa.com <mailto:fabiolasantos@udipesquisa.com> ) immediately. Thank you for your cooperation. --- Este email foi escaneado pelo Avast antivírus. https://www.avast.com/antivirus

1 0

Problemas com NA's
by Mateus Dias Nunes 03 Abr '17

03 Abr '17

Olá Estou com um problema para abrir meus dados. Segui o procediemnto correto para a Abertura dos mesmo, porém um erro insiste em aparecer. library(openair) mydata<-read.csv("abrir.csv",header=T) mydata$date <- as.POSIXct(strptime(mydata$date, format = "%d/%m/%Y %H:%M", tz = "GMT")) mean(mydata$direta, na.rm = TRUE) Warning message: In mean.default(mydata$direta, na.rm = TRUE) : argumento não é numérico nem lógico: retornando NA Os dados faltantes estão em branco no meu vetor de dados, isso que causa estranhesa. abaixo o dado a ser trabalhado. Desde já, agradeço. dados https://www.dropbox.com/s/7p5pw9wk0mz5lmq/abrir.csv?dl=0 _______________________________________________________________________________________ MATEUS DIAS NUNES Postgraduate Program in Meteorology - PGMET - CPTEC / INPE Center for Weather Forecasting and Climate Studies / National Institute for Space Research Cachoeira Paulista - SP - Brasil FONE: +55 (53) 98112-5154 SKYPE: omsndomateus(a)hotmail.com ____________________________________________ ___________________________________________ <https://mailtrack.io/> Sent with Mailtrack <https://mailtrack.io/install?source=signature&lang=en&referral=nunes.mateus…>

3 2

[R-Br] Alguém sabe sobre o suporte a vetores longos (> 2 GB)?
by Maurício Camargo 03 Abr '17

03 Abr '17

Olá a todos. Preciso alocar um vetor com mais de 2 GB. Alguém sabe como anda o suporte a LongVectors do R? Obrigado. Maurício Camargo

3 3

Data frame não reconhece uma variavel
by Edmar Caldas 31 Mar '17

31 Mar '17

Pessoal outra duvida. Criei uma variável data_hoje no meu data.frame - segue comando: Banco$data_hoje <-c(Sys.Date()) funcionou. Está na ultima linha da tabela. Agora quero a diferença entre: data_hoje - datanasc. está com esse erro. data_hoje not found ???????????? Edmar

4 6

Calibratio Plot
by Ary Serpa Neto 31 Mar '17

31 Mar '17

3 2

Agregate
by Edmar Caldas 31 Mar '17

31 Mar '17

Pessoal ainda não consegui. Preciso colocar nas colunas a primeira e ultima data no data.frame. só funciona com o valor_int. Edmar

2 1

Distância euclidiana, centroide
by Marcus Diogo 30 Mar '17

30 Mar '17

Boa noite, Vcs sabem onde posso achar exemplos de implementacoes de Distância euclidiana, Classificador Naive Bayes, Classificador de máxima verossimilhança. É uma tabela do uci repositório.300x7 Agradeço Marcus

1 0

Análise Exploratória dos Dados da Pesquisa Origem e Destino 2007 na Região Metropolitana de São Paulo utilizando a Linguagem R
by Vagner Sanches Vasconcelos 30 Mar '17

30 Mar '17

Projeto realizado para a curso Analista de dados da Udacity https://www.linkedin.com/pulse/an%C3%A1lise-explorat%C3%B3ria-dos-dados-da-… Abraços,

1 0

GLM
by ana paula coelho madeira 29 Mar '17

29 Mar '17

Prezados colegas, Bom dia! Gostaria de saber como faço para extrair a estatística de Qui quadrado e respectivos graus de liberdade de um modelo de GLM significativo. Desde já agradeço. Ana Paula

3 3

Re: [R-br] Predição espacial de dado tipo "linha"
by Elias T. Krainski 28 Mar '17

28 Mar '17

Você quer 'espacializar' a atenuação do sinal? Essa atenuação no sinal seria inversamente proporcional à distância à torre emitente?

1 0

Predição espacial de dado tipo "linha"
by Wagner Wolff 28 Mar '17

28 Mar '17

Olá pessoal da lista! Pretendo trabalhar com um tipo de dado que representa a atenuação de microondas entre duas antenas celulares. Sendo assim, estes dados representam uma "linha" (link) no espaço e não um ponto, como geralmente são feitas as análises geoestatśticas. Gostaria de saber o que muda na análise em comparação com a tradicional por pontos. Por exemplo, ao invés de usar a distância euclidiana entre as amostrar usaria uma distância ponderada pelo comprimento dessas "linhas"? Att -- *Wagner Wolff, **PhD* "*Luiz de Queiroz**" College of Agriculture,* University of São Paulo Pádua Dias avenue11 | 13418-900| Piracicaba-SP| Brazil Phone: +55 19 982385582 http://orcid.org/0000-0003-3426-308X https://github.com/wwolff7 http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4463141A1

2 1

Função cut
by Edmar Caldas 28 Mar '17

28 Mar '17

Bom Dia! Preciso de outra ajuda, preciso quebrar a variável idade em faixas. Estou usando essa função, porém está retornando mensagem de errro. cut(lazer1$idade, breaks = c(25,35,45), labels=c("A","B","C"), right= FALSE)) Error: unexpected ')' in "cut(lazer1$idade, breaks = c(25,35,45), labels=c("A","B","C"), right= FALSE))" outra forma > cut(idade, breaks = c(25,35,45), labels=c("Jovem","Meia Idade","Experiente"), right= FALSE))Error: unexpected ')' in "cut(idade, breaks = c(25,35,45), labels=c("Jovem","Meia Idade","Experiente"), right= FALSE))" Edmar

4 3

tabela de interação
by fúlvia cristina oliveira 27 Mar '17

27 Mar '17

Alguém pode me dizer como eu monto a tabela de desdobramento de interação?A saída foi essa modelo3<-aov(burst~tratamento+tempo+tratamento*tempo)> summary(modelo3) Df Sum Sq Mean Sq F value Pr(>F) tratamento 3 0.0020 0.00066 0.054 0.9833 tempo 2 0.3517 0.17586 14.386 4.81e-06 *** tratamento:tempo 6 0.1968 0.03280 2.683 0.0202 * Residuals 78 0.9535 0.01222 --- Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1 20 observations deleted due to missingness medias<-model.tables(model4,"means",cterms = "tratamento") > medias Tables of means Grand mean 0.3885556 tratamento T1 T2 T3 T4 0.3845 0.3881 0.3965 0.3852 rep 22.0000 24.0000 22.0000 22.0000 model5<-aov(burst~tratamento+tempo+tratamento:tempo) > medias <-model.tables(model5,"means",cterms = "tempo") > medias Tables of means Grand mean 0.3885556 tempo 0 3 6 0.4316 0.4248 0.2907 rep 33.0000 31.0000 26.0000 model6<-aov(burst~tratamento+tempo+tratamento:tempo) > medias<-model.tables(model6,"means",cterms = "tratamento:tempo") > medias Tables of means Grand mean 0.3885556 tratamento:tempo tempo tratamento 0 3 6 T1 0.432 0.422 0.294 rep 7.000 8.000 7.000 T2 0.410 0.403 0.341 rep 10.000 7.000 7.000 T3 0.518 0.370 0.270 rep 8.000 8.000 6.000 T4 0.372 0.501 0.249 rep 8.000 8.000 6.000 TukeyHSD(model4,"tratamento") Tukey multiple comparisons of means 95% family-wise confidence level Fit: aov(formula = burst ~ tratamento + tempo + tratamento:tempo) $tratamento diff lwr upr p adj T2-T1 0.0035833333 -0.08209100 0.08925766 0.9995222 T3-T1 0.0119545455 -0.07556246 0.09947155 0.9840895 T4-T1 0.0007272727 -0.08678973 0.08824427 0.9999962 T3-T2 0.0083712121 -0.07730312 0.09404554 0.9940338 T4-T2 -0.0028560606 -0.08853039 0.08281827 0.9997577 T4-T3 -0.0112272727 -0.09874427 0.07628973 0.9867430 TukeyHSD(model5,"tempo") Tukey multiple comparisons of means 95% family-wise confidence level Fit: aov(formula = burst ~ tratamento + tempo + tratamento:tempo) $tempo diff lwr upr p adj 3-0 -0.006810191 -0.07288363 0.05926324 0.9671358 6-0 -0.140958351 -0.21023025 -0.07168645 0.0000175 6-3 -0.134148160 -0.20439788 -0.06389844 0.0000546 > TukeyHSD(model6,"tratamento:tempo") Tukey multiple comparisons of means 95% family-wise confidence level Fit: aov(formula = burst ~ tratamento + tempo + tratamento:tempo) $`tratamento:tempo` diff lwr upr p adj T2:0-T1:0 -0.021357143 -0.20493664 0.16222236 0.9999998 T3:0-T1:0 0.086392857 -0.10640418 0.27918989 0.9334458 T4:0-T1:0 -0.060232143 -0.25302918 0.13256489 0.9957673 T1:3-T1:0 -0.009357143 -0.20215418 0.18343989 1.0000000 T2:3-T1:0 -0.028428571 -0.22754850 0.17069135 0.9999979 T3:3-T1:0 -0.062232143 -0.25502918 0.13056489 0.9944118 T4:3-T1:0 0.069267857 -0.12352918 0.26206489 0.9865875 T1:6-T1:0 -0.138142857 -0.33726278 0.06097707 0.4599804 T2:6-T1:0 -0.091142857 -0.29026278 0.10797707 0.9236524 T3:6-T1:0 -0.162023810 -0.36927440 0.04522678 0.2795024 T4:6-T1:0 -0.183023810 -0.39027440 0.02422678 0.1359708 T3:0-T2:0 0.107750000 -0.06895140 0.28445140 0.6553978 T4:0-T2:0 -0.038875000 -0.21557640 0.13782640 0.9998325 T1:3-T2:0 0.012000000 -0.16470140 0.18870140 1.0000000 T2:3-T2:0 -0.007071429 -0.19065093 0.17650807 1.0000000 T3:3-T2:0 -0.040875000 -0.21757640 0.13582640 0.9997280 T4:3-T2:0 0.090625000 -0.08607640 0.26732640 0.8496727 T1:6-T2:0 -0.116785714 -0.30036521 0.06679378 0.5941646 T2:6-T2:0 -0.069785714 -0.25336521 0.11379378 0.9792305 T3:6-T2:0 -0.140666667 -0.33303479 0.05170146 0.3778522 T4:6-T2:0 -0.161666667 -0.35403479 0.03070146 0.1875552 T4:0-T3:0 -0.146625000 -0.33288463 0.03963463 0.2699473 T1:3-T3:0 -0.095750000 -0.28200963 0.09050963 0.8477318 T2:3-T3:0 -0.114821429 -0.30761847 0.07797561 0.6875834 T3:3-T3:0 -0.148625000 -0.33488463 0.03763463 0.2516698 T4:3-T3:0 -0.017125000 -0.20338463 0.16913463 1.0000000 T1:6-T3:0 -0.224535714 -0.41733275 -0.03173868 0.0095108 T2:6-T3:0 -0.177535714 -0.37033275 0.01526132 0.1002695 T3:6-T3:0 -0.248416667 -0.44960006 -0.04723327 0.0043573 T4:6-T3:0 -0.269416667 -0.47060006 -0.06823327 0.0012754 T1:3-T4:0 0.050875000 -0.13538463 0.23713463 0.9987106 T2:3-T4:0 0.031803571 -0.16099347 0.22460061 0.9999907 T3:3-T4:0 -0.002000000 -0.18825963 0.18425963 1.0000000 T4:3-T4:0 0.129500000 -0.05675963 0.31575963 0.4565814 T1:6-T4:0 -0.077910714 -0.27070775 0.11488632 0.9673298 T2:6-T4:0 -0.030910714 -0.22370775 0.16188632 0.9999930 T3:6-T4:0 -0.101791667 -0.30297506 0.09939173 0.8606043 T4:6-T4:0 -0.122791667 -0.32397506 0.07839173 0.6541132 T2:3-T1:3 -0.019071429 -0.21186847 0.17372561 1.0000000 T3:3-T1:3 -0.052875000 -0.23913463 0.13338463 0.9981734 T4:3-T1:3 0.078625000 -0.10763463 0.26488463 0.9555010 T1:6-T1:3 -0.128785714 -0.32158275 0.06401132 0.5194542 T2:6-T1:3 -0.081785714 -0.27458275 0.11101132 0.9539596 T3:6-T1:3 -0.152666667 -0.35385006 0.04851673 0.3220981 T4:6-T1:3 -0.173666667 -0.37485006 0.02751673 0.1584271 T3:3-T2:3 -0.033803571 -0.22660061 0.15899347 0.9999826 T4:3-T2:3 0.097696429 -0.09510061 0.29049347 0.8594125 T1:6-T2:3 -0.109714286 -0.30883421 0.08940564 0.7814023 T2:6-T2:3 -0.062714286 -0.26183421 0.13640564 0.9954629 T3:6-T2:3 -0.133595238 -0.34084583 0.07365535 0.5743224 T4:6-T2:3 -0.154595238 -0.36184583 0.05265535 0.3475958 T4:3-T3:3 0.131500000 -0.05475963 0.31775963 0.4324767 T1:6-T3:3 -0.075910714 -0.26870775 0.11688632 0.9729806 T2:6-T3:3 -0.028910714 -0.22170775 0.16388632 0.9999965 T3:6-T3:3 -0.099791667 -0.30097506 0.10139173 0.8755146 T4:6-T3:3 -0.120791667 -0.32197506 0.08039173 0.6767045 T1:6-T4:3 -0.207410714 -0.40020775 -0.01461368 0.0240599 T2:6-T4:3 -0.160410714 -0.35320775 0.03238632 0.1991828 T3:6-T4:3 -0.231291667 -0.43247506 -0.03010827 0.0111757 T4:6-T4:3 -0.252291667 -0.45347506 -0.05110827 0.0034935 T2:6-T1:6 0.047000000 -0.15211992 0.24611992 0.9996700 T3:6-T1:6 -0.023880952 -0.23113154 0.18336964 0.9999998 T4:6-T1:6 -0.044880952 -0.25213154 0.16236964 0.9998565 T3:6-T2:6 -0.070880952 -0.27813154 0.13636964 0.9909705 T4:6-T2:6 -0.091880952 -0.29913154 0.11536964 0.9379661 T4:6-T3:6 -0.021000000 -0.23607410 0.19407410 1.0000000 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

estrutura do arquivo
by Edmar Caldas 27 Mar '17

27 Mar '17

Como faço para visualizar o arquivo por inteiro. Quero visualizar as 31 linhas e 17 colunas. Estou usando Rstudio. Isso me ajuda em partes. head(append1, n=10)tail(append1, n=10)str(append1) Edmar

2 1

Ajuda na Função Append
by Edmar Caldas 27 Mar '17

27 Mar '17

Bom Dia! Estou precisando de uma ajuda na função append. Quando tenho as mesmas variáveis nos dois conjuntos de dados o código abaixo funciona. Mas quando tenho uma variável num conjunto e não tenho no outro, o mesmo código não funciona. 1) Situação que funciona. Arquivo A possui as seguintes variáveis. id nome sexo1 Edmar M2 Aline F Arquivo B possui as seguintes variáveis. id nome sexo.1 Adriana F2 Russo M código: append <-rbind(A,B) ok funciona. Mas quando tenho a seguinte situação. Me traz uma mensagem de erro. Arquivo A. Id nome sexo trabalho Arquivo B Id nome sexo lazer As variáveis trabalho e lazer cada uma está em arquivos diferentes e com isso não consigo unir pela função append simples. tem algum forma de corrigir isso? estou tentando usar esse código e não funciona. append <-rbind(lazer1,by.= "trabalho", lazer2, by = "lazer") Agradeço pela ajuda. Edmar

2 1

Fwd: Predição espacial de dado tipo "linha"
by Wagner Wolff 26 Mar '17

26 Mar '17

Olá pessoal da lista! Pretendo trabalhar com um tipo de dado que representa a atenuação de microondas entre duas antenas celulares. Sendo assim, estes dados representam uma "linha" (link) no espaço e não um ponto, como geralmente são feitas as análises geoestatśticas. Gostaria de saber o que muda na análise em comparação com a tradicional por pontos. Por exemplo, ao invés de usar a distância euclidiana entre as amostrar usaria uma distância ponderada pelo comprimento dessas "linhas"? Abaixo esta uma figura exemplificando o que mencionei acima. Att -- *Wagner Wolff, **PhD* "*Luiz de Queiroz**" College of Agriculture,* University of São Paulo Pádua Dias avenue11 | 13418-900| Piracicaba-SP| Brazil Phone: +55 19 982385582 <+55%2019%2098238-5582> http://orcid.org/0000-0003-3426-308X https://github.com/wwolff7 http://buscatextual.cnpq.br/buscatextual/visualizacv.do?id=K4463141A1

1 0

Parcela dividida - coletas de amostras dentro da parcela e no tempo
by Fabiola Santos 26 Mar '17

26 Mar '17

Boa tarde Tenho um experimento da seguinte forma: 3 tratamentos (tipo de manejo diferente) com 4 repeticoes em bloco. Parcela de 1ha As avaliações foram feitas em subparcelas onde escolhia 5 pontos aleatorios e coletava informações dentro dessas parcelas(de 1ha) em 10 plantas avaliando o numero de lagartas e especie. Porem foram feitas semanalmente as avaliações sendo em varios estagios do milho Gostarida de saber se o modelo a seguir estaria correto: O que eu nao estou conseguindo e fazer os teste de Tukey ou SK no proprio R. mod1 <- aov(PeqS~Rep+Trat*Est + Ponto + Error(Rep:Est, Rep:Ponto), data=x) Error: Rep:Est Df Sum Sq Mean Sq F value Pr(>F) Rep 3 8.010e-27 2.669e-27 0.746 0.541 Est 5 1.815e-26 3.630e-27 1.015 0.443 Residuals 15 5.365e-26 3.577e-27 Error: Within Df Sum Sq Mean Sq F value Pr(>F) Trat 2 0 0.0 0 1 Ponto 1 648 648.0 10689 <2e-16 *** Trat:Est 10 0 0.0 0 1 Residuals 3563 216 0.1 Obrigada Abraços Ótimo final de semana Vídeo Institucional https://www.youtube.com/watch?v=u1Z3_YZ2WpY <https://www.facebook.com/udipesquisaedesenvolvimento/> https://www.facebook.com/udipesquisaedesenvolvimento/ <https://www.linkedin.com/company/10401615?trk=tyah&trkInfo=clickedVertical% 3Acompany%2CclickedEntityId%3A10401615%2Cidx%3A2-1-2%2CtarId%3A1479812006594 %2Ctas%3AUDI%20PESQUISA%20E%20DESENV> https://www.linkedin.com/company/10401615?trk=tyah&trkInfo=clickedVertical%3 Acompany%2CclickedEntityId%3A10401615%2Cidx%3A2-1-2%2CtarId%3A1479812006594% 2Ctas%3AUDI%20PESQUISA%20E%20DESENV Atenciosamente / Best Regards / Mit freundlichen Grüßen / Saludos www.udipesquisa.om Esta mensagem, incluindo anexos, é confidencial, privilegiada e para o único benefício da pessoa devida ou endereçada. Caso seja recebido erroneamente, favor apagá-la do sistema, não faça cópias, não exponha seu conteúdo a uma terceira parte ou utilize-a para seu benefício pessoal ou de outra pessoa. Favor avisar ( <mailto:fabiolasantos@udipesquisa.com> fabiolasantos(a)udipesquisa.com) do seu recebimento, o mais rápido possível. Agradecemos a sua colaboração. The information contained in this e-mail and in the attached files are for the exclusive use of the addressee herein nominated, and may contain trade secrets, privileged and other confidential information, protected by the applicable laws. In case you are not the right addressee, you are hereby notified that any reviewing, reading, copying and/or distributing of this e-mail's content is strictly prohibited and unauthorized. Please, delete the e-mail's content and notify the sender ( <mailto:fabiolasantos@udipesquisa.com> fabiolasantos(a)udipesquisa.com) immediately. Thank you for your cooperation. --- Este email foi escaneado pelo Avast antivírus. https://www.avast.com/antivirus

1 0

[OFF-TOPIC]
by Julimar Pinto 25 Mar '17

25 Mar '17

Agradeço a oportunidade por esse tempo de learning-by-doing com a lista R-br.

1 0

Função aggregate.
by Edmar Caldas 24 Mar '17

24 Mar '17

Estou com um arquivo neste formato. Estou querendo agregar por paciente. Informando a média da coluna valor_int e quero também pegar a menor e a maior data da coluna data_int.Quero que o arquivo fique assim; Desenvolvi esse codigo. agora estou tentando colocar as datas. Edmar

3 2

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

grafico de distribuicao acumulada de uma variavel continua
by Nilson Guiotoku 18 Mar '17

18 Mar '17

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.

2 1

Box's M Test
by Pollyanna Alves 17 Mar '17

17 Mar '17

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

1 0

Teste de Normalidade Multivariada
by Pollyanna Alves 12 Mar '17

12 Mar '17

Olá pessoal, sou iniciante no programa R e estou tentando realizar um teste de normalidade multivariada. Dei uma pesquisada na internet e encontrei o comando > mshapiro.test(M), no pacote mvnormtest. Tenho uma planinha no formato .csv com dados numéricos apenas (castorianb). Digitei o comando > mshapiro.test (Castoria). Está correto? O comando gerou erros... Vocês teriam algum script deste teste para compartilhar? Erro gerado: Error in mshapiro.test(castorianb) : U[] is not a matrix with number of columns (sample size) between 3 and 5000 Grata 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

4 7

Re: [R-br] Digest R-br, volume 75, assunto 8
by sznelwar＠uol.com.br 12 Mar '17

12 Mar '17

1 0

Re: [R-br] Digest R-br, volume 75, assunto 8
by Maurício Sangiogo 12 Mar '17

12 Mar '17

Obrigado Walmes, Obrigado pela indicação. Mas acabei optando por realizar o cálculo dos quadrados médios combinados e graus de liberdade combinados manualmente mesmo, era isso que eu estava na dúvida... e de fato deu resultado diferente do método que estava de exemplo no pacote agricolae. Acredito eu que está certo como fiz. Se a alguém interessar, ou se disponibilizar de conferir, mando os dados e comandos: # Análise de um experimento em faixas com interação significativa: dados<- structure(list(local = c(1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 4, 1, 2, 3, 4), trat = c(0, 0, 0, 0, 30, 30, 30, 30, 60, 60, 60, 60, 0, 0, 0, 0, 30, 30, 30, 30, 60, 60, 60, 60, 0, 0, 0, 0, 30, 30, 30, 30, 60, 60, 60, 60), bloco = c(1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 2, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3, 3), yield = c(150, 163, 164, 171, 188, 174, 182, 182, 220, 170, 170, 182, 147, 141, 133, 140, 157, 166, 144, 143, 173, 165, 152, 170, 146, 140, 149, 150, 172, 158, 142, 160, 204, 163, 150, 175)), .Names = c("local", "trat", "bloco", "yield" ), row.names = c(NA, 36L), class = "data.frame") model<-with(dados,strip.plot(bloco, trat, local, yield)) # Se interação não significativa: comparison1<-with(dados,HSD.test(yield,trat,model$gl.a,model$Ea));comparison1 comparison2<-with(dados,HSD.test(yield,local,model$gl.b,model$Eb)); comparison2 # Se tiver interação significativa: QMea<- model$Ea ;Dfea<-model$gl.a QMeb<- model$Eb; Dfeb<-model$gl.b QMec<- model$Ec ; Dfec<-model$gl.c a <-3; b <-4 #Erro medio a e c QMeac <-(QMea+(b-1)*QMec)/b ##Graus de liberdade a e c Glac<-((QMea+(b-1)*QMec)^2) / ( ((QMea^2)/Dfea) + ((((b-1)*QMec)^2)/Dfec) ) #Erro medio b e c QMebc<-(QMeb+(a-1)*QMec)/a ##Graus de liberdade a e c Glbc<-((QMeb+(a-1)*QMec)^2) / ( ((QMeb^2)/Dfeb) + ((((a-1)*QMec)^2)/Dfec) ) #--------------------------------------------------------------------------------- # Comparação de médias por Tukey: ## Tratamentos dentro de locais: L1<-subset(dados,dados$local=="1") t1<-HSD.test(L1$yield,L1$trat,Glac,QMeac) t1 L2<-subset(dados,dados$local=="2") t2<-HSD.test(L2$yield,L2$trat,Glac,QMeac) t2 L3<-subset(dados,dados$local=="3") t3<-HSD.test(L3$yield,L3$trat,Glac,QMeac) t3 L4<-subset(dados,dados$local=="4") t4<-HSD.test(L4$yield,L4$trat,Glac,QMeac) t4 ## Locais dentro de tratamentos: A1<-subset(dados,dados$trat=="0") ta1<-HSD.test(A1$yield,A1$local,Glbc,QMebc) ta1 A2<-subset(dados,dados$trat=="30") ta2<-HSD.test(A2$yield,A2$local,Glbc,QMebc) ta2 A3<-subset(dados,dados$trat=="60") ta3<-HSD.test(A3$yield,A3$local,Glbc,QMebc) ta3 # FIM ? Abraços! .................. Me. 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: domingo, 12 de março de 2017 12:00 Para: r-br(a)listas.c3sl.ufpr.br Assunto: Digest R-br, volume 75, assunto 8 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: Desdobramento de um experimento em faixas (Walmes Zeviani) ---------------------------------------------------------------------- Message: 1 Date: Sat, 11 Mar 2017 17:44:03 -0300 From: Walmes Zeviani <walmeszeviani(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] Desdobramento de um experimento em faixas Message-ID: <CAFU=EkYuWZnFRExk8eNhcLSq3u1=OG0Yg=ZQUNeGoxQQc8obzQ(a)mail.gmail.com> Content-Type: text/plain; charset="utf-8" Verifique se essa postagem do histórico da lista serve: http://r-br.2285057.n4.nabble.com/R-br-Anova-em-faixas-td4655899.html. Caso R-br - [R-br] Anova em faixas<http://r-br.2285057.n4.nabble.com/R-br-Anova-em-faixas-td4655899.html> r-br.2285057.n4.nabble.com [R-br] Anova em faixas. Srs, Preciso de uma ajuda. Preciso fazer uma ANOVA de um experimento com três fatores, porém o problema está no delineamento que foi em faixas. Tenho o quadro da análise... não funcione, envie no corpo da mensagem os resultados que você tem como sendo os certos e o resultado que você está conseguindo. Não esqueça de disponibilizar acesso aos dados usados no corpo da mensagem (use dput()). À disposição. Walmes. -------------- Próxima Parte ---------- Um anexo em HTML foi limpo... URL: <http://listas.inf.ufpr.br/pipermail/r-br/attachments/20170311/02ec2ff1/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 8 ****************************************

1 0

Desdobramento de um experimento em faixas
by Maurício Sangiogo 11 Mar '17

11 Mar '17

Bom dia pessoal, Mais uma vez venho recorrer a vocês em relação a uma dúvida... Estou trabalhando com um experimento em DBC em faixas, porém tenho encontrado poucos exemplos em relação a comparação de médias no caso de interação. Encontrei um exemplo em uma aula na internet e um exemplo com os comandos no pacote agricolae (http://127.0.0.1:25234/library/agricolae/html/strip.plot.html) porém mesmo usando os mesmos valores a comparação de médias não está batendo entre o exemplo da aula e quando rodo usando os comandos do exemplo do agricolae. Alguém ai poderia conferir se esses comandos que envio logo a baixo estão corretos? Principalmente na parte referente a comparação de médias. (Seria um experimento em blocos completos no esquema em faixas.) Desde já muito obrigado! Abraços # Yield library(agricolae) data(huasahuasi) YIELD<-huasahuasi$YIELD market <- YIELD$y1da + YIELD$y2da non_market <- YIELD$y3da yield <- market + non_market model<-with(YIELD,strip.plot(block, clon, trt, yield)) comparison<-with(YIELD,LSD.test(yield,clon,model$gl.a,model$Ea)) comparison<-with(YIELD,LSD.test(yield,trt,model$gl.b,model$Eb)) # simple effects A<-model$data a<-nlevels(A$clon) b<-nlevels(A$trt) r<-nlevels(A$block) Ea<-model$Ea; Eb<-model$Ec; Ec<-model$Ec; gla<-model$gl.a; glb<-model$gl.b; glc<-model$gl.c; B <-tapply.stat(A[,4],A[,2:3],mean) std<-tapply.stat(A[,4],A[,2:3],function(x) sd(x)/sqrt(length(x))) B<-data.frame(B[,1:2],yield=B[,3],std=std[,3]) cmab<-(b-1)*Ec + Ea cmba<-(a-1)*Ec + Eb # order.group # Tukey ta <- qtukey(0.95,a,gla) tb <- qtukey(0.95,b,glb) tc <- qtukey(0.95,a*b,glc) tab<- ((b-1)*Ec*tc + Ea*ta)/cmab tba<- ((a-1)*Ec*tc + Eb*tb)/cmba # Comparison of clones by treatment groups<-by(B,B[,2], function(x) order.group(x$clon,x$yield,N=b*r, cmab,tab, std.err=x$std,parameter=0.5,console=FALSE)) groups # Comparison of treatments by clon groups<-by(B,B[,1], function(x) order.group(x$trt,x$yield,N=a*r,cmba,tba, std.err=x$std,parameter=0.5,console=FALSE)) groups .................. Me. Eng. Agr. Maurício Sangiogo TEL:(55) 96822030

2 1

Erro em fazendo pacote. .Random.seed e caracteres ASCII
by Pedro Emmanuel Alvarenga Americano do Brasil 10 Mar '17

10 Mar '17

Amigos de R, Submeti recentemente um pacote ao CRAN, e ha erros que preciso consertar antes da publicaçao. Mas dois deles eu nao estou sabendo como. https://win-builder.r-project.org/incoming_pretest/170309_202247_ems_0311/0… O primeiro deles é que o check acusa que no pacote há um banco chamado .Random.seed. Esse banco nao deveria estar la. Eu nao sei ao certo de onde veio, nem onde ele está. A unica coisa que sei é que ele nao está na pasta dos dados e nem está visivel no namespace do pacote. Imagino que retira-lo da arvore do pacote resolveria, mas nao estou sabendo como. O segundo problema é que o banco de dados de exemplo utilizado no pacote, possui níveis de variaveis categoricas que ficaram "desconfiguradas", por conta dos "ç" ou "ã". Isso faz com que caracteres fiquem irreconheciveis. Imagino que em algum momento entre importar editar e salvar, o encoding dos caracteres mudou, e retornar ao encoding original resolveria. Como faço pra voltar ao encoding original? Ou há alguma maneira mais inteligente de resolver e eliminar o erro de presença de caracteres nao ASCII? Alguem tem um bizu pra eu saber como arrumar essas questoes? Abraço forte, Pedro Brasil (:)=

1 0

Teste de Normalidade Multivariada
by Gustavo Agudelo 09 Mar '17

09 Mar '17

Olá Pollyanna, Tente usar as ferramentas do pacote MVN (Multivariate Normality Tests). Achei o material do pacote muito instrutivo e relativamente fácil de se seguir, tendo várias funções para testar normalidade multivariada e ferramentas gráficas úteis. Abraços, Gustavo. *====================================* *Gustavo A. Agudelo-Cantero, M.Sc.*Ph.D. Student in Physiology Laboratory of Ecophysiology and Evolutionary Physiology Institute of Biosciences University of São Paulo Rua do Matão 101, Travessa 14, Room 219. CEP 05508-090 (+55-11) 3091-7609 (+55-11) 94281-6827 São Paulo - SP Brazil ====================================

1 0

Silenciar os warnings
by Andre Oliveira 08 Mar '17

08 Mar '17

Boa tarde,sempre em um relatório cito os pacotes que usei nas referências, no entanto require(biotools) deixa a aviso solto no meio de meu relatório de um .rnw. Estou usando: ```{r,include=F} # Esconde avisos e comentários e códigos R require(knitr) opts_chunk$set(echo = F, comment="", message=F, warnings = F) suppressPackageStartupMessages(require(biotools,quietly = TRUE)) ```Alguém sabe como posso excluir tal aviso? André Oliveira Souza. Graduação em Matemática, mestrado em estatística aplicada.Instituto Federal de Educação, Ciência e Tecnologia do Espirito Santo. IFES

3 3

Re: [R-br] Silenciar os warnings
by Fernando Antonio de souza 08 Mar '17

08 Mar '17

Acho que isso deve funcionar Vc não falou como é a mensagem, mas caso seja de funções mascaradas por outra, isso deve resolver. require(pacote, warn.conflicts=FALSE) Att Em 08/03/2017 15:24, "Walmes Zeviani via R-br" <r-br(a)listas.c3sl.ufpr.br> escreveu: Eu uso suppressMessages(pacote). À disposição. Walmes. _______________________________________________ R-br mailing list R-br(a)listas.c3sl.ufpr.br https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br Leia o guia de postagem (http://www.leg.ufpr.br/r-br-guia) e forneça código mínimo reproduzível.

1 0