Hi Robert:
Thanks for your response, as well.
I'm sorry. As you've discovered, I made some capitalization errors when
posting my R code. The actual R code does use "lme", "data", and "random".
The dataframe is indeed named "emiss" and each item in the formula is a
column in the dataframe. I used the following R code to read in the
comma-delimited file (first line contains headings):
emiss <- read.table("data.csv", header=TRUE, sep=",") # NA for
missing data (default)
A review of the dataframe "emiss" shows that all items, except "STUDY",
"VEHICLE" and "NEW", are of type "num". The items "STUDY", "VEHICLE",
and "NEW" are of type "factor".
I will work on creating some "fake data". However, since I don't have
SAS v6.12 anymore, I can't provide the so-called "correct results".
As for the R hang, your wording is correct: "R gets caught in a
processing loop that produces no errors or warnings", even after 15
minutes on an 8-core Mac Pro. As I recall, the SAS code produced an
answer in 15 seconds or less on a single core Mac II back in 1999.
Thanks again!
--Dennis
On 8/11/17 9:27 AM, Robert Baer wrote:
On 8/10/2017 8:34 AM, Dennis F. Kahlbaum wrote:
-- snip --
I don't have real help, but I'll remind you that R is case sensitive,
and it looks like that will be at least one problem in the solution
your are working on below:
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 remaining settings. I've tried
RANDOM = ~1 + rv + t5 + t9 + ar + ol + ox + su + bz | new)
but R hangs
Are the items in your random formula columns in a dataframe named
emiss? Do they have data types? Even if the data are proprietary
some fake data can make the problem more concrete.
You are saying "gets caught in a processing loop that produces no
errors or warnings"???
and never produces a result. Therefore, what is the equivalent code
for the SAS RANDOM?
Thanks!
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
______________________________________________
R-help@r-project.org mailing list -- To UNSUBSCRIBE and more, see
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
