I'm a beginner trying my best to fit a multi-level model using complex survey data (American Community Survey microdata). My outcome is a binary variable at the person level, and my level two variable is metropolitan area (MSA). Ultimately, I hope to incorporate MSA-level predictors, but have not tried that yet.    The following code runs, but I get the error " ". I've gotten that error with a few specifications, which is why I'm using "FASTQUAD," "qpoints=1" and "inititer=10" > these changes were all based on other guidance I found online on how to address this error. Clear that "FASTQUAD" is ignored, though.   proc glimmix data=tmp2.usa1y2021_diss inititer=10 method=quadrature(fastquad qpoints=1)INITGLM empirical=classical; class MET2013; model doubledup (descending) =  /dist=binary cl solution  link=logit obsweight=PERWT; random intercept  / type=un subject=MET2013 weight=MET2013wgt; run;   WARNING: Sandwich covariance estimator is not available with the FASTQUAD suboption.The FASTQUAD  suboption is ignored.   NOTE: The EMPIRICAL option for METHOD=QUAD adjusts the covariance matrix of the fixed effects and the covariance parameters. It also affects the prediction covariance matrix for the random effectssolutions.   NOTE: The GLIMMIX procedure is modeling the probability that doubledup='1'. WARNING: The initial estimates did not yield a valid objective function.     --   I also found an example of someone using microdata with the following code, but this doesn't work for me. I thought it might be useful to go this route with something simpler, given the error I was receiving. This ran, but the results don't seem to be weighted to the population. Any help would be much appreciated!    proc glimmix data=tmp2.usa1y2021_diss; where MET2013 notin(0); weight perwt; class MET2013; model doubledup (desc)= /dist=binary link=logit ddfm=bw solution; random intercept/subject= MET2013 solution; run; ... View more

Hi, How can i add count rows by groups through proc sql AND THEN (this is the key thing I can't find instructions on) create a new "total" variable that is the same for each member of the group -- NOT just the Last member of the group. In conceptual terms, I want each group member to know how many members of its group have x=1.    Have dataset x y n 0 a 0 1 a 1 0 b 0 1 b 1 1 b 2 0 c 0 1 c 1 1 c 2 1 c 3   wanted output  x y n n2 0 a 0 1 1 a 1 1 0 b 0 2 1 b 1 2 1 b 2 2 0 c 0 3 1 c 1 3 1 c 2 3 1 c 3 3     Grouping variable is y, the variable to be counted/totaled is x.   i know  how to do it with data step: Proc sort data = have ;  by y; run;   data want; set have; by y; if first.y then n =0; n + x; run;     ... View more