search for: jsgaard

...hing close to that..). The dataset is so large, that exporting it as an Excel file from SAS is not feasible (more than 65000 lines). I am reluctant to ask her to go through all the data base steps (then she'll just stick to SAS...). Can anyone help me out on that one? Thanks in advance Søren Højsgaard [[alternative HTML version deleted]]

...:stats> does not really enlighten me. How can I get to see the implementation (I guess it should be possible according to the general philosophy of the R project)? Thanks in advance S??ren ============================================================================================= S??ren H??jsgaard, PhD, Head of Research Unit Phone: +45 8999 1703 Biometry Research Unit, Fax: +45 8999 1300 Danish Institute of Agricultural Sciences E-mail: sorenh at agrsci.dk Research Centre Foulum, DK-8830 Tjele, Denmark Homepage : http://www.jbs.agrsci.dk/~sorenh/

Thanks. I wonder if there is a general way of extracting var(u) and var(e), which would be needed to simulate u and e. Clearly, one can get the estimated parameters, but is there a clever way of 'setting up' the matrices?? Best S??ren On 6/19/05, S??ren H??jsgaard <Soren.Hojsgaard at agrsci.dk> wrote: > I would like to draw a sample from a linear mixed model y=Xb+Zu+e which has been fitted with lme(), i.e. a model y ~ N(Xb, C), where C=Z cov(u) Z' + cov(e). > I've tried to figure out how to extract C from an lme object, because that would...

...people around the World, but perhaps some other people from the R-community do and would be able to/interested in establishing such a connection... Just and idea; comments appreciated :-) S?ren ============================================================================================= S?ren H?jsgaard, MSc, PhD (Statistics), Phone: +45 8999 1703 Head of Statistics and Decision Theory Research Group, Fax: +45 8999 1300 Danish Institute of Agricultural Sciences sorenh at agrsci.dk Research Centre Foulum, DK-8830 Tjele, Denmark http://genetics.agrsci.dk/...

# The result I am after is the result after a substitution in an expression, such as substitute(expression(a+b+c), list(a=1)) expression(1 + b + c) # However, the way I want to do it is for a an expression "stored as a variable" as (expr <- expression(a+b+c)) expression(a + b + c) # a) The following does not work (expr2 <- substitute(expr, list(a=1))) expr # b) - whereas this

Dear all, One of my students have installed R1.9.0 on windows, and gets the fatal error 'invalid HOMEDRIVE' Can anyone help her/me out on that one? Thanks in advance Søren Højsgaard [[alternative HTML version deleted]]

...like to have both the current version of R and R-devel installed. I've searched the net for info on doing this but so far without luck. Can anyone help me getting this up and running? A related question: How do I upgrade the current version of R to the patched versions? Best regards S?ren H?jsgaard

Hi, SAS transport files created with the xport engine in SAS can be read using read.xport. However, the xport engine only works with SAS version 6, and consequently long variable names are not allowed... Can anyone tell me how to get SAS data (ver 8) into R (easily)? Thanks in advance S?ren H?jsgaard sorenh at agrsci.dk http://www.jbs.agrsci.dk/~sorenh ------------------------------------------------------------------------- Time flies like an arrow, fruit flies like a banana. ------------------------------------------------------------------------- -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-...

...call.ole.method release.ole.object for communicating with other programs which work as a COM server (on Windows). Is it possible to do something similar in R (I've studied the 'connections' facilities, but they do not seem to work). ========================================== S?ren H?jsgaard, PhD, Senior Scientist Biometry Research Unit Danish Institute of Agricultural Sciences Research Centre Foulum, DK-8830 Tjele, Denmark Phone: +45 8999 1703 E-mail : sorenh at agrsci.dk Homepage : http://www.jbs.agrsci.dk/~sorenh/

Dear all, In SAS (GLM and MIXED) estimable functions (linear functions of the parameters) can be specified in the ESTIMATE and CONTRAST statements. Has anyone written a similar "utility" for use in connection with lm? Thanks in advance S?ren H?jsgaard -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- 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...

...ta from the package the first four columns are read in as numeric variables. This is consistent with the documentation of read.table - but it is not what I want! I can of course change the coding of the variables, but there ought to be another way. Can anyone help me on that? Best regards S??ren H??jsgaard

Dear all, Starting from a recent version of R, the $ became "unusable" on atomic vectors, e.g. > x <- c(a=1,b=2) > x$a NULL Warning message: In x$a : $ operator is invalid for atomic vectors, returning NULL I can of course do > x['a'] - but that requires more typing (5 characters rather than 2). Apologies if I've missed a an announcement regarding this, but

I've defined the function getFunNames <- function(FUN){ if (!is.list(FUN)) fun.names <- paste(deparse(substitute(FUN)), collapse = " ") else fun.names <- unlist(lapply(substitute(FUN)[-1], function(a) paste(a))) fun.names } which gives what I want : > getFunNames(mean) [1] "mean" > getFunNames(ff) [1] "ff" >

Dear all; I would like to draw a graph with vertices and edges, and I guess tcl/tk would be appropriate. I don't know tcl/tk but have googled for a 10-page (or so) introduction to 'getting started with tcl/tk in R' but without any luck. - Does anyone know of the existence of such a document or any other web-based material on the subject? - Does anyone have an (informative) piece of

...really huge data set (the kind that loglm() was designed to handle, after all) do you really want to be carting around several copies of it around in memory? That said, do you have a suggested patch that might accommodate this refinement? Bill Venables. -----Original Message----- From: S?ren H?jsgaard [mailto:Soren.Hojsgaard at agrsci.dk] Sent: Thursday, October 30, 2003 8:10 AM To: r-help at stat.math.ethz.ch Subject: [R] loglm() uses only a reference to data, and not data itself - is that on purpose?? Dear all, Consider the following: library(MASS); data(HairEyeColor) l1 <- loglm(~ 1+2+...

I would like to see how the matrix multiplication operator %*% is implemented (because I want to see which external Fortran/C routines are used). How can I do so? Best S??ren

In xyplot, I would like to get the "data into the panel function" in the following sense: Consider xyplot(scc~time|cowidp, data=cow.s,type=c("l"), panel=function(x,y,subscripts,...){ panel.xyplot(x,y,...) vvv<-cow.s[which(!is.na(cow.s[subscripts,"mastreat"])),"time"] panel.abline(v=vvv,col="red",lwd=2) } ) If I want to use a

Hi, I would like to create a plot of y1,y2,y3,y4 against x for several subjects such that y1 and y2 are plotted against x in one panel and y3 and y4 against x in another panel. Thus if there are 3 subjects I should end up with 6 panels. Is there a simple way of doing so (i.e. without calling xyplot() several times, and then padding the results together)?? Regards S?ren

I would like to draw a sample from a linear mixed model y=Xb+Zu+e which has been fitted with lme(), i.e. a model y ~ N(Xb, C), where C=Z cov(u) Z' + cov(e). I've tried to figure out how to extract C from an lme object, because that would solve my problem when also using the predict() function, but without any luck. Can anyone help on that?

...rce files installing indices Error in load(zfile, envir = envir) : input has been corrupted, with LF replaced by CR Execution halted make[2]: *** [indices] Error 1 make[1]: *** [all] Error 2 ============================================================================================= S?ren H?jsgaard, MSc, PhD (Statistics), Phone: +45 8999 1703 Head of Statistics and Decision Theory Research Group, Fax: +45 8999 1300 Danish Institute of Agricultural Sciences sorenh at agrsci.dk Research Centre Foulum, DK-8830 Tjele, Denmark http://genetics.agrsci.dk/...