Simple and multiple logistic regression

Question

I have almost two questions. I need a single covariate logistic regression (LR) for each of my variables. Should I do it manually in SPSS, selecting each variable and do logistic regression? Is there a "for each" cycle to do it? I should switch to R language to have what I want.

In the multivariables (multi covariates) LR, could I have missing values?

Thanks!!

Answer 1

If I understand you correctly, you want to fit two successive simple logistic regression model. I don't know if there's a specific instruction in SPSS that allows to switch the covariate of interest or cycle through them, but I guess you can run the two models in succession. In R, if your data are organized in a matrix or data.frame, this is easily done as

X <- replicate(2, rnorm(100))  # two random deviates
y <- rnorm(100)
apply(X, 2, function(x) lm(y ~ x))

About your second question, models like this are generally estimated using listwise deletion: any individuals having one or more missing observations on the covariates are deleted before estimating model parameters. Again, in R:

X[2,2] <- NA
summary(lm(y ~ X))

shows that one observation has been deleted, yielding 96 DF (instead of 97).

Answer 2

Here is some SPSS code to loop through each of your models.

*creating a simulated dataset.
input program.
loop #i = 1 to 100.
compute yvar = RV.BERNOULLI(.5).
compute xvar1 = RV.NORMAL(100,10).
compute xvar2 = RV.NORMAL(100,10).
compute xvar3 = RV.NORMAL(100,10).
compute listwise = RV.BERNOULLI(.1).
end case.
end loop.
end file.
end input program.
execute.

*using these variables you can run the regressions one equation at a time.
LOGISTIC REGRESSION VARIABLES YVAR
  /SELECT=LISTWISE EQ 0
  /METHOD=ENTER XVAR1.

*or you can create a macro and loop through your independent variable list.
DEFINE !logit_loop (dep = !TOKENS(1)
                             /ind = !CMDEND )
!DO !i !IN !ind
LOGISTIC REGRESSION VARIABLES !dep
  /SELECT=listwise = 0 
  /METHOD=ENTER !i.
!DOEND.
!ENDDEFINE.
!logit_loop dep = yvar ind = xvar1 xvar2 xvar3.

Note the loop does a selection on a variable named listwise = 0 , otherwise missing values for any variables would just be dropped (and potentially produce equations with different sets). While looping through your list is pretty simple in SPSS, where I think R is more convienant is the fact that the elements of your different objects are readily available in R (and so it is easier to access them and manipulate them). I'll update with an example later, but basically all the SPSS code does is produce alot of output which you have to read to find the values your interested in. In R you produce objects that have attributes you can extract. You can do the same thing in SPSS, but it is IMO more difficult and requires going into the output XML and parsing the info you want (which I don't know how to do).

So long story short if you are alright with just viewing the output this is easily accomplished in SPSS (or PASW now). If you want to produce a program to summarize specific elements of those models, I think it is easier to learn a solution in R.

Answer 3

SPSS has some scripting facilities (syntax, sax basic, pyhton, ...). I myself have so far only used syntax. Maybe the link below can help you to construct a loop that does the job, it points to the excellent website of the UCLA. The problem solved there is only slightly different from yours so chances are high that you can modify it according to your needs.

http://www.ats.ucla.edu/stat/spss/faq/looping_parallel_lists.html

psj