Junho 2014 - R-br - DINF UFPR

[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

RCurl: passando parâmetros e executando função javascript
by Éder Comunello 04 Dez '14

04 Dez '14

Colegas, boa tarde! Estou precisando de ajuda com relação ao pacote *RCurl*. Tomando por exemplo o endereço < http://hidroweb.ana.gov.br/Estacao.asp?Codigo=2254001>, ao acessá-lo via *browser *pode-se notar dois botões que executam trechos de código em *javascript *e geram links para o *download *de arquivos de dados: * Botão [Arquivo Access] remete ao comando "javascript:criarArq(2254001,2)" * Botão [Arquivo Texto] remete ao comando "javascript:criarArq(2254001,1)" Estou tentando usar o *RCurl *para passar diretamente os parâmetros para criação do *link *para o "Arquivo Texto". Isto feito, posteriormente pretendo usar o *link *gerado para *download* dos dados. Contudo, não estou tendo sucesso e ficarei muito grato por qualquer ajuda. Abaixo segue o código que tentei aplicar: ### <BEGIN> ### library(RCurl) myURL <- 'http://hidroweb.ana.gov.br/Estacao.asp?Codigo=2254001' ### URL resultante da operação via browser: ### myURL2 <- ' http://hidroweb.ana.gov.br/Estacao.asp?Codigo=2254001&CriaArq=true&TipoArq=1 ' myCurl <- getCurlHandle(cookiefile = "myCookies.txt", useragent="Mozilla/5.0", verbose = T, autoreferer = T, followlocation = T) myParams <- list(Codigo=2254001, CriaArq='true', TipoArq=1) ### Codigo: ID da estação ### CriaArq: true (executa função?), ### TipoArq: 1 (tipo texto) e 2 (Access) postForm(myURL, .params=myParams, curl=myCurl) getURL(myURL, curl=myCurl) ### Operação falha, retornando página original! ### Em caso de sucesso embutiria no código um link no formato: ### <a href="ARQ/A20131008-164720-14/CHUVAS.ZIP"> ### <END> ### -- Éder Comunello <c <comunello.eder(a)gmail.com>omunello.eder(a)gmail.com> Dourados, MS - [22 16.5'S, 54 49'W]

4 3

Efeito fixos e aleatórios
by Alisson Lucrecio 03 Jul '14

03 Jul '14

Caro Colegas da r-br e Walmes, Bom dia a todos. Eu posso incluir em um modelo uma variável independete no efeito aleatório e no efeito fixo assim como esta o modelo abaixo? Obrigado. > dput(estirpe_r_br)structure(list(Tratamento = structure(c(4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 7L, 7L, 7L, 7L, 3L, 3L, 3L, 3L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 7L, 7L, 7L, 7L, 3L, 3L, 3L, 3L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 7L, 7L, 7L, 7L, 3L, 3L, 3L, 3L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 7L, 7L, 7L, 7L, 3L, 3L, 3L, 3L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 7L, 7L, 7L, 7L, 3L, 3L, 3L, 3L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 7L, 7L, 7L, 7L, 3L, 3L, 3L, 3L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 7L, 7L, 7L, 7L, 3L, 3L, 3L, 3L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 7L, 7L, 7L, 7L, 3L, 3L, 3L, 3L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 7L, 7L, 7L, 7L, 3L, 3L, 3L, 3L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 4L, 4L, 4L, 4L, 5L, 5L, 5L, 5L, 6L, 6L, 6L, 6L, 7L, 7L, 7L, 7L, 3L, 3L, 3L, 3L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L), .Label = c("T_CN", "T_SN", "CIAT 899", "UFLA 02-100", "UFLA 02-127", "UFLA 02-68", "UFLA 04-195"), class = "factor"), Bloco = structure(c(1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L, 1L, 2L, 3L, 4L ), .Label = c("1", "2", "3", "4"), class = "factor"), Local = structure(c(3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 7L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 5L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 9L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L, 6L), .Label = c("Lavras", "Lavras 2", "Patos de Minas", "Patos de Minas 2", "PO", "PO 2", "Bambui", "GM", "Luminarias", "Pitangui"), class = "factor"), MSPA = c(22.96, 20.75, 28.56, 22.95, 24.91, 15.2, 12, 20.65, 18.39, 16.18, 14.45, 17.14, 13.02, 26.95, 8.45, 11.47, 23.62, 10.62, 10.35, 19.78, 34.52, 34.16, 35.99, 47.89, 15.32, 12.61, 10.45, 10.45, 24.93, 33.3, 52.25, 19.67, 53.04, 48.77, 34.64, 57.85, 52.43, 36.78, 37.15, 24.28, 23.48, 32.41, 42.63, 32.35, 40.17, 25.84, 31.77, 46.27, 53.27, 29.05, 54.03, 43.19, 28.18, 49.61, 47.09, 42.65, 72.1, 78.8, 58.2, 63.9, 62.5, 90.8, 61.6, 35, 61.4, 98, 64, 68.6, 73.3, 99.3, 72.3, 47.3, 74.4, 45.8, 79.7, 66.6, 140.2, 105.2, 111.8, 127.8, 79.1, 100.8, 91.6, 69.1, 67.8, 76.6, 82.7, 76.4, 81.2, 75.1, 82.3, 73.5, 81, 79.4, 112.1, 73.3, 56.6, 91.8, 86.2, 79.4, 65.2, 85.8, 75.5, 66, 136.2, 152.5, 126.3, 129.8, 68.4, 71.3, 95.5, 64.9, 66.64, 43.91, 72.9, 53.02, 46.5, 50.22, 65.15, 35.87, 49.14, 49.84, 92.12, 34.62, 81.94, 51.85, 72.23, 72.92, 67.5, 46.71, 48.24, 58.98, 66.78, 56.53, 59.89, 38.13, 54.03, 52.78, 51.5, 62.31, 20.06, 26.91, 11.52, 51.73, 24.31, 12.99, 17.35, 26.46, 8.59, 14.13, 25.43, 26.04, 15.53, 16.8, 19.27, 27.85, 12.36, 19.09, 11.05, 41.24, 12.03, 20.47, 17.5, 26.95, 13.49, 16.04, 16.72, 14.69, 42.73, 42.28, 87.11, 54.69, 61.6, 90.5, 38.43, 101.82, 48.37, 72.41, 53.5, 60.12, 41.68, 44.1, 46.95, 43.33, 47.92, 36.53, 33.5, 42.78, 62.4, 75.84, 51.5, 62.83, 44.6, 42.08, 44.29, 44.25, 29.81, 18.11, 21.24, 40.49, 22.49, 26.99, 33.88, 25.28, 22.88, 21.6, 26.87, 19.58, 21.38, 22.74, 18.34, 33.79, 29.69, 25.85, 19.9, 40.38, 18.96, 29.66, 27.93, 14.67, 20.34, 31.48, 22.15, 28.91, 80.38, 98.7, 96.36, 115.54, 56.14, 127.9, 76.04, 65.84, 40.4, 54.64, 36.2, 31.96, 30.76, 72.56, 51.02, 23, 69.26, 75, 51.82, 74.38, 57, 106.94, 80.16, 35.02, 37.36, 38.92, 49.76, 23.4, 29.56, 36.61, 56.75, 51.24, 42.03, 23.76, 46.05, 34.88, 46.99, 37.68, 30.95, 37.45, 37.92, 46.35, 46.75, 37.85, 27.32, 30.34, 56.91, 47.88, 29.41, 33.5, 33.94, 36, 40.16, 33.35, 35.55, 35.3), Irrigacao = structure(c(2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L), .Label = c("Sem", "Com"), class = "factor"), NT = structure(c(2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L ), .Label = c("2", "3", "4"), class = "factor"), Safra = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 3L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L), .Label = c("Seca", "Aguas", "Inverno"), class = "factor"), Prep_Solo = structure(c(1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 2L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L, 1L), .Label = c("CO", "PD" ), class = "factor")), .Names = c("Tratamento", "Bloco", "Local", "MSPA", "Irrigacao", "NT", "Safra", "Prep_Solo"), row.names = c(NA, -280L), class = "data.frame") > lmer1 <- lmer(MSPA ~ Tratamento + (1 + Tratamento | Local ), data = estirpe_r_br,+ control=lmerControl(optCtrl=list(maxfun=50000)))> > summary(lmer1)Linear mixed model fit by REML ['lmerMod'] Formula: MSPA ~ Tratamento + (1 + Tratamento | Local) Data: estirpe_r_br Control: lmerControl(optCtrl = list(maxfun = 50000)) REML criterion at convergence: 2253.7 Scaled residuals: Min 1Q Median 3Q Max -2.8159 -0.4992 -0.1138 0.4992 3.8454 Random effects: Groups Name Variance Std.Dev. Corr Local (Intercept) 1535.5 39.19 TratamentoT_SN 390.8 19.77 -0.90 TratamentoCIAT 899 573.7 23.95 -0.92 0.89 TratamentoUFLA 02-100 716.9 26.77 -0.80 0.71 0.95 TratamentoUFLA 02-127 635.4 25.21 -0.83 0.75 0.92 0.95 TratamentoUFLA 02-68 370.3 19.24 -0.91 0.98 0.92 0.77 0.84 TratamentoUFLA 04-195 496.5 22.28 -0.88 0.98 0.94 0.79 0.79 0.97 Residual 150.7 12.28 Number of obs: 280, groups: Local, 10 Fixed effects: Estimate Std. Error t value (Intercept) 60.399 12.543 4.815 TratamentoT_SN -17.636 6.828 -2.583 TratamentoCIAT 899 -15.847 8.056 -1.967 TratamentoUFLA 02-100 -9.145 8.901 -1.027 TratamentoUFLA 02-127 -10.860 8.430 -1.288 TratamentoUFLA 02-68 -15.496 6.676 -2.321 TratamentoUFLA 04-195 -16.052 7.562 -2.123 Correlation of Fixed Effects: (Intr) TrT_SN TCIAT8 TUFLA02-10 TUFLA02-12 TUFLA02-6 TratmntT_SN -0.857 TrtmCIAT899 -0.888 0.837 TUFLA02-100 -0.784 0.679 0.900 TUFLA02-127 -0.812 0.716 0.878 0.909 TrUFLA02-68 -0.868 0.903 0.854 0.729 0.792 TUFLA04-195 -0.851 0.913 0.884 0.755 0.759 0.898 -- Alisson Lucrecio da Costa

