Outubro 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

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

Desdobrar interação dupla e tripla com glm
by Maurício Lordêlo 07 Nov '14

07 Nov '14

Saudações a todos! No mês passado, solicitei um auxílio sobre como analisar os dados de um experimento realizado com explantes de bananeiras considerando três fatores. Na ocasião descrevi o problema mais ou menos da seguinte forma (ocorreram algumas mudanças): "*Experimento realizado com explantes de bananeiras considerando três fatores:* *1. Genótipo, com dois níveis;* *2. Meio, com sete níveis;* *3. Tipo de corte, com oito níveis.* *O objetivo é verificar qual tipo de corte induziu melhor o explante na formação do embrião considerando os setes meios e os dois genótipos.* *Para isso foram preparadas duas placas (caracterizando o número de repetições) com uma quantidade de explantes idênticos de acordo com o corte e o genótipo (esta quantidade de explantes em cada placa variou de 6 a 11). Cada placa com estes explantes foi introduzida nos sete diferentes meios. As respostas foram obtidas registrando-se, em relação ao total de cada placa, quantos explantes oxidaram, quantos formaram calos, quantos não desenvolveram, etc.* *Os dados com uma variável resposta está neste link:* *https://www.dropbox.com/s/ci4u8c4vn4cmk4l/dados.txt?dl=0 <https://www.dropbox.com/s/ci4u8c4vn4cmk4l/dados.txt?dl=0> *" Obtive retorno do Walmes que contribui com valiosas sugestões. *O que preciso agora é realizar o desdobramento de interações duplas e triplas utilizando "glm"*. No CRM abaixo, apenas a interação dupla foi significativa. Porém, tenho outras variáveis neste experimento cuja interação tripla foi significativa e será necessário também desdobrar. Na busca que fiz, encontrei apenas comandos de desdobramentos para modelos ajustados com "aov". Antecipadamente, agradeço quem possa contribuir. Maurício dados = read.table("dados.txt",header=TRUE, sep="\t") str(dados) attach(dados) # Modelo fatorial completo, glm (quasi)binomial para variável resposta "n_ex_ox" model1 = glm(cbind(good=n_ex_ox, bad=n_explant-n_ex_ox)~genotipo*meio*tipo, data=dados, family=quasibinomial) anova(model1,test="F") summary(model1)

2 4

Dúvida
by Vanessa Felix 04 Nov '14

04 Nov '14

Alguém poderia me ajudar a passar esse comando de SAS pra R, tem que usar mtos if's e não tô conseguindo acertar. SELECT; WHEN(V0602=2) DO; *(FREQUENTA A ESCOLA); SELECT; WHEN(V6003=1) DO; *(REGULAR DE PRIMEIRO GRAU); SELECT; WHEN (V6030=1) DO;*(DURAÇÃO DO ENSINO FUNDAMENTAL = 8 ANOS); SELECT; WHEN(V0605=1) EDUCA3=0; WHEN(V0605=2) EDUCA3=1; WHEN(V0605=3) EDUCA3=2; WHEN(V0605=4) EDUCA3=3; WHEN(V0605=5) EDUCA3=4; WHEN(V0605=6) EDUCA3=5; WHEN(V0605=7) EDUCA3=6; WHEN(V0605=8) EDUCA3=7; OTHERWISE EDUCA3=100; Desde já obrigada, -- ------------------------------------------------------------------------------------------ Vanessa Felix do Nascimento Sergio Graduada em Estatística - ICE/UFJF Mestranda em Saúde Coletiva - IMS/UERJ Assistente de pesquisa I - Instituto de Pesquisa Econômica Aplicada – IPEA (32) 8854-4881 / 9142-8428 (21) 97656-3448 -------------------------------------------------------------------------------------------

3 4

Legenda mapa erro padrão predição krigagem
by Helder Gramacho 31 Out '14

31 Out '14

