Novembro 2014 - 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

Estatísticas do R
by Andre Oliveira 26 Abr '16

26 Abr '16

Pessoal,preciso de estatísticas do R tais como nº de bibliotecas oficiais, pacotes que o R já vem na instalação. 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

8 8

Re: [R-br] Regressão logística ordinal e uso de pesos amostrais
by Luciane Maria Pilotto 23 Dez '14

23 Dez '14

Não consegui usar o comando svyolr, está dando os erros abaixo. # Utilizando as variáveis criadas q11 e q13 que foram transformadas em fator: > id$q11 <- as.factor(id$q11) > id$q13 <- as.factor(id$q13) > dstrat<-svydesign(ids=~1, probs=NULL, strata = NULL, variable = NULL, fpc=NULL, + data = NULL, nest = FALSE, weights=~peso,pps=FALSE) > dstrat Independent Sampling design (with replacement) svydesign(ids = ~1, probs = NULL, strata = NULL, variable = NULL, fpc = NULL, data = NULL, nest = FALSE, weights = ~peso, pps = FALSE) > svyolr(q13~q11,design=dstrat, method="logistic") Erro em svyolr.survey.design2(q13 ~ q11, design = dstrat, method = "logistic") : response must be a factor # Utilizando as variáveis quest_11 e quest_13 que são fatores no banco: > svyolr(quest_13~quest_11,design=dstrat, method="logistic") Erro em tapply(1:NROW(x), list(factor(strata)), function(index) { : arguments must have same length ========================================================== Eu tb poderia usar os pesos amostrais nas tabelas de contingencia, no entanto, tb não consigo usar os pesos nestas. Utilizei o comando crosstab (crosstab(q11, q13, peso, expected = TRUE,chisq = TRUE, plot = FALSE)). Dicas com outros comandos que utilizam os pesos amostrais são bem vindas!!! __________________________________________________ Luciane Maria Pilotto Mestre e Doutoranda em Saúde Bucal Coletiva - FO/UFRGS -------------------------------------------- Em ter, 25/11/14, Leonardo Ferreira Fontenelle <leonardof(a)leonardof.med.br> escreveu: Assunto: Re: [R-br] Regressão logística ordinal e uso de pesos amostrais Para: r-br(a)listas.c3sl.ufpr.br Data: Terça-feira, 25 de Novembro de 2014, 17:09 Não serve a http://r-survey.r-forge.r-project.org/survey/html/svyolr.html? Leonardo Ferreira Fontenelle Em Ter 25 nov. 2014, às 15:36, lutipilotto(a)yahoo.com.br escreveu: Olá, estou tentando rodar regressão logística ordinal e não estou conseguindo usar os pesos amostrais, está dando o seguinte erro: > m11 <- polr(q13 ~ q11 , method="logistic", weights = peso, data=id, Hess=TRUE) Aviso: glm.fit: fitted probabilities numerically 0 or 1 occurred Erro em optim(s0, fmin, gmin, method = "BFGS", ...) : valor inicial em vmmin não é finito *A variável q13 possui 5 ctegorias (1,2,3,4,5) e a variável q11 possui 4 categorias (1,2,3,4). Tem como solicitar o qui-quadrado, o número de observações e o pseudo R2 no mesmo comando utilizado para reg. log. ord.? Como utilizar o peso amostral nas tabelas de contingência? Estou usando os seguintes comandos: table(q11,q13) rowPercents(table(q11,q13)) ou CrossTable(q11, q13, digits=2, expected=TRUE, prop.r=TRUE, prop.c=TRUE, prop.t=FALSE,prop.chisq=FALSE, chisq = FALSE, format=c("SPSS"), dnn = c("local atend","satisfação")) Atenciosamente, Luciane M. Pilotto Doutoranda em Saúde Bucal Coletiva/UFRGS Enviado do Email do Windows _______________________________________________ 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. -----Anexo incorporado----- _______________________________________________ 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 15

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

