Dezembro 2019 - 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

componente da variância no modelo de efeitos aleatórios
by Luiz Leal 12 Jan '20

12 Jan '20

Prezados, boa tarde. No modelo de efeitos aleatórios estimamos a componente da variância conforme abaixo. y <-c(2370, 1687, 2592, 2283, 2910, 3020, 1282, 1527, 871, 1025, 825, 920, 562, 321, 636, 317, 485, 842, 173, 127, 132, 150, 129, 227, 193, 71, 82, 62, 96, 44)x <- as.factor(c("A", "A", "A", "A", "A", "A", "B", "B", "B", "B", "B", "B", "C", "C", "C", "C", "C", "C", "D", "D", "D", "D", "D", "D", "E", "E", "E", "E", "E", "E"))summary(aov(y~x)) Mas não sei como calcular quando utilizo mínimos quadrados generalizados. y <-c(2370, 1687, 2592, 2283, 2910, 3020, 1282, 1527, 871, 1025, 825, 920, 562, 321, 636, 317, 485, 842, 173, 127, 132, 150, 129, 227, 193, 71, 82, 62, 96, 44) x <- as.factor(c("A", "A", "A", "A", "A", "A", "B", "B", "B", "B", "B", "B", "C", "C", "C", "C", "C", "C", "D", "D", "D", "D", "D", "D", "E", "E", "E", "E", "E", "E"))library(nlme) modelo <- gls(y~x, weights = varExp() , method = "ML") Alguém sabe como calculo as componentes da variância neste caso? Desde já agradeço Luiz

2 4

Tinn-R project: new version (6.01.01.03) released
by Jose Claudio Faria 23 Dez '19

23 Dez '19

Hi, A new version of Tinn-R project (6.01.01.03) was released today. Download: - https://nbcgib.uesc.br/tinnr/en/download - https://sourceforge.net/projects/tinn-r/files/Tinn-R setup/ What is new: - https://nbcgib.uesc.br/tinnr/en/download#patch - https://sourceforge.net/projects/tinn-r/files/Tinn-R setup/6.1.1.3/ Best, ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ Jose Claudio Faria UESC/DCET/Brasil joseclaudio.faria at gmail.com Telefones: 55(73)3680.5545 - UESC 55(73)99966.9100 - VIVO ///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\///\\\ If you have software to deal with statistics, you have arms; if you have good software, you have arms and legs; if you have software like R, you have arms, legs and wings... the height of your flight depends only on you!

1 0

Erro predição - Função predict
by Andre Oliveira 17 Dez '19

17 Dez '19

Boa noite a todos! Estou procurando obter a métrica RMSE e R2 entre os objetos pred e diameter no modelo abaixo, mas está ocorrendo um erro. Parece que reclama de missing, mas ainda não consegui resolver, mesmo depois de muitas pesquisas. O interesse é pegar o predito levando em consideração os dois efeitos do modelo. Alguém teria alguma dica e sugestão? # CMDR rm(list=ls(all=TRUE))require(lme4)data(Penicillin, package="lme4")attach(Penicillin) library(MCMCglmm)model <- MCMCglmm(diameter ~ 1, random = ~ plate + sample, verbose= T, nitt=200000,burn=10000,thin=10, rcov = ~ units, data=Penicillin) Penicillin$pred=predict(model, data=Penicillin, marginal = NULL) Error in object$Sol[, missing.random, FALSE] <- rnorm(nrow(object$Sol) * : incorrect number of subscripts Pedir os valores preditos como abaixo resolve, mas, demora muito tempo em caso de bases de dados grandes. Penicillin$pred=predict(model, data=Penicillin, marginal = NULL, interval="prediction") att,.André

1 0

Manter a estrutura do objeto da funcao randomForest
by Nilson Guiotoku 10 Dez '19

10 Dez '19

Pessoal estou com uma dúvida quanto o pacote randomForest do R. Usei a função randomForest para modelar o meu conjunto de dados. Como irei aplicar o modelo já treinado para meu conjunto teste, salvei como .rda. Na próxima etapa iria subir o treino (.rda) para rodar o banco de dados independente. Porém, a forma como salvei está dando erro quando aplico a função predict. É como se o load ou a extensão .rda não conservasse o formato do banco de dados. E o erro é que o objeto "model" deveria ser um randomForest e ele fica como "character". Eu sei que eu poderia usar o que foi salvo na memoria do R, o modelCT, esse se usar no predict não daria erro. Mas para cada novo conjunto de dados independente, não gostaria de ficar treinando a todo o momento, gostaria de passar somente os parâmetros já treinados. Como um código executável para exemplificar o problema usei o da Iris. library(randomForest) modelCT <- randomForest::randomForest(Sepal.Length ~ ., data = iris, importance = TRUE) # modelagem com randomForest str(modelCT) save(modelCT, file = "model.rda") # salvei como rda model = load("model.rda") # Faço a chamada do arquivo treinado .rda str(model) predValid <- randomForest::predict(model, iris, type = "class") # o banco de dados deveria ser independente, mas aqui é somente para mostrar o erro do formato como está o model predValid <- predict(model, iris, type = "class")

2 1

Comparar elementos em 2 dataframes
by sérgio marques 06 Dez '19

06 Dez '19

Boa tarde, Alguma dica para iniciante no R, para se obter o seguinte. Tenho 2 tabelas (dataframes) com várias linhas e várias colunas, em cada célula dessas tabelas está um valor (cada valor é único, mas pode estar repetido ou não na tabela 2), quero obter uma 3.ª tabela (dataframe) onde conste cada valor da tabela 1 que não esteja presente na Tabela 2. Obrigado! Sérgio ________________________________________________________________________________

5 5

Re: [R-br] Validação cruzada
by sznelwar＠uol.com.br 05 Dez '19

05 Dez '19

4 4

Tamanho do efeito para Wilcoxon
by Elias Carvalho 05 Dez '19

05 Dez '19

Boa tarde Tenho duas amostras independentes com distribuição nao-normal e usei o método abaixo para ver a diferença entre elas: coin::wilcoxsign_test(data.amostra1 ~ data.amostra2, distribution="exact") Alguem pode me ajudar a calcular o tamanho do efeito? -- *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 2

Ajuda com download de csv pela internet
by Olympio Neto 04 Dez '19

04 Dez '19

Pessoal, bom dia. Estou com o seguinte problema. Preciso fazer automaticamente um download pela internet de um arquivo em CSV no site do TSE: http://www.tse.jus.br/eleitor/biometria/biometria-atual-do-dia, mas o link é gerado a partir de um comando em javascrit, algo parecido com *javascript:window.location.href=apex.server.url({p_request: 'FLOW_EXCEL_OUTPUT_R69876138832319545_pt-br'},2); null;* Alguém já teve esse problema? Anteriormente o site não limitava o número de registros e conseguia fazer o web scraping via html , mas a partir desta semana o site foi atualizado e está limitando a 30 por página. Alguma dica? Muito Obrigado Olympio Neto

2 3