3 3

Função 'getManyCsv"
by Leonardo Monteiro 01 Jul '14

01 Jul '14

Boa tarde, Estou tentando usar a função "getManyCsv", a qual estaria dentro do pacote "ChemoSpec", porém estou recebendo a mensagem: "Error, could not find function "getManyCsv" " Alguem saberia alguma solucao? obrigado, Leonardo.

2 1

Problema em abrir conexão com source() do Dropbox/Public
by walmes . 01 Jul '14

01 Jul '14

Saudações, Muitas das funções que disponibilizo em cursos eu hospedo na pasta pública do Dropbox (Dropbox/Public). O procedimento para mim e também para os alunos sempre funcionou usando o link que o Dropbox dá ao clicar com o botão direito do mouse "Dropbox > Copy Public Link". O link vem com https e o "s" tem que ser removido. Isso sempre funcionou mas desse a última sexta-feira pra cá não funcionou mais. O erro que dá é mostrado no código abaixo. Googlei e os erros apontam problemas de proxy (praticamente em windows). Não creio que seja meu caso. Tentei algumas soluções lidas mas não deu. Alguém experimentou algo semelhante e/ou tem alguma indicação? > url <- "http://dl.dropboxusercontent.com/u/48140237/aac.R" > browseURL(url) > download.file(url, basename(url)) trying URL 'http://dl.dropboxusercontent.com/u/48140237/aac.R' Error in download.file(url, basename(url)) : *cannot open URL* 'http://dl.dropboxusercontent.com/u/48140237/aac.R' > source(url) Error in file(filename, "r", encoding = encoding) : *cannot open the connection* > sessionInfo() R version 3.1.0 (2014-04-10) Platform: x86_64-pc-linux-gnu (64-bit) locale: [1] LC_CTYPE=en_US.UTF-8 LC_NUMERIC=C [3] LC_TIME=pt_BR.UTF-8 LC_COLLATE=en_US.UTF-8 [5] LC_MONETARY=pt_BR.UTF-8 LC_MESSAGES=en_US.UTF-8 [7] LC_PAPER=pt_BR.UTF-8 LC_NAME=C [9] LC_ADDRESS=C LC_TELEPHONE=C [11] LC_MEASUREMENT=pt_BR.UTF-8 LC_IDENTIFICATION=C attached base packages: [1] stats graphics grDevices utils datasets methods base loaded via a namespace (and not attached): [1] grid_3.1.0 lattice_0.20-29 nlme_3.1-117 > Grato. Walmes. ========================================================================== Walmes Marques Zeviani LEG (Laboratório de Estatística e Geoinformação, 25.450418 S, 49.231759 W) Departamento de Estatística - Universidade Federal do Paraná fone: (+55) 41 3361 3573 skype: walmeszeviani homepage: http://www.leg.ufpr.br/~walmes linux user number: 531218 ==========================================================================

