similar to: Schafer's CAT for MI

Displaying 20 results from an estimated 2000 matches similar to: "Schafer's CAT for MI"

2003 Apr 01
1
Shafer's MI software for S-plus
Greetings folks, Shafer's S-plus package "norm" for multiple imputation of missing values in multivariate normal data has been most kindly and usefully ported to R by Alvaro A. Novo. Shafer's website http://www.stat.psu.edu/~jls/ lists four S-plus packages in all: NORM - multiple imputation of multivariate continuous data CAT - multiple imputation of multivariate
2003 Apr 29
1
Shafer's MIX: Query on code
Thanks to Fernando Tusell and especially to Brian Ripley for their work on 'mix', leading to an apparently good package mow available on CRAN. Going through the R code for the function prelim.mix, I am wondering why the following method of calculation is used at one point: umd <- as.integer(round(exp(cumsum(log(d))))) (d is a vector containing, in effect, the numbers of levels of
2003 Sep 20
3
conditional function definition?
Hi Folks, What is the best way to avoid a function being read in anew (and masking an exiting function) when a definition of it has already been established in R? Reason: Fernando Tusell and I are working up Schafer's 'CAT' for R (basically done now, just needs some cosmetic tidying up). This uses a function 'slice.index', present in S but not in the versions of R we were
2002 Apr 08
4
Missing data and Imputation
Hi Folks, I'm currently looking at missing data/imputation methods (including multiple imputation). S-Plus has a "missing data library". What similar resources are available within R? Or does one roll one's own? Best wishes to all, Ted. -------------------------------------------------------------------- E-Mail: (Ted Harding) <Ted.Harding at nessie.mcc.ac.uk>
2006 Mar 18
0
No subject
To estimate the covariance matrix of e you could use the sample covariance matrix of the residuals. If desired, use its cholesky decomposition to transform to make the error approximately uncorrelated, then refit (and back-transform the coefficient matrix). Stacking the columns of Y and replicating X won't do what you write; it forces each univariate regression to have the same coefficients.
2003 Oct 20
0
aliases
How about: nis.na <- complete.cases --- From: <Ted.Harding at nessie.mcc.ac.uk> Hi Folks, My recent response to Laura Quinn's query about matrix subsetting reminded of a question. I wrote: iDir <- ((Winds[,20]<45)|(Winds[,20]>315))&(!is.na(Winds[,20])) Now, I find "!is.na" a bit awkward to type, so I might prefer to type it as "nis.na".
2004 Dec 06
0
What is the most useful way to detect nonlinearity in lo
> -----Original Message----- > From: r-help-bounces at stat.math.ethz.ch > [mailto:r-help-bounces at stat.math.ethz.ch] On Behalf Of > Ted.Harding at nessie.mcc.ac.uk > Sent: Sunday, December 05, 2004 7:14 PM > To: r-help at stat.math.ethz.ch > Subject: Re: [R] What is the most useful way to detect > nonlinearity in lo > > > On 05-Dec-04 Peter Dalgaard wrote:
2006 Jul 21
0
connection to X11 problem: problem fixed
Hi, I finally managed to make it work (I just needed to have a X11 window open). Thank you very much for your help. Agnes -----Original Message----- From: Ted Harding [mailto:Ted.Harding at nessie.mcc.ac.uk] Sent: Friday, July 21, 2006 3:34 PM To: Paquet, Agnes Cc: r-help at stat.math.ethz.ch Subject: RE: [R] connection to X11 problem On 21-Jul-06 Paquet, Agnes wrote: > Dear List, >
2004 Apr 02
1
tan(mu) link in GLM
Hi Folks, I am interested in extending the repertoire of link functions in glm(Y~X, family=binomial(link=...)) to include a "tan" link: eta = (4/pi)*tan(mu) i.e. this link bears the same relation to the Cauchy distribution as the probit link bears to the Gaussian. I'm interested in sage advice about this from people who know their way aroung glm. >From the surface, it looks
2003 Sep 16
1
"Old" libraries with new R?
Hi Folks, I'm currently installing R-1.7.1 off CRAN. As it happens, I have a CD (kindly made for me by Linux Emporium) containing all the libraries which were on CRAN early this year when I installed R-1.6.1. This is highly convenient, since the alternative would be several hours on-line. While a recent library will on installation announce the fact should it need a newer version of R than
2003 Jul 06
1
Conditional Distribution of MVN variates
Hi Folks, Given k RVs with MVN distribution N(mu,S) (S a kxk covariance matrix), let (w.l.o.g.) X1 denote the first r of them, and X2 the last (k-r). Likewise, let mu1 and mu2 denote their respective expectations. Then, of course, the expectation of X2 given X1=x1 is mu2 + S21*inv(S22)*(x1 - mu1) and the covariance matrix of X2 given X1=x2 is S22 - S21*inv(X11)*S12 where Sij is the
2002 Dec 20
1
Printing correlation matrices (lm/glm)
Hi Folks, I'm analysing some data which, in its simplest aspect, has 3 factors A, B, C each at 2 levels. If I do lm1 <- lm(y ~ A*B) say, and then summary(lm1, corr=T) I get the correlation matrix of the estimated coeffcients with numerical values for the correlations (3 coeffs in this case). Likewise with 'glm' instead of 'lm'. However, if I do lm2 <- lm(y ~
2003 Jan 16
1
Multivariate regression in R
Hi Folks, I want to do multivariate regression in R, i.e. basically (but with a complication -- see below): given an Nxp matrix Y of p-variate responses, and an Nxk matrix X of covariates, to fit the model Y = X*B + e with estimation of the kxp matrix of coefficients B and estimation of the pxp matrix of covariances between the p variates in Y. I haven't managed to find a
2003 Feb 03
1
Lattice not plotting within loop
Something I don't understand ... (!) With the lattice library loaded, I have a loop for( Z in ... ) { ... xyplot(y~x | t, xlab=..., ylab=... ) } and no plot appears on the R graphics device. However, when I run the commands within {...} separately for each instance of Z, I get the plot displayed each time. So it looks as though "xyplot" is not outputting to the graphics
2003 Jun 18
1
Query: Sd2Rd and nroff macros in S docs
Documentation for S3 functions is apparently written in troff markup with macro tags like .BG .FN .TL .DN .CS ... Inspection of S3 documentation source files gives a pretty clear idea of what these mean, semantically (and Sd2Rd is a perl script which converts this markup into the Rd format, providing further semantic information along the way). My query is: Can anyone point to troff
2003 Jul 15
2
printf and friends in R?
Hi folks Does R have anything straightforwardly resembling the commands fprintf, sprintf, printf (derived from C, and present in octave and matlab)? As in printf(format_string, ... ) where "format_string" defines the print format (including any fixed text) and "..." is a list of variables whose values are to be inserted into the line. Example: printf("Case
2003 Jul 15
0
Multivariate regression method
Hi Folks, Thanks to several people's suggestions and clarifications, I think I have implemented a function which computes the conditional mean and covariance matrix of a subset of the dimensions of an MV-normal variable, given the values on the other dimensions (these conditioning value can be presented as a matrix, to deal with several cases at once). The code is below, for anyone who would
2003 Aug 04
0
Specifying weird models
[re-sending this one since it apparently didn't get through yesterday] Hi Folks, I'm pondering the following type of question in the context of specifying a linear model formula. Basically, it's a matter of specifying a "non-homogeneous" model. The following example (not a real case, and over-simplified, but it illustrates the point cleanly) shows what I mean. There are 3
2003 Oct 08
1
Saving workspace image
Hi folks, On quitting R with q(), is it possible to save the workspace to a directory other than the one R was started from? (I sometimes have a project "master" directory with the major R code and data in that directory, but divisions of the project having their specific stuff in sub-directories. So if I quit while running a sub-project, I'd like to save the workspace back into
2003 Oct 28
1
Loading a "sub-package"
Hi Folks, The inspiration for this query is described below, but it prompts a general question: If one wants to use only one or a few functions from a library, is there a way to load only these, without loading the library, short of going into the package source and extracting what is needed (including of course any auxiliary functions and compiled code they may depend on)? What prompted this