Outubro 2018 - R-br

[OFF -Tópic] Cotas em concursos para negros e índios!
by Andre Oliveira 12 Jun '49

12 Jun '49

Pessoal, tem um colegar meu que vou fazer um concurso ele é descendente de negro ele é pardo. Ele veio com a pergunta como provar que sou negro, pois, tem vagas reservada para negros e índios? Eu sinceramente não sei e não tenho ideia pois nunca tive esta curiosidade. O que de fato comprova sua raça? Algum documento? O que são usados em concursos para provar? obrigado André Oliveira Souza

5 6

Legenda eixo X
by Luciano F. Sgarbi 10 Jun '49

10 Jun '49

Prezados, Estou tentando centralizar a legenda do eixo X em uma figura com vários gráficos. Estou conseguindo centralizar a legenda da linha inferior como no exemplo, mas as linhas do meio não. par(mfrow=c(2,2)) plot(rnorm(100),xlab="",ylab="y1") plot(rnorm(100),xlab="",ylab="y2") mtext(side=1,cex=2,line=-1.8,text="Variável Explanatória 1",outer=TRUE) ### o problema está aqui, pois a legenda desce para a parte inferior da figura. plot(rnorm(100),xlab="",ylab="y1") plot(rnorm(100),xlab="",ylab="y2") mtext(side=1,cex=2,line=-1.8,text="Variável Explanatória 2",outer=TRUE)## Esse dá certo Alguém podia dar uma luz de como eu posso fazer isso? Abraços -- *Luciano F. Sgarbi* Mestrando em Ecologia e Evolução Laboratório de Ecologia de Insetos (sl. 222) Departamento de Ecologia Instituto de Ciências Biológicas - ICB 1 Universidade Federal de Goiás, campus II Goiânia-GO 74001-970 Brazil

2 2

Re: [R-br] Diferenças entre sistemas operacionais (geoR)
by Yury Duarte 07 Nov '18

07 Nov '18

Boa noite Paulo, como vai? Agradeço desde já pela atenção. Tentei remover ao máximo as partes que não afetam na função e anexei o script no arquivo texto que segue. Att Em qui, 4 de out de 2018 19:23, Paulo Justiniano Ribeiro Junior < paulojus(a)ufpr.br> escreveu: > Olá! > > krige.conv() só utiliza simulações caso seja indicado em output.options > o número de simulações da distribuição preditiva. > Neste caso são geradas simulações. > As simulações só são geradas se: > 1. especificado explicitamente > 2. se for utilizada alguma transformação com lambda diferente de 0, 1/2 > > O seu exemplo se enquadra em algum caso? > > A msg tb indica que voce colocou algum em micro.scale e seria bom saber o > que deseja com isto. > > Se dados e códigos são grande Envie o comando completo da chamada da > krige.conv() para vermos as opções > > P.J. > > > > > ------------------------------ > *From: *"Yury Duarte" <yurynepomuceno(a)gmail.com> > *To: *"Thiago V. dos Santos" <thi_veloso(a)yahoo.com.br> > *Cc: *"a lista Brasileira oficial de discussão do programa R." < > r-br(a)listas.c3sl.ufpr.br>, paulojus(a)ufpr.br, "p diggle" < > p.diggle(a)lancaster.ac.uk> > *Sent: *Thursday, 4 October, 2018 6:08:32 PM > *Subject: *Re: [R-br] Diferenças entre sistemas operacionais (geoR) > > Não sei se existe um gerador de números aleatórios dentro da função, mas > estou passando pra ela todos os parâmetros já pré definidos antes dela ser > chamada. > Vou colocar um set.sead() no início do código, conforme sugestão do César > para então comparar os valores gerados nos diferentes sistemas operacionais > e ver o que acontece. > > Em qui, 4 de out de 2018 17:45, Thiago V. dos Santos < > thi_veloso(a)yahoo.com.br> escreveu: > >> Da maneira que eu vejo, a unica explicacao plausivel seria se a >> krige.conv usasse um algoritmo gerador de numeros aleatorios. Caso voce >> tenha experiencia com esmagamento de bugs, vc pode conferir isso no codigo >> dela. Ou entao esperar o feedback dos autores. >> >> Greetings, >> -- Thiago V. dos Santos >> >> Postdoctoral Research Fellow >> Department of Climate and Space Science and Engineering >> University of Michigan >> >> >> On Thursday, October 4, 2018, 4:39:08 PM EDT, Yury Duarte < >> yurynepomuceno(a)gmail.com> wrote: >> >> >> Sim, Thiago. >> Até o último passo os resultados são iguais. A mesma grade com o mesmo >> número de pontos é gerada. As mesmas informações são passadas como >> parâmetro para a função, mas os resultados não batem. >> Nem mesmo em casos onde a função não quebra, os valores gerados no Ubuntu >> não batem com os valores gerados em Windows. Isso está me intrigando muito. >> >> Em qui, 4 de out de 2018 17:23, Thiago V. dos Santos < >> thi_veloso(a)yahoo.com.br> escreveu: >> >> Nao deveria ser assim. >> >> O seu codigo produz resultados identicos ate o momento em que voce roda a >> krige.conv? >> >> Em caso positivo, acho que seria uma boa ideia reportar esse erro ao >> Paulo (paulojus(a)ufpr.br) ou ao Peter (p.diggle(a)lancaster.ac.uk) >> >> Greetings, >> -- Thiago V. dos Santos >> >> Postdoctoral Research Fellow >> Department of Climate and Space Science and Engineering >> University of Michigan >> >> >> On Thursday, October 4, 2018, 3:42:35 PM EDT, Yury Duarte via R-br < >> r-br(a)listas.c3sl.ufpr.br> wrote: >> >> >> Boa tarde a todos os coleas! >> >> Fiz um código simples para interpolar valores que tenho em uma grade que >> construo dentro do próprio código. Para isso, estou utilizando um modelo >> esférico de interpolação por krigagem através da função krige.conv(). >> Estou achando muito estranho o fato de que, quando rodo o script em >> ambiente Windows, não tenho problemas e meu script gera as saídas >> exatamente como o esperado, entretanto, quando rodo o mesmo script (com o >> mesmo banco de dados) em ambiente Linux (Ubuntu), a função krige.conv() >> corrompe e o erro gerado é: >> krige.control: micro.scale must be in the interval [0, nugget] >> >> Não encontrei nada em outros fóruns que pudesse me ajudar para explicar >> isso. >> Caso alguém tenha interesse em testar exatamente o mesmo conjunto de >> dados que estou utilizando, informem que passo um link no google drive, já >> que o código e os dados de entrada são grandes demais para serem anexados >> no corpo do e-mail/nos anexos. >> >> Desde já, agradeço a atenção e a ajuda de todos! >> >> Yury Duarte >> Engenheiro Agrônomo - ESALQ/USP >> _______________________________________________ >> 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. >> >> >