Numero de dias acima de um valor
by Jefferson Ferreira-Ferreira 02 Dez '14

02 Dez '14

Olá amigos!! Estou com uma questão que não estou sabendo como resolver. Eu tenho uma série de dados de nivel d'água (3 vezes por dia). O que eu quero saber é em qual data começa a enchente. Como o nível sofre algumas "pequenas enchentes" que de fato não representam o real início da enchente, o que eu preciso é saber em qual data começa o seguinte evento: nível > 2cm por 45 dias consecutivos. Isso porque se o nível permanecer por menos de 45 dias consecutivos acima de 2cm eu não estou considerando isso como o início da enchente, mas como uma subida intermitente da água. Outro complicador é que preciso saber essa informação para cada um dos 10 equipamento (datalogger), discriminado na culona que chamei de logger. Ou seja, a data (em dia juliano) do início da enchente para cada logger. Em anexo envio os dados de apenas um dos meus aparelhos (logger == 1) Alguma ideia?

4 8

Raiz Quadrada do Erro Médio Quadrático e Teste de Médias
by Helder Gramacho 30 Nov '14

30 Nov '14

Olá Pessoal, boa noite! A seguir tem as explicações do problema e o código para reproduzir os resultados, procurei detalhar bem para facilitar a vida de quem se dispõe a ajudar. No trabalho disponível no link abaixo: http://www.lsie.unb.br/rbc/index.php/rbc/article/download/433/428 os autores validam o resultado de uma interpolação utilizando pontos de controle e calculando a Raiz Quadrada do Erro Médio Quadrático (REMQ) dado pela seguinte equação: [image: Imagem inline 2] Onde VMLTC é o Valor do Metro Linear de Testada Corrigida. Repeti o procedimento para 4 diferentes superfícies de interpolação e calculei o REMQ para cada uma como mostrado abaixo: #Download do arquivo no dropbox links <- c("https://www.dropbox.com/s/1grlssvxpihtcmt/valid.csv") tokens <- gsub("^.*/s/","",dirname(links)) fileNames <- basename(links) newLinks <- file.path("http://dl.dropbox.com/s", tokens, fileNames); newLinks for (a in newLinks) { tryCatch(download.file(a, dest=basename(a), mode='wb'), error=function(...) print("Falha no download!"))} valid<-read.table(file="valid.csv",sep=",",header=T,dec=".") # Função que calcula o REMQ (ve=valor estimado vr=valor real) remq<-function(ve,vr) { r<-ve-vr r2<-r^2 soma<-sum(r2) n<-length(r2) remq<-sqrt((1/(n-1))*soma) result<-c("Raiz do Erro Médio Quadrático"=remq) return(result) } # Calculando o REMQ remq.orig<-remq(valid[,2],valid[,1]) remq.mod.1.1<-remq(valid[,3],valid[,1]) remq.mod.2.1<-remq(valid[,4],valid[,1]) remq.mod.2.1.sar<-remq(valid[,5],valid[,1]) resumo.remq<-rbind(orig=remq.orig,mod.1.1=remq.mod.1.1,mod.2.1=remq.mod.2.1,mod.sar=remq.mod.2.1.sar);resumo.remq Agora gostaria de saber se as diferenças entre os REMQ são significativas estatisticamente e é aqui que está a minha dúvida, gostaria que avaliassem se o procedimento abaixo está correto, utilizo as diferenças entre os pontos de controle e os pontos obtidos nas superfícies para realizar a análise de variância e em seguida o Teste de Tukey. # Preparando para Teste de Médias m1<-valid[,2]-valid[,1] m2<-valid[,3]-valid[,1] m3<-valid[,4]-valid[,1] m4<-valid[,5]-valid[,1] dados.mod<-cbind(m1,m2,m3,m4);dados.mod dados.org <-data.frame(dif.mod=c(dados.mod[,1],dados.mod[,2],dados.mod[,3],dados.mod[,4]),Tipo=factor(c(rep("m1",20),rep("m2",20),rep("m3",20),rep("m4",20))),observ=factor(rep(1:20,4))) # Realizando a análise de variância tapply(dados.org$dif.mod,dados.org$Tipo,mean) ajuste<-aov(dados.org$dif.mod~dados.org$Tipo+dados.org$observ) summary(ajuste) posthoc<-TukeyHSD(x=ajuste,'dados.org$Tipo',conf.level=0.95);posthoc E então? Está correto utilizar as diferenças entre os pontos de controle e os pontos obtidos na superfície para realizar o teste de médias? Observe que os REMQ são bem diferentes, mas no teste médias as diferenças não são significativas, é aí que está a dúvida. Não sei se da forma como realizei o teste de médias está realmente avaliando as diferenças entre os REMQ. Desde já agradeço toda ajuda, *Hélder Gramacho * Recife-PE / *agrohelder(a)gmail.com <agrohelder(a)hotmail.com>*

