Respondendo à minha própria questão... m <- svyglm(formula = y ~ x, design = subset(d, complete.cases(x, y)), family = quasibinomial) sii <- svycontrast(m, quote((exp(`(Intercept)` + x) / (exp(`(Intercept)` + x) + 1)) - (exp(`(Intercept)`)/(exp(`(Intercept)`) + 1)))) Abraços, Leonardo Ferreira Fontenelle[1] Em Qua 10 dez. 2014, às 15:53, Leonardo Ferreira Fontenelle escreveu:
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[2]
Links: 1. http://lattes.cnpq.br/9234772336296638 2. http://lattes.cnpq.br/9234772336296638