Dear All,
I am relatively new to R but I am very interested in using some of the
available packages for astronomical applications.
While I have been able to understand and use other packages like cenken, I have
some difficulties with CensMixReg.
Here are the specifics (with info on R and package versions at the end of the
email), I would really appreciate your help.
1. First I tried to run the example in the CensMixReg.pdf file that accompanies
the download.
While I could reproduce the " left mixture censored Student-t model fit to the
data” (first example), the “left mixture censored Normal model fit to the data”
fails with the following error message:
Error in `row.names<-.data.frame`(`*tmp*`, value = value) :
invalid 'row.names' length
Is this issue documented somewhere?
2. Also in the CensMixReg.pdf file it is suggested to use the imm.fmcr function
to compute the uncertainties. However, the function does not seem to be part of
the CensMixReg or CensRegMod libraries (IFmixRegT <- imm.fmcr(y,x1,cc, fitT) ->
Error: could not find function "imm.fmcr"). I searched the R site but could not
find the imm.fmcr function. I am wondering if this is now osbolete and if there
is another function to compute uncertainties from the CensMixReg output.
3. Finally (and this is maybe the most important part), I am missing a working
example with clearly laid out inputs and outputs to use confidently the
CensMixReg package. Does anybody have such an example? More specifically, I was
not sure of the transformation of x into x1 in the CensMixReg example (where
there is an extra column of 1 in x1), is this to have the constant to the
fitted line stored in beta0? Are the other betas the slopes of respectively
age, educ and hours/1000? Where are the uncertainties on the slope and
intercept stored?
Thank you very much in advance.
With best wishes,
Ilaria
== Info on R and package versions:
> sessionInfo()
R version 3.3.0 (2016-05-03)
Platform: x86_64-apple-darwin13.4.0 (64-bit)
Running under: OS X 10.10.5 (Yosemite)
locale:
[1] C
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] CensRegMod_1.0 CensMixReg_0.7
loaded via a namespace (and not attached):
[1] tools_3.3.0 mvtnorm_1.0-5 mixsmsn_1.1-1
************************************
Dr. Ilaria Pascucci
Associate Professor
Department of Planetary Sciences
The University of Arizona
ilariapascucci.com
************************************
[[alternative HTML version deleted]]
______________________________________________
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.
