search for: mcculloch

Hello, I'm just beginning to learn R. What books/online learning modules with datasets would you suggest? Thank you! Best wishes, Michael ____________________________________ Michael McCulloch Pine Street Clinic Pine Street Foundation 124 Pine Street, San Anselmo, CA 94960-2674 tel 415.407.1357 fax 415.485.1065 email: mm at pinest.org web: www.pinest.org www.pinestreetfoundation.org www.medepi.net/meta

...al reviewer says this is not any use. I cannot afford to pay for colour prints. Any ideas on what is the best way to distinguish 10 groups on an xy scatter plot? If all else fails?I can just remove the graph and give them a table of regression coefficients. Thanks. Yours Sincerely Andrew McCulloch

Hi , I have a data frame (Data1) which has a area variable (AREA) with codes (RA,RG,LT etc.). I would like to split the data frame into individual data frames called RA, RG, LT conatining only the observations from those areas. Thanks yours sincerely Andrew McCulloch -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send "info", "help", or "[un]subscribe" (in the "body", not the subject !) To: r-help-request at stat...

...I had been using does not have functions for doing a weighted Cox analysis. Any advice pointing me in the direction of resources to help probability-of-treatment-weighted Cox proportional hazards would be most appreciated. Best wishes, Michael ____________________________________ Michael McCulloch Pine Street Clinic Pine Street Foundation 124 Pine Street, San Anselmo, CA 94960-2674 tel 415.407.1357 fax 415.485.1065 email: mm at pinest.org web: www.pinest.org www.pinestreetfoundation.org www.medepi.net/meta

Hi all, sorry for my english, but I don't speak yours language. I'm trying to use bart() and rbf(). The package I'm using now is "BayesTree" and "neural", respectively. I could create the models, but I can't predict my test data. Does anyone have such an experience? Any advice is appreciated! Thank you in advanced!. Andr? -- View this message in

hi All, thanks for the responses. After being dropped into the # Filesystem repair prompt, ( on account of "inode 27344909 has illegal blocks" ) following warm reboot (via "reboot") after finding (SAN ) filesystem in read-only mode yesterday morning (possibly because of HBA fault on SAN) , I ran fsck -r /data (Linux version 2.6.18-92.1.18.el5 , Red Hat 4.1.2-42 , ext3

...in the R community know of an R function (or S-plus or Mathematica or Matlab or C++ or even Fortran) to do this. (Note that bkde2D, ash2 etc will not work with binned data input and using 4 x 10^7 raw samples on a 0(n^2) algorithm is probably not feasible) Thanks in advance. Regards, James McCulloch Post-Doc Fellow University of Technology Sydney DISCLAIMER\ ================================================...{{dropped}}

...trom uses the "full maximum likelihood" by approximating the integral via Gauss- Hermite quadrature. However, glmmML is only valid for the random intercept model and the binomial family must be represented only as binary data. Does the lmer do the work? The reason for asking that is: McCulloch (1997): "Unfortunately, PQL methods can perform poorly relative to Maximum likelihood" Agresti (2002) p. 524 "........but where possible we recommend using ML rather than PQL". Thanks a lot, Abderrahim [[alternative HTML version deleted]]

hi All, Re below - thanks for the replies to the fsck question. I have run fsck -r /data and this eventually completed with several changes made. Now I am getting the following error from the qla2xxx driver on boot . . . qla2xxx 0000:05:0d.0: Mailbox command timeout occured. Scheduling ISP abort. qla2xxx 0000:05:0d.0: Mailbox command timeout occured. Scheduling ISP abort. . etc . It looks

As I understand it lme (in R v1.9.x) estimates random effect variances on a log scale, constraining them to be positive. Whilst this seems sensible, it does lead to apparently biased estimates if the variance is actually zero - which makes our simulation results look strange. Whilst we need to think a bit deeper about it - I still haven't got my head around what a negative variance could

Hi, I am wondering whether anyone of you already implemented in R or S+ the methods of estimation the parameters of stable distribution. If there is, i will be happy to get the copy of the codes. Regards Dedi

Hi, I have very simple balanced randomized block design where I total have 48 observations of a measure of weights of a product, the product was manufactured at 4 sites, so each site has 12 observations. I want to use lme() from nlme package to estimate the standard error of the product weight. So the data look like: MW site 1 54031 1 2 55286 1 3 54396 2 4 52327 2 5 55963

hi all We have a rh linux server connected to two HP SAN controllers, one an HSV200 (on the way out), the other an HSV400 (on the way in). (Via a Qlogic HBA). /etc/multipath.conf contains this : device { vendor "(COMPAQ|HP)" product "HSV1[01]1|HSV2[01]0|HSV300|HSV4[05]0" getuid_callout "/sbin/scsi_id -g -u -s /block/%n"

CRAN (and crantastic) updates this week New packages ------------ * bdoc (1.0) Michael Anderson http://crantastic.org/packages/bdoc This package contains a function that will classify DNA barcodes as well as a few test and reference data sets. * bdsmatrix (1.0) Terry Therneau http://crantastic.org/packages/bdsmatrix This is a special case of sparse matrices, used by coxme and

hi All, A fault on our SAN dropped us down to a read-only filesystem and after reboot, we have an "Unexpected Inconsistency" and I am being instructed by the boot to run fsck manually without -a or -p (this was after I think processing around 15% of the filesystem) The specific message is "inode 27344909 has illegal blocks" I recall running fsck some years ago on smaller

The reason for the behaviour observed below turned out to be that the device entry in /etc/multipath.conf was inadvertently appended *after* the devices section , rather than inside it - so that we had #devices { # device { # blah blah # } (file has a bunch of defaults commented out) # etc #} # # device { our settings } *rather than*