Olá pessoal boa noite, Na figura estão representados o erro padrão de predição de duas krigagens obtidas a partir dos seguintes comandos: # mapa da esquerda image(kr.ln.vti.ha,val=sqrt(kr.ln.vti.ha$krige.var),ylim=c(-10,110),xlim=c(1,80),col=rev(heat.colors(10)),ylab="Coordenada N",xlab="Coordenada E") legend.krige(x.leg=c(0,80),y.leg=c(-8,0),val=sqrt(kr.ln.vti.ha$krige.var),col=rev(heat.colors(10))) # mapa da direita image(kr.vp.mod.2.1.sar,val=sqrt(kr.vp.mod.2.1.sar$krige.var),ylim=c(-10,110),xlim=c(1,80),col=rev(heat.colors(10)),ylab="Coordenada N",xlab="Coordenada E") legend.krige(x.leg=c(0.81,83.74),y.leg=c(-8,0),val=sqrt(kr.vp.mod.2.1.sar$krige.var),col=rev(heat.colors(10))) Como dá para perceber, os intervalos de representação das cores são diferentes, como o objetivo é compará-las o ideal é que o intervalo fosse o mesmo para ambas. É possível fixar o intervalo de representação para geração dos dois mapas? Como? [image: Imagem inline 1] Desde já agradeço a ajuda, *Hélder Gramacho * Recife-PE / *agrohelder(a)gmail.com <agrohelder(a)hotmail.com>*

1 0

Teste de hipoteses para assimetria e curtose.
by Andre Oliveira 30 Out '14

30 Out '14

Pessoas,estou buscando como testar as hipóteses para os coeficientes de assimetria e curtose. Aqui eu consigo os coeficientes mas não tenho o teste. require(e1071) x=c(72, 70, 77, 60, 67, 69, 68, 66, 65, 71, 69) kurtosis(x,type=1)+3 skewness(x,type=1) Para assimetria gostaria de testar se o coeficiente B1 difere de 0 (Um teste de normalidade o suficiente?) Para curtose preciso testar se B2 difere de 3. Alguém tem alguma sugestões? E só mesmo para colocar em minha notas de aulas. 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

2 1

Instalação plotKML
by Felinto COSTA 30 Out '14

30 Out '14

Colegas. Não consigo executar o pacote <plotKML>. O console retorna a mensagem: ------------------------------------------------------------------------ > library(plotKML) Error in loadNamespace(j <- i[[1L]], c(lib.loc, .libPaths()), versionCheck = vI[[j]]) : there is no package called ‘Rcpp’ Erro: package or namespace load failed for ‘plotKML’ ------------------------------------------------------------------------ Essa é a versão que tenho do R: ------------------------------------------------------------------------ > sessionInfo() R version 3.1.1 (2014-07-10) Platform: x86_64-w64-mingw32/x64 (64-bit) locale: [1] LC_COLLATE=Portuguese_Brazil.1252 LC_CTYPE=Portuguese_Brazil.1252 [3] LC_MONETARY=Portuguese_Brazil.1252 LC_NUMERIC=C [5] LC_TIME=Portuguese_Brazil.1252 attached base packages: [1] stats graphics grDevices utils datasets methods base ------------------------------------------------------------------------ Já tentei a instalação do mesmo com <install.packages("plotKML", dependencies=T)> sem dar certo também. Mesmo problema ao tentar executar o pacote <ggplot2>. Agradeço a cooperação.

2 2

Função qnorm()
by Jodavid Ferreira 29 Out '14

29 Out '14

Alguém sabe como funciona a função qnorm()? estou precisando implementar alguém sabe a fórmula utlizada? em que eu tenho o valor da probabilidade e me retorne o valor de Z. por exemplo em dou a entrada de 0.5 e me retorne 0 ou dou entrada de 0.05 e me retorne -1.64 -- Jodavid Ferreira Bacharelando em Estatística - Universidade Federal da Paraíba - www.de.ufpb.br Estagiário no Laboratório de Estatística Aplicada ao Processamento de Imagens e Geoprocessamento - http://de.ufpb.br/~leapig Site Pessoal - http://blog.jodavidferreira.com Twitter - @JodavidFerreira <http://twitter.com/JodavidFerreira>

2 1