similar to: Memory getting eaten up with XML

Dear All, I have some code that works in S-Plus for writing saving a graphics file to disk :- graphsheet(type = "auto", format = "WMF", file = "G:\\north0l.wmf", pages = "auto", print.background = F, orientation="landscape", color.style="color") plot(x,y) dev.off() This works fine in S-Plus. I have tried playing

I have had the same problem and I wrote this function rmulti <- function(n, size, p) { NrDim <- length(p) if(NrDim<2) stop("The simulated variabel has to be at least 2-dimensional") res <- matrix(data=NA, nrow=n, ncol=NrDim) p <- p/sum(p) TempSize <- size for(i in 1:NrDim) { TempP <- p[i]/sum(p[i:NrDim]) TempBin <- rbinom(n=n, size=TempSize,

Dear R users I am using GLMMpql to analyse some nested negative binomial response data. How do I summarise the significance of my main effects? For example, in a standard linear mixed model (lme), I would use anova.lme to obtain an F statistic and P value for each of my main effects: how do I achieve a similar goal using GLMMpql? many thanks Sarah Richardson Sarah Richardson Plant

Dear all, I an running R on a Windows 2000 machine (1.5Gb RAM) and am trying to load the lme4 package, however, whenever I attempt to load the library "lme4" I get the following error message I have installed lme4 and Matrix locally from the zip file with no problems. R version 2.6.2 (2008-02-08) Copyright (C) 2008 The R Foundation for Statistical Computing ISBN 3-900051-07-0

Dear All, I am trying to put some Confidence intervals on some regressions from a linear model with no luck. I can extract the fitted values using 'predict', but am having difficulty in getting at the confidence intervals, or the standard errors. Any suggestions would be welcome Cheers Guy Using Version 2.1.0 (2005-04-18) on a PC vol.mod3 <-

Dear All, Using:- R : Copyright 2005, The R Foundation for Statistical Computing Version 2.1.1 (2005-06-20), ISBN 3-900051-07-0 and Jose Pinheiro, Douglas Bates, Saikat DebRoy and Deepayan Sarkar (2005). nlme: Linear and nonlinear mixed effects models. R package version 3.1-65. on a WINDOWS 2000 machine I am trying to run the scripts from the Mixed Models book and am running into some

Hello together, Is there a tool to test the statistical differences between parameter estimates of a nlsList fit? I fitted degradation data using the nlsList method and want to find out whether derived rate constants are significantly different from each other at the grouping factors soil and temperature. Regards Frank Scherr

Hi, It feels like I should be able to do something like: R CMD INSTALL lib='/usr/lib64/R/library' repos='http://proxy.url/cran' package We have a bunch of servers (compute nodes in a Rocks cluster) in an isolated subnet, there is a basic pass-through proxy set up on the firewall (the head node) which just passes HTTP requests through to our nearest CRAN mirror. when using

Given a matrix B, where B=A'A, how can I find A? In other words, if I have a matrix B which I know is another matrix A times its transpose, can I find matrix A? Thanks, Mike

Hello, From time to time Samba 3.0.2 performance suddenly fails bellow what is acceptable. One smbd process eats between 30% and 100% of CPU usage, and for the machine associated with the smbd process i got dozens of (debug level 2): [2004/02/25 17:59:18, 0] smbd/connection.c:register_message_flags(220) register_message_flags: tdb_fetch failed [2004/02/25 17:59:19, 0]

Hi Can some one please point out where i am wrong. I am trying to position set of nodes column-wise in cytoscape using RCytoscape A----D B----E C----F ------------------- g <- new ('graphNEL', edgemode='undirected') cw <- CytoscapeWindow ('smallExample', graph=RCytoscape::makeSimpleGraph()) layout (cw, 'jgraph-spring') redraw(cw) nodesFr = c('A',

In trying to install the pbatR package, I was greeted with the error Error: package 'tcltk' does not have a name space Execution halted Directly installing the package tcltk2 returned the following error: Loading required package: tcltk Error in firstlib(which.lib.loc, package) : Tcl/Tk support is not available on this system I have seen from previous posts that tcl/tk must be

Hello again, Another stumble here that is defeating me. I try: a<-readLines(url("http://feeds.feedburner.com/grokin")) t<-XML::xmlTreeParse(a, ignoreBlanks=TRUE, replaceEntities=FALSE, asText=TRUE) elem<- XML::getNodeSet(XML::xmlRoot(t),"/rss/channel/item")[[1]] And I get: Start tag expected, '<' not found Error: 1: Start tag expected, '<' not

Hi list, In a function I am writing, I need to extract the dimension names of an array. I know this can be acheived easily using dimnames() but my problem is that I want my function to be robust when the number of dimensions varies. Consider the following case: foo <- array(data = rnorm(32), dim = c(4,4,2), dimnames=list(letters[1:4], LETTERS[1:4], letters[5:6])) # What I want is to extract

Hi, I'm trying to use the XML package to read an RSS feed. To get started, I was trying to use this post as an example: http://www.r-bloggers.com/how-to-build-a-dataset-in-r-using-an-rss-feed-or-web-page/ I can replicate the beginning section of the post, but when I try to use another RSS feed I have an issue. The RSS feed I would like to use is: > URL <-

I am a relatively new user to R, and I am trying to learn more about converting data in an XML document into "2-dimensional format" such as a table or array. I might eventually wish to export this data into a relational database such as SQL, and/or to work with this data within the R package. My sample XML document is located at "

The offset argument used in glm and other functions seems to have been removed from the argument list for coxph. I am wondering if there is a reason for this and if there is a possible work-around in order to produce a cox-ph object without fitting coefficients? Thanks, Mike

I have created a data frame using the read.table command. I want to be able to access the rows by the row name, or a vector of row names. I know that you can access columns by using the data.frame.name$col.name. Is there a way to access row names in a similar manner? [[alternative HTML version deleted]]

I am trying to preprocess a large dataset of affymetrix data. Creating an affybatch is not possible with the computer I am running it on, so I have used the justRMA command to run RMA. I have read the affy document describing the justRMA command and the help documentation but I am unclear as to whether this command uses median polish after normalization. I assume this is the case but would like

I am analyzing data from an experiment with two factors: Carbon (+/-) and O3 (+/-), with 4 replicates of each treatment, and 4 harvests over a year. The treatments are assigned in a block design to individual Rings. I have approaches this as a repeated measures design. Fixed factors are Carbon, O3 and Harvest, with Ring assigned as a random variable. I have performed repeated measures analysis