Dezembro 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

Discretizando uma variável contínua no chisq.test
by ASANTOS 07 Jan '15

07 Jan '15

Caros Listeiros, Gostaria de comparar valores entre parcelas observadas e estimadas para saber a similaridade ou não entre dados estimados vs observados. No entanto, as variáveis observadas são dados de contagem e os dados estimados são variáveis contínuas, então minha dúvida é se violo algum pressuposto estatístico ao discretizar os dados estimados com a função round() para realizar o teste de Chi-quadrado? Segue CRM: ## Comparação entre parcelas estimadas e observadas # #Dados das parcelas observadas - Dados de contagem data.obs<-c(11, 8, 8, 8, 7, 14) # #Medias estimadas para as parcelas observadas - Dados contínuos data.est<-c(14.6,16.0, 19.1, 3.4, 9.4, 15.8) # #Distretizando as médias estimadas tr.data.est<-round(data.est) # #Comparando parcelas observadas e estimadas chisq.test(x=tr.data.est, y=data.obs, simulate.p.value=TRUE) Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ====================================================================== --- Este email foi escaneado pelo Avast antivírus. http://www.avast.com

2 3

ANOVA medidas repetidas no R
by Roberto Guevara 03 Jan '15

03 Jan '15

Caros usuários de R, minha dúvida na verdade é mais sobre estatística do que o R em si. Estou utilizando uma ANOVA de medidas repetidas com três fatores sendo que há repetição somente para dois destes fatores. O código e o resultado que estou obtendo segue o modelo abaixo: > summary(aov(dependente ~ fator1 * fator2 * fator3 + + Error(individuo/(fator2 * fator3)))) Error: individuo Df Sum Sq Mean Sq F value Pr(>F) fator1 1 80616 80616 0.260 0.645 fator2 1 809787 809787 2.611 0.205 Residuals 3 930281 310094 Error: individuo:fator2 Df Sum Sq Mean Sq F value Pr(>F) fator2 2 181649 90824 1.491 0.282 fator1:fator2 2 186811 93405 1.533 0.273 Residuals 8 487309 60914 Error: individuo:fator3 Df Sum Sq Mean Sq F value Pr(>F) fator3 2 474691 237345 14.526 0.00502 ** fator1:fator3 2 138023 69012 4.224 0.07163 . fator2:fator3 2 184239 92120 5.638 0.04189 * Residuals 6 98035 16339 --- Signif. codes: 0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1 Error: individuo:fator2:fator3 Df Sum Sq Mean Sq F value Pr(>F) fator2:fator3 4 58378 14594 0.374 0.824 fator1:fator2:fator3 4 198072 49518 1.270 0.323 Residuals 16 624020 39001 Error: Within Df Sum Sq Mean Sq F value Pr(>F) Residuals 51 2676429 52479 Todas as cinco tabelas são imprescindíveis para fins de avaliar o efeito dos fatores sobre a variável dependente? Como diferencio o que cada uma expõe? Por exemplo, o resultado da primeira tabela indica o efeito dos fatores 1 e 2 isoladamente ao retirar o efeito dos indivíduos? Qual a diferença entre os resultados expostos para o fator2 na primeira e na segunda tabela? Desde já, muito obrigado pela atenção. -- Roberto Guevara Ferreira Lima Biólogo - CRBio 6ª Região: 73146/06-D Mestre em Zoologia Laboratório de Ecologia e Zoologia de Vertebrados Instituto de Ciências Biológicas Universidade Federal do Pará Rua Augusto Corrêa, nº 01 - Bairro: Guamá. CEP: 66075-110 Currículo Lattes: http://lattes.cnpq.br/7363382159007600

1 1

Espaço entre texto e imagens - Sweave
by Andre Oliveira 30 Dez '14

30 Dez '14

Pessoas, boa noite. Depois de muita pesquisas eu recorro a vocês. O que posso colocar no preâmbulo de meu .rnw para reduzir os espaços entre imagens e texto? obrigado 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

4 7

Selecionando níveis dentro de fator
by ASANTOS 29 Dez '14

29 Dez '14

Caros listeiros, Tenho um tratamento constituído de 26 letras e gostaria de criar um vetor arbitrário para separar alguns níveis dentro deste tratamento, porém não quero utilizar a função & para cada nível (pois no meu banco de dados original ficaria muito longo), mas especificar diretamente um lista com eles sendo: samp.niv<-c("a","k","o","p","r","s","t")##Níveis escolhidos dentro de trat Porém, segundo o meu CRM se faço uma relação de igualdade dentro do banco de dados, no meu caso da, não funciona, alguém poderia me ajudar, segue CRM: ##Criando um banco artificial rm(list=ls(all=TRUE)) da <- expand.grid(trat=rep(paste(letters[1:26]),100), tempo=1:10) da$y1 <- rnorm(da$trat) samp.niv<-c("a","k","o","p","r","s","t")##Níveis escolhidos dentro de trat da.1<-da[da[,1]==samp.niv,]## Selecionando os níveis de interesse Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ====================================================================== --- Este email foi escaneado pelo Avast antivírus. http://www.avast.com

2 2

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

22 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

Problema em biblioteca do R para uso com Rstudio no Ubuntu
by ASANTOS 20 Dez '14

20 Dez '14

Caros listeiros, Estou usando o R version 2.15.0 (2012-03-30) no Ubuntu em função de pacotes pessoais e outros não atualizados para as versões 3.x, no entanto, ao tentar instalar o Rstudio para utilizar o R, esse não funciona com o erro R shared library (/usr/local/lib/R/lib/libR.so) not found. Minha biblioteca do R esta em /usr/local/bin/R e gostaria de saber o que devo fazer para o Rstudio funcionar, Obrigado, -- ====================================================================== Alexandre dos Santos Proteção Florestal IFMT - Instituto Federal de Educação, Ciência e Tecnologia de Mato Grosso Campus Cáceres Caixa Postal 244 Avenida dos Ramires, s/n Bairro: Distrito Industrial Cáceres - MT CEP: 78.200-000 Fone: (+55) 65 8132-8112 (TIM) (+55) 65 9686-6970 (VIVO) e-mails:alexandresantosbr@yahoo.com.br alexandre.santos(a)cas.ifmt.edu.br Lattes: http://lattes.cnpq.br/1360403201088680 ====================================================================== --- Este email foi escaneado pelo Avast antivírus. http://www.avast.com

3 3

simetria e curtose, limites - distribuição normal
by Luiz Roberto Martins Pinto 17 Dez '14

17 Dez '14

Caros, Considerem x=rnorm(1000). As funções skewness(x) # retorna valor simetria. kurtosis(x) # retorna valor curtose. A distribuição normal teórica é simetrica, com skewness(x)=0, e mesocúrtica, com kurtosis(x)=0. Preciso saber os intervalos para skewness(x) e kurtosis(x), de modo que o vetor x seja considerado 'normalmente distribuído'. Luiz Roberto. Luiz Roberto Martins Pinto Prof. Pleno/DCET/UESC Laboratório de Estatística Computacional Universidade Estadual de Santa Cruz Ilhéus-Bahia luizroberto.uesc(a)gmail.com skype: lrmpinto http://lattes.cnpq.br/2732314327604831

2 3