Displaying 5 results from an estimated 5 matches for "ln_thc".
Did you mean:
en_tc
2017 Aug 10
4
PROC MIXED RANDOM equivalence in R nlme
...SAS v6.12 code is provided below:
------------------------------------------------------------------
/* READ DATA */
DATA emiss;
INFILE 'data.tab' LRECL=8000 FIRSTOBS=2 DLM='09'X MISSOVER DSD;
INPUT study $ vehicle $ thc rv t5 t9 ar ol ox su bz;
/* CREATE NEW VARIABLES */
ln_thc = log (thc);
new = study||vehicle;
/* PERFORM ANALYSIS */
PROC MIXED DATA=emiss MAXITER=1000 CONVH=1E-8 METHOD=REML NOCLPRINT
NOITPRINT;
CLASS new;
MODEL ln_thc = rv t5 t9 ar ol ox su bz
/S DDFM=RES;
RANDOM int rv t5 t9 ar ol ox su bz
/SUB=new;
RUN;
-...
2017 Aug 11
0
PROC MIXED RANDOM equivalence in R nlme
...gt; ------------------------------------------------------------------
> /* READ DATA */
> DATA emiss;
> INFILE 'data.tab' LRECL=8000 FIRSTOBS=2 DLM='09'X MISSOVER DSD;
> INPUT study $ vehicle $ thc rv t5 t9 ar ol ox su bz;
>
> /* CREATE NEW VARIABLES */
> ln_thc = log (thc);
> new = study||vehicle;
>
> /* PERFORM ANALYSIS */
> PROC MIXED DATA=emiss MAXITER=1000 CONVH=1E-8 METHOD=REML NOCLPRINT
> NOITPRINT;
> CLASS new;
>
> MODEL ln_thc = rv t5 t9 ar ol ox su bz
> /S DDFM=RES;
>
> RANDOM int rv t5...
2017 Aug 11
0
PROC MIXED RANDOM equivalence in R nlme
...ow:
lme not LME
data not DATA
random = RANDOM
> ------------------------------------------------------------------
>
> The R code I've devised for the PROC MIXED statement is shown below:
>
> ------------------------------------------------------------------
> FitTHC <- LME(ln_thc ~ rv + t5 + t9 + ar + ol + ox + su + bz,
> DATA = emiss,
> RANDOM = ??????? )
> ------------------------------------------------------------------
>
> As indicated, the problem I'm having is in constructing the equivalent
> code for the RANDOM and any rema...
2017 Aug 11
2
PROC MIXED RANDOM equivalence in R nlme
...> random = RANDOM
>> ------------------------------------------------------------------
>>
>> The R code I've devised for the PROC MIXED statement is shown below:
>>
>> ------------------------------------------------------------------
>> FitTHC <- LME(ln_thc ~ rv + t5 + t9 + ar + ol + ox + su + bz,
>> DATA = emiss,
>> RANDOM = ??????? )
>> ------------------------------------------------------------------
>>
>> As indicated, the problem I'm having is in constructing the
>> equivalent code for...
2017 Aug 11
0
PROC MIXED RANDOM equivalence in R nlme
...>> ------------------------------------------------------------------
>>>
>>> The R code I've devised for the PROC MIXED statement is shown below:
>>>
>>> ------------------------------------------------------------------
>>> FitTHC <- LME(ln_thc ~ rv + t5 + t9 + ar + ol + ox + su + bz,
>>> DATA = emiss,
>>> RANDOM = ??????? )
>>> ------------------------------------------------------------------
>>>
>>> As indicated, the problem I'm having is in constructing the equivalent...