2 2

Krigagem lognormal
by Helder Gramacho 27 Nov '14

27 Nov '14

Bom dia, Quando numa variável foi aplicada a transformação logarítmica como é possível realizar a krigagem lognormal no geo-r? e após a obtenção da estimativa como deve ser feita a transformação para a escala original? Alguém possui um exemplo de como isso pode ser realizado? Desde já muito obrigado, *Hélder Gramacho * Recife-PE / *agrohelder(a)gmail.com <agrohelder(a)hotmail.com>*

3 4

Arrumando um banco de dados
by Simone D. Sartorio 26 Nov '14

26 Nov '14

Pessoal, boa noite! Gostaria de uma ajuda, talvez seja bem fácil para aqueles que programam... 1) Tenho um banco de dados ("dados_prod") grande (110mil linhas) e preciso arrumar uma variavel que irei utilizar...a variavel é "dia_ciclo". O problema é que preciso q esta variavel (dia_ciclo), dentro de cada "n_injeção" esteja entre 1 a 14 (Ou seja, da injeção de n. 0 a 17, cada injeção tem 14 dias). O problema seria fácil se todas as minhas 240 vacas tivessem o mesmo número de tempos observados no inicio do experimento...mas não temos...alguma sugestão para resolver isso? 2) Deste mesmo banco de dados, um outro problema, seria retirar todas as injeções "n_injecao= -1". Estive pensando aqui, resolvendo o problema 2) primeiro o 1) fica bem facil... Alguém pode me ajudar nisso?Desde já, muito obrigadaSimone Exemplo: Dados de 1 vaca!> (vaca1731<- subset(dados_prod, dados_prod$n_vaca==1731)); dim(vaca1731) n_vaca n_lactacao trat pri_multiparas n_injecao dia_CICLO dia_exp producao saude1 1731 3 1 2 -1 -7 -21 46.0 healthy2 1731 3 1 2 -1 -6 -20 49.7 healthy3 1731 3 1 2 -1 -5 -19 50.2 healthy4 1731 3 1 2 -1 -4 -18 54.3 healthy5 1731 3 1 2 -1 -3 -17 47.4 heat6 1731 3 1 2 -1 -2 -16 49.3 healthy7 1731 3 1 2 -1 -1 -15 46.7 healthy8 1731 3 1 2 -1 0 -14 45.1 healthy9 1731 3 1 2 0 -13 -13 47.8 healthy10 1731 3 1 2 0 -12 -12 45.5 healthy11 1731 3 1 2 0 -11 -11 45.2 healthy12 1731 3 1 2 0 -10 -10 47.4 enrollment13 1731 3 1 2 0 -9 -9 48.3 enrollment14 1731 3 1 2 0 -8 -8 46.4 enrollment15 1731 3 1 2 0 -7 -7 47.9 enrollment16 1731 3 1 2 0 -6 -6 50.3 enrollment17 1731 3 1 2 0 -5 -5 50.8 enrollment18 1731 3 1 2 0 -4 -4 47.9 enrollment19 1731 3 1 2 0 -3 -3 46.2 heat20 1731 3 1 2 0 -2 -2 50.8 healthy21 1731 3 1 2 0 -1 -1 45.8 healthy22 1731 3 1 2 0 0 0 50.5 healthy23 1731 3 1 2 1 1 1 49.0 healthy24 1731 3 1 2 1 2 2 48.8 healthy25 1731 3 1 2 1 3 3 53.9 healthy26 1731 3 1 2 1 4 4 45.4 healthy27 1731 3 1 2 1 5 5 38.6 healthy28 1731 3 1 2 1 6 6 41.3 healthy29 1731 3 1 2 1 7 7 50.2 healthy30 1731 3 1 2 1 8 8 50.2 healthy31 1731 3 1 2 1 9 9 50.4 healthy32 1731 3 1 2 1 10 10 49.9 healthy33 1731 3 1 2 1 11 11 50.1 healthy34 1731 3 1 2 1 12 12 52.0 healthy35 1731 3 1 2 1 13 13 47.8 healthy36 1731 3 1 2 1 0 14 50.0 healthy37 1731 3 1 2 2 1 15 49.7 healthy38 1731 3 1 2 2 2 16 50.4 healthy39 1731 3 1 2 2 3 17 48.9 healthy40 1731 3 1 2 2 4 18 49.5 healthy41 1731 3 1 2 2 5 19 48.8 healthy42 1731 3 1 2 2 6 20 56.3 healthy43 1731 3 1 2 2 7 21 52.1 healthy44 1731 3 1 2 2 8 22 51.6 healthy45 1731 3 1 2 2 9 23 50.5 healthy46 1731 3 1 2 2 10 24 46.0 healthy47 1731 3 1 2 2 11 25 53.7 healthy48 1731 3 1 2 2 12 26 51.7 healthy49 1731 3 1 2 2 13 27 49.1 healthy50 1731 3 1 2 2 0 28 51.9 healthy51 1731 3 1 2 3 1 29 51.4 healthy52 1731 3 1 2 3 2 30 43.6 healthy...259 1731 3 1 2 17 13 237 19.9 lameness260 1731 3 1 2 17 0 238 18.8 lameness[1] 259 9 ***---***---***---***---***---***---***---***---***---***---***---***---**** Simone Daniela Sartorio* Professora Adjunta I no CCA/UFSCar, Araras/SP.* Doutora e Mestre em Estatística e Experimentação Agronômica - ESALQ/USP.* Licenciada em Matemática - UNESP/Rio Claro.* Contatos: * tel: (19)3543-2614, ramal 2664* cel: (19)98182-0586 - Tim***---***---***---***---***---***---***---***---***---***---***---***---*** Tenha um bom dia! ;)

2 1

(sem assunto)
by Thaís Brenda Martins 26 Nov '14

26 Nov '14

Boa noite prezad@s, Desenvolvi uma interface gráfica no R com a biblioteca Tcl/Tk. Gostaria que o usuário inserisse os elementos de uma amostra populacional na caixa de entrada (tkentry) ou importasse esses dados de um txt por exemplo. A primeira parte já fiz, funciona corretamente. Mas a segunda não consegui. O que tentei fazer foi: - criei um botão pra que o usuário escolha o arquivo. - Na escopo da função desse botão precisaria jogar os valores do arquivo em forma de uma string na caixa de entrada (onde o usuário digita os elementos quando não quer importar). Para quando o usuário clicar no botão de realizar teste não seja impedido de prosseguir porque o campo esta vazio? Alguém que trabalhe com está biblioteca poderia me dar uma ajuda? Att., Thaís Brenda Martins

3 3