2 2

Diferença resultados agricolae::HSD.test e lsmeas
by Fernando Antonio de souza 05 Nov '18

05 Nov '18

Pessoal, poque o resultado dos testes HSD.test e lsmeans diferem neste exemplo? O resultado fornecido pelo lsmeas está em acordo com o material de onde tirei esse exemplo e que havia sido rodado no SAS. Já o HSD.test difere tanto no valor das médias estimadas como na significancia. install.packages(c("agricolae",lsmeans")) library(agricolae) library(lsmeans) COVARIANCIA<-read.csv(" https://www.dropbox.com/s/0uvwqnnvsdsvphc/Covariancia.csv?raw=1",head=TRUE) modeloCOVAR<-lm(GANHO~PESOINIC+DIETA,data=COVARIANCIA) anova(modeloCOVAR) HSD.test(modeloCOVAR,"DIETA",console=TRUE) lsmeans(modeloCOVAR,pairwise~DIETA,adjust="tukey") -- ========================================= Fernando Souza Zootecnista, DSc. Produção e Alimentação Animal Celular: (31)99796-8781 (Vivo) E-mail:nandodesouza@gmail.com <e-mail%3Anandodesouza(a)gmail.com> Lattes: http://lattes.cnpq.br/6519538815038307 Blog: https://producaoanimalcomr.wordpress.com/ ==========================================

3 2

Análise de experimentos em faixas (medidas repetidas no tempo)
by Andre Oliveira 31 Out '18

31 Out '18

bom dia,de todos os métodos análises que busquei em livros, no google e aqui na lista r-br, sempre achei análises de experimento em blocos, achei nas bibliotecas ExDes e na agricalae tal análise, mas, preciso de script para o DIC. Alguém tem alguma dica de material para eu dar uma estudada? Sei que posso usar a lme4, mas, não estou seguro nas testes a posteriori. Grato

2 2

Teste qui-quadrado
by Lucas Stempkowski 30 Out '18

30 Out '18

Boa tarde, Preciso realizar o teste qui-quadrado para comparação de frequências de classes de severidade de doenças, porém não tenho ideia de como fazê-lo utilizando o R. Eu tenho cinco classes de severidade (1 a 5), 4 doses de N e duas cultivares. A intenção é comparar as frequências entre doses e cultivares. Se alguém puder me ajudar agradeço. -- ----------------------------------------------------- Lucas Antonio Stempkowski Eng. Agronômo - UFFS Mestrando em Produção Vegetal CAV/UDESC +55 54 99153-0124/ 98448-4633 --------------------------------------------------

2 1

[off-topic] Divulgação VIII Encontro dos Alunos em Estatística - ESALQ/USP
by Eduardo Junior 25 Out '18

25 Out '18

Prezados(as) redistas, O Programa de Pós-graduação em Estatística e Experimentação Agronômica da ESALQ está organizando o *VIII Encontro dos Alunos em Estatística e Experimentação*. O evento ocorrerá nos dias *21, 22 e 23 de novembro* e contará com minicurso, conferências, mesa redonda e comunicações orais. Nessa edição, o minicurso será "*Modelos de regressão não linear: teoria e aplicações*", ministrado pelo Prof. Dr. *Walmes Zeviani* (LEG/UFPR). As conferências serão proferidas por excelentes pesquisadores da área: Profa. Dra. *Clarice Demétrio* (ESALQ/USP), Prof. Dr. *Afrânio Vieira* (UFSCar), Profa. Dra. *Vera Tomazella* (UFSCar) e Profa. Dra. *Luzia Trinca* (UNESP/Botucatu). Confira os detalhes da programação e demais informações para submissão de trabalhos inscrições no endereço *https://esalq-ppgeea.github.io/encontro2018*/ <https://esalq-ppgeea.github.io/encontro2018/>. Save the date! [image: cartaz.png] -- Eduardo Elias Ribeiro Junior <https://jreduardo.github.io/> Master's candidate in Statistics and Agricultural Experimentation (ESALQ-USP) Statistics and Geoinformation Laboratory (LEG) URL: http://leg.ufpr.br/~eduardojr

1 0

Re: [R-br] COMPARANDO CURVAS [DÚVIDA]
by sznelwar＠uol.com.br 23 Out '18

23 Out '18

1 0

COMPARANDO CURVAS [DÚVIDA]
by Diego Vieira 22 Out '18

22 Out '18

Boa tarde pessoal, Estou querendo fazer a comparação de duas curvas ajustadas utilizando a função "nlsLM" do pacote "minpack.lm". No entanto, não consigo ajustar as duas curvas em um único nlsLM, para depois usar a anova para realizar um teste da razão de verossimilhança para testar a restrição em A, B e C. Alguém pode me ajudar a solucionar esse problema. Desde já agradeço muito a todas as sugestões. Além disso, preciso usar esse pacote por causa do algoritmo de Levenberg-Marquardt. Segue o script reproduzível. if(!require("minpack.lm")){install.packages("minpack.lm",dependencies=T);require("minpack.lm")} if(!require("lattice")){install.packages("lattice",dependencies=T);require("lattice")} d <- expand.grid(parA=c(1,1.99,3,4,5,6,6), sistema = factor(c("TRAT1","TRAT2"))) d$diversidade <- c(55.000, 19.572, 10.231, 7.225, 6.004, 5.4022, 5.0566, 62.000, 31.6952, 21.4178, 17.6125, 15.7654, 14.7072, 14.0334) #ajusta as duas curvas em um unico nlsLM n0 <- nlsLM(diversidade ~ A[sistema]/(1+B[sistema]*exp(-C[sistema]*parA)), data = d, start = list(A = c(-2.790,11.106), B = c(-1.050,-0.821), C = c(-0.090, 0.248))) #Mensagem de erro: # Error in dimnames(x) <- dn : # length of 'dimnames' [2] not equal to array extent # In addition: Warning message: # In matrix(out$hessian, nrow = length(unlist(par))) : # data length [9] is not a sub-multiple or multiple of the number of rows [6] #ajusta o modelo em que A, B e C seja comum aos dois tratamentos n1 <- nlsLM(diversidade ~ A/(1+B*exp(-C*parA)), data = d, start = list(A = 4.158, B = -0.9355, C = 0.079)) anova(n0,n1) -- Diego dos Santos Vieira Engenheiro Florestal - UFRA Mestre em Ciência Florestal - UFVJM Doutorando em Ciência Florestal - UFVJM

2 2

Anova fatorial com medidas repetidas (e aninhamento?)
by Thallita Grande 17 Out '18

17 Out '18

Olá pessoal. Desejo comparar a variável "PC" (variável dependente contínua) entre 3 fatores (variáveis independentes categóricas: rede, ambiente e sensibilidade). Como tenho mais de uma medida para a mesma paisagem, creio que devo fazer uma ANOVA fatorial com medidas repetidas. No entanto, creio que haja um aninhamento entre os níveis da variável "sensibilidade", pois o tipo "sensível" é sempre um subconjunto do tipo "tolerante" (está contido em). Como considero esse aninhamento no teste? Por enquanto, meu código é: resultado=aov(PC ~ rede + ambiente + sensibilidade + Error(factor(paisagem))) paisagem PC rede ambiente sensibilidade 35 238.00 Sim florestal sensivel 35 247.00 Sim florestal tolerante 35 6.76 Sim savanico tolerante 35 228.00 Não florestal sensivel 35 228.00 Não florestal tolerante 35 5.15 Não savanico tolerante 36 113.00 Sim florestal sensivel 36 123.00 Sim florestal tolerante 36 144.00 Sim savanico sensivel 36 160.00 Sim savanico tolerante 36 87.50 Não florestal sensivel 36 87.50 Não florestal tolerante 36 127.00 Não savanico sensivel 36 127.00 Não savanico tolerante Obrigada. ----------------------------- Thallita Grande, M.Sc. Doutorado em Ecologia Instituto de Ciencias Biologicas Universidade de Brasilia - Campus Darcy Ribeiro Skype: thallita.grande CRBio 080256/04-D

2 4