- R-br - DINF UFPR

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

parallel no R
by Fernanda De Bastiani 17 Dez '14

17 Dez '14

Caros, alguém pode me indicar algum material (por ex. tutorial) que tenha exemplos de como fazer paralelização de função no R, além do manual do pacote parallel? Desde já muito obrigada,Fernanda

4 3

Mestrado em Modelagem Matemática e Computacional
by Wagner Tassinari 15 Dez '14

15 Dez '14

Gostaria de anunciar, que o Programa em Modelagem Matemática e Computacional (PPGMMC) está com as inscrições abertas para o Processo Seletivo de ingresso para o ano de 2015 . As inscrições são gratuitas e estão abertas do dia 19 de novembro de 2014 a 13 de janeiro de 2015. Conforme consta no edital, o endereço de inscrição é: Secretaria do PROPPG - UFRRJ – Campus Seropédica, sala 115, Prédio P1 BR 465 Km 7, Seropédica, RJ CEP.: 23897-000 Tel./Fax: (21) 26814707 As inscrições podem ser efetuadas pessoalmente (das 9:00 às 11:30 horas e das 13:30 às 16:00 horas) ou através de postagem. Quaisquer dúvidas ou informações os candidatos podem entrar em contato através dos emails: ppgmmc.ufrrj(a)gmail.com ppgmmc(a)ufrrj.br Site: Home Telefones: (21) 2681-4707 (21) 2681-4600 | | | | | | | | | HomeNotícias do MMC-UFRRJ | | | | Visualizar em r1.ufrrj.br | Visualizado por Yahoo | | | | |

1 0

Cálculo do deviance das variáveis explicativas em modelo de Poisson inflacionado de zeros
by ASANTOS 15 Dez '14

15 Dez '14

Bom dia Pessoal, Estou querendo saber como realizar o cálculo do deviance das variáveis explicativas em um modelo de Poisson inflacionado de zeros com a função zeroinfl() do pacote pscl, pois para um modelo do tipo glm eu faria uma anova(modelo, test="Chi"), mas um objeto criado pela função zeroinfl() não permite anova() e o teste de Vuong serve apenas para comparação entre modelos, existe alguma opção? Segue CRM: #------------------------------------------------------------------ # Definições da sessão. rm(list=ls()) require(pscl) #------------------------------------------------------------------ # Dados artificiais. da <- expand.grid(trat=gl(4,1), tempo=1:10) X <- model.matrix(~trat+tempo, da); ncol(X) betas <- c(0.1,0.9,0.6,0.3,0.7) eta <- X%*%betas y1 <- rpois(da$trat, lambda=exp(eta)) y2 <- rbinom(y1, size=1, prob=0.7) da$y <- y1*y2 str(da) #------------------------------------------------------------------ # Ajuste do modelo. m0 <- zeroinfl(y~trat+tempo|trat, data=da) summary(m0) #------------------------------------------------------------------ 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

Numeração de figuras Sweave na classe artigo
by Andre Oliveira 14 Dez '14

14 Dez '14

Pessoal,bom dia. Minhas figuras estão saindo assim. Figura 1Figura 2 Gostaria que fosse impressa assim. Figura 1.1Figura 1.2 Estou trabalhando assim. \SweaveOpts{width=15,height=12}\setkeys{Gin}{width=0.8\textwidth}\begin{figure}[H]\centering<<fig=T,keep.source=TRUE>>= par(mfrow=c(1,2))curve(dnorm(x,mean=1,sd=sqrt(2)),lwd=2,from=-6,to=17,ylab="Densidade")curve(dnorm(x,mean=10,sd=sqrt(2)),col=2,lwd=2,add=T)legend('topright',legend=c(expression(mu==1),expression(mu==10)),text.col=c(1,2),cex=1.5,bty="n")title("Comparando a distribuição normal com médias diferentes")curve(dnorm(x,mean=4,sd=sqrt(3)),lwd=2,from=-15,to=25,ylab="Densidade")curve(dnorm(x,mean=4,sd=sqrt(20)),col=2,lwd=2,add=T)legend('topright',legend=c(expression(sigma^2==3),expression(sigma^2==20)),text.col=c(1,2),cex=1.5,bty="n")title("Comparando a distribuição normal com variâncias diferentes")@\caption{Curva normal}\label{curvas}\end{figure} Tem como acertar isto? 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

1 0

Contastes não lineares
by Leonardo Ferreira Fontenelle 12 Dez '14

12 Dez '14

Boa tarde a todos! Procurei na web, mas ainda não achei a resposta para esta dúvida. Tenho um objeto *d* criado pela função *svydesign()*, do pacote *survey*. Duas variáveis desse objeto são *x* (uma variável contínua de 0 a 1) e *y* (uma variável dicotômica, 0 ou 1). Depois de ajustar o seguinte modelo: m <- svyglm(formula = y ~ x, design = subset(d, complete.cases(x, y)), family = quasibinomial) Eu gostaria de estimar a proporção (não as chances) de *y* quando *x* é igual a 0; o mesmo para quando *x* é igual a 1; e então calcular a diferença ou a razão entre essas proporções. Pelo que li, *svycontrast()* é a função mais flexível para estimar contrastes, inclusive não lineares. No entanto, ainda não consegui descobrir qual é a sintaxe para eu convencer o *svycontrast()* a estimar o que estou procurando. Alguém poderia me ajudar com isso? Obrigado! Leonardo Ferreira Fontenelle[1] Links: 1. http://lattes.cnpq.br/9234772336296638

1 1

Re: [R-br] Barplot múltiplo com dados de proporção de variáveis categóricas
by Ivan Bezerra Allaman 11 Dez '14

11 Dez '14

Há relação entre as variáveis? Se houver você pode fazer: tab <- ftable(d) mosaicplot(tab) O resto da confecção você tem capacidade de fazer! (s,f,p) Allaman

1 0

Barplot múltiplo com dados de proporção de variáveis categóricas
by ASANTOS 10 Dez '14

10 Dez '14

Bom dia Pessoal, Tenho 3 variáveis categóricas a, b e c, sendo a com três fatores, b com dois e c com quatro. Eu consigo representar a proporção delas em um barplot individualmente, porém gostaria de colocar as três variáveis em um mesmo barplot, isso é possível? Segue CRM: d <- data.frame(a=factor(sample(c('doente','sadio','indeterminado'), 100, rep=TRUE)), b=factor(sample(c('sim','nao'), 100, rep=TRUE)), c=factor(sample(c('ruim','regular','bom','otimo'), 100, rep=TRUE))) a.tb <- prop.table(table(d$a)) barplot(a.tb) b.tb <- prop.table(table(d$b)) barplot(b.tb) c.tb <- prop.table(table(d$c)) barplot(c.tb) 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

1 0