Re: [R] NANs in test Breslow-Day for svytable??

2012-09-02 Thread Diana Marcela Martinez Ruiz
12 12:34:11 +1200 > Subject: Re: [R] test Breslow-Day for svytable?? > From: tlum...@uw.edu > To: dwinsem...@comcast.net > CC: dianamm...@hotmail.com; r-help@r-project.org; m.steven...@massey.ac.nz > > On Sat, Sep 1, 2012 at 4:27 AM, David Winsemius > wrote: > > > > O

[R] test Breslow-Day for svytable??

2012-08-31 Thread Diana Marcela Martinez Ruiz
Hi all, I want to know how to perform the test Breslow-Day test for homogeneity of odds ratios (OR) stratified for svytable. This test is obtained with the following code: epi.2by2 (dat = daty, method = "case.control" conf.level = 0.95, units = 100, homogeneity = "breslow.day", verbose =

[R] sensitivity and specificity in svyglm??

2012-08-15 Thread Diana Marcela Martinez Ruiz
Hello, As obtained from a table svyglm clasificaion, sensitivity and specificity. The funtion ConfusionMatrix () of the library (caret) gives these results but not how to apply it to svyglm. thanks [[alternative HTML version deleted]] _

Re: [R] summary(svyglm) Pr (> | t |) ?

2012-08-02 Thread Diana Marcela Martinez Ruiz
will produce a Wald test for each > term in a model fit by svyglm(). > > I hope this helps, > John > > > -Original Message- > > From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] > > On Behalf Of Diana Marcela Martinez Ruiz > > S

[R] summary(svyglm) Pr (> | t |) ?

2012-08-02 Thread Diana Marcela Martinez Ruiz
Hello I want to know if the summary of the logistic model with survey Pr (> | t |) to test if the coefficient of the model is significant, ie is the p_valor wald test for the model coefficients, for I am interested to know if the three levels of the variable educational level are significant

[R] svyglm AIC: NA??????'

2012-07-25 Thread Diana Marcela Martinez Ruiz
hello I want to know why when I use the function "svyglm" for a logistic regression I get the AIC: NA. The code and the result is mestran below: mod2<-svyglm(APES_DICOT~Nivel_Educativo+Ocupacion_principal+Afiliacion_salud+Tiene_cuidador+Presencia_enfer_cronica+ Consumo_tabaco+Consumo_alcohol+

[R] (no subject)

2012-07-25 Thread Diana Marcela Martinez Ruiz
hello I want to know why when I use the function "svyglm" for a logistic regression I get the AIC: NA. The code and the result is mestran below: mod2<-svyglm(APES_DICOT~Nivel_Educativo+Ocupacion_principal+Afiliacion_salud+Tiene_cuidador+Presencia_enfer_cronica+ Consumo_tabaco+Consumo_alcohol+P

[R] stepwise in svyglm???

2012-07-24 Thread Diana Marcela Martinez Ruiz
Hello, I want to know how to perform stepwise elimination of variables to svyglm thanks [[alternative HTML version deleted]] __ R-help@r-project.org mailing list https://stat.ethz.ch/mailman/listinf

[R] test parallel slopes with svyolr

2012-07-07 Thread Diana Marcela Martinez Ruiz
Hello, I would like to know how to test the assumption of proportional odds or parallel lines or slopes for an ordinal logistic regression with svyolr Thanks [[alternative HTML version deleted]] __

[R] Help with multinomial regression to "survey"

2012-07-06 Thread Diana Marcela Martinez Ruiz
I would like to know which function to use for modeling a multinomial logistic regression model in the package "survey". thanks [[alternative HTML version deleted]] __ R-help@r-project.org mailing lis

[R] Help with the proportional odds model (Colombia)

2012-06-29 Thread Diana Marcela Martinez Ruiz
Good evening, I'm trying to mount the proportional odds model (ordinal polytomous regression model) and most of my independent variables are qualitative (and are transformed into dummy). I have an inconvenience that I would like to help me resolve, below I present part of the code I'm developi

[R] SVY: variance inflation factor VIF with complex survey

2012-06-28 Thread Diana Marcela Martinez Ruiz
Hello, Seeking a way to get the variance inflation factor VIF for a model of regression in complex survey, I have understood that without this package (SURVEY) RGui VIF obtained as follows: fit <- lm(mpg~disp+hp+wt+drat, data=mtcars) vif(fit) But I want to kn

[R] V-Cramer in svy

2012-06-10 Thread Diana Marcela Martinez Ruiz
Hello forum, I want to find in survey how to test V-cramer and coefficient Goodman-Kruskal gamma, to test the collinearity between varibles to enter a model. For the V-Cramer function is used with the package vcd assocstats but without considering the survey package. I appreciate the cooper

Re: [R] svychisq??

2012-05-31 Thread Diana Marcela Martinez Ruiz
Hello forum, I want to do a test of independence with svychisq, but I get an error, then this my code: am18 <- read.spss("C:/Users/diana/Dropbox/Semestre 10/Tesis 10/Tesis Diana/AMcomuna18-29MAR2012.sav", use.value.labels=TRUE, max.value.labels=Inf, to.data.frame=TRUE) b<-matrix(c(am