5 6

combinar Cluster e Ordination (PCA)
by Jackson Rodrigues 01 Jul '14

01 Jul '14

Olá todos, Gostaria de uma ajuda para plotar Principal Components Analysis (PCA) e clusters juntos. Possuo dados ambientais de 25 locais espalhados em 5 diferentes ecossistemas que quando agrupados em 5 clusters, locais de ecossitemas diferentes se juntam no mesmo grupo. Logo, uma bagunça. Assim, quero plotar os 2 juntos (PCA e Clusters) de maneira que, locais pertencentes ao mesmo ecossistema possuam as mesmas CORES, e locais agrupados no mesmo cluster, possuam mesma FORMA. Por exemplo: 1 local da Floresta de Araucária e 1 local da Floresta Decidual (logo com cores diferentes azul e verde respectivamente) estão no mesmo cluster (logo, representados pela mesma figura, triângulo). Como fazer essa combinação? Abaixo segue os comandos que estou usando para fazer minhas combinação, no entanto. Nesse código, consigo juntar PCA e Cluster. no entanto sem diferenciar as cores de cada local de acordo com eu ecossistema. #PCA mydata.pca<-rda(mydata) #Cluster mydata.w <- hclust(dist(Sqchord.mydata), "ward") plot(mydata.w, hang=-1) rect.hclust(mydata.w, 5) # Cut the dendrogram to yield 5 groups gr <- cutree(mydata.w, k=5) grl <- levels(factor(gr)) # Get the site scores, scaling 1 sit.sc1 <- scores(mydata.pca, display="w", scaling=3) # Plot the sites with cluster symbols and colours (scaling 3) p <- plot(mydata.pca, display="wa", scaling=3, type="n", main="Mydata PCA and clusters") abline(v=0, lty="dotted") abline(h=0, lty="dotted") for (i in 1:length(grl)) { points(sit.sc1[gr==i,], pch=(14+i), cex=2.5, col=i+1)} text(sit.sc1, row.names(Pre_euro_veg.1.all), cex=0.7, pos=3) # Add legend interactively legend(locator(1), paste("Cluster", c(1:length(grl))), pch=14+c(1:length(grl)), col=1+c(1:length(grl)), pt.cex=2) Desde já, agradeço por qualquer ajuda. Jackson M. Rodrigues

2 1

Re: [R-br] Questão de concurso verossimilhança...
by Leal 01 Jul '14

01 Jul '14

Está errado o gabarito. É uma distribuição uniforme. Lambada emv = máximo, neste caso 24. Conectado pela Motorola andrebvs(a)bol.com.br escreveu: >_______________________________________________ >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.

4 3

Re: [R-br] Anosim_dúvidas (Silma Leite Rocha)
by Marcos Vital 30 Jun '14

30 Jun '14

Olá, Silma A anosim é um teste de aleatorização de matrizes (mais ou menos parecido com o teste de Mantel, que é bem mais conhecido) que "pergunta" se em uma matriz de distância (ou de similaridade, tanto faz), as distâncias dentro de uma classe é diferente das distâncias entre diferentes classes. Em que contexto você precisa usar a análise? Ela é muito usada em ecologia de comunidades, quando temos uma matriz de distância entre unidades amostrais com dados de presença ou abundância de espécies. A matriz é normalmente construída com uma métrica apropriada para este tipo de dados (Jaccard, Bray-Curtios, etc), e as unidades amostrais pertencem a diferentes classes (estações, tipos de ambiente, etc). A análise pode ser feita pelo comando anosim() do pacote vegan. Você irá precisar de dois objetos: a matriz, que poderá ser criada com um comando como dist() ou vegdist() (este último contem boa parte das métricas típicas da ecologia de comunidades); e um vetor que contenha todas as classes, para cada unidade amostral. Veja o exemplo a seguir, com dados inventados, para ter uma ideia de como funciona: ####Início do Script#### library(vegan) dados<-structure(list(UA = 1:8, Ambiente = structure(c(1L, 1L, 1L, 1L, 2L, 2L, 2L, 2L), .Label = c("A", "B"), class = "factor"), sp1 = c(2L, 0L, 3L, 1L, 3L, 1L, 3L, 0L), sp2 = c(1L, 4L, 1L, 4L, 0L, 3L, 4L, 5L), sp3 = c(1L, 0L, 5L, 0L, 4L, 4L, 3L, 1L), sp4 = c(1L, 3L, 2L, 3L, 0L, 0L, 3L, 5L), sp5 = c(5L, 2L, 3L, 4L, 2L, 5L, 4L, 4L)), .Names = c("UA", "Ambiente", "sp1", "sp2", "sp3", "sp4", "sp5"), class = "data.frame", row.names = c(NA, -8L)) dados #Apenas para visualizar matriz<-vegdist(dados[,3:7], method="bray") classes<-dados$Ambiente anosim(matriz, classes) ####Fim do Script#### Recomendo que você dê uma olhada no material do pacote vegan que está no CRAN, pois ele apresenta um manual bem detalhado. Espero ter ajudado Atenciosamente Marcos Em 30 de junho de 2014 11:17, <r-br-request(a)listas.c3sl.ufpr.br> escreveu: > Enviar submissões para a lista de discussão R-br para > r-br(a)listas.c3sl.ufpr.br > > Para se cadastrar ou descadastrar via WWW, visite o endereço > https://listas.inf.ufpr.br/cgi-bin/mailman/listinfo/r-br > ou, via email, envie uma mensagem com a palavra 'help' no assunto ou > corpo da mensagem para > r-br-request(a)listas.c3sl.ufpr.br > > Você poderá entrar em contato com a pessoa que gerencia a lista pelo > endereço > r-br-owner(a)listas.c3sl.ufpr.br > > Quando responder, por favor edite sua linha Assunto assim ela será > mais específica que "Re: Contents of R-br digest..." > > > Tópicos de Hoje: > > 1. Anosim_dúvidas (Silma Leite Rocha) > 2. Re: Questão de concurso verossimilhança... (Leal) > 3. Re: Questão de concurso verossimilhança... (Leal) > 4. Re: Questão de concurso verossimilhança... (Paulo Justiniano) > 5. Função 'getManyCsv" (Leonardo Monteiro) > 6. Resíduos agrupados dataset: Iris (Helder Gramacho) > 7. Re: Resíduos agrupados dataset: Iris > (Leonardo Ferreira Fontenelle) > 8. combinar Cluster e Ordination (PCA) (Jackson Rodrigues) > 9. Re: Resíduos agrupados dataset: Iris (Vinicius Brito Rocha) > 10. Re: Resíduos agrupados dataset: Iris (Helder Gramacho) > > > ---------------------------------------------------------------------- > > Message: 1 > Date: Sun, 29 Jun 2014 15:39:05 -0300 > From: Silma Leite Rocha <silmalrocha22(a)gmail.com> > To: r-br(a)listas.c3sl.ufpr.br > Subject: [R-br] Anosim_dúvidas > Message-ID: > <CAKsbjo59LB33j= > TsL2cdtU5enRReUETTMonf2Sz+ZpHK+e41-g(a)mail.gmail.com> > Content-Type: text/plain; charset="utf-8" > > Boa tarde, > > preciso aprender a fazer uma análise - Anosim. Dei uma olhadinha no grupo e > vi que já fizeram uma postagem sobre isso mas não entendi muito. Por > exemplo: uma pessoa cita que são necessárias duas tabelas, o que seria essa > outra tabela? > Alguém teria o script para me ajudar? > Nunca trabalhei com este tipo de análise. Fiz os meus testes utilizando o > G-test em um artigo, e o referee solicitou que eu substituísse a análise > pela Anosim. > > Desde já agradeço o auxílio. > -------------- Próxima Parte ---------- > Um anexo em HTML foi limpo... > URL: < > http://listas.inf.ufpr.br/pipermail/r-br/attachments/20140629/001bc91c/atta… > > > -- Marcos Vinícius Carneiro Vital Universidade Federal de Alagoas Instituto de Ciências Biológicas e da Saúde Setor de Biodiversidade e Ecologia

2 1

Resíduos agrupados dataset: Iris
by Helder Gramacho 30 Jun '14

30 Jun '14

Bom dia pessoal, Estava realizando um tutorial que utilizava o dataset: Iris. Após a análise de regressão o gráfico dos resíduos x valores ajustados mostrava que os pontos estavam dois grupos separados, como na figura abaixo: [image: Imagem inline 2] Gostaria de saber quando este gráfico apresenta este padrão, quais seriam as causas e quais as possíveis soluções. Segue o cmr: ### <code r> data(iris) attach(iris) iris.lm<-lm(Petal.Width ~ Petal.Length, data=iris) iris.lm summary(iris.lm) confint(iris.lm) par(mfrow=c(2,2),xpd=F) plot(iris.lm) ### </code> Desde já agradeço as contribuições, *Hélder Gramacho * Recife-PE / *agrohelder(a)gmail.com <agrohelder(a)hotmail.com>*

3 3