similar to: predict function

Hello. A quick question about incorporating variation due to study in the metafor package. I'm working with a particular data set for meta-analysis where some studies have multiple measurements. Others do not. So, let's say the effect I'm looking at is response to two different kinds of drug treatment - let's call their effect sizes T1 and T2. Some studies have multiple

I have got following error in drawing wavelet fitting. can some one help? > library(faraway) > data(lidar) > newlidar<-lidar[c(1:128),] > library(wavethresh) > wds <- wd(newlidar$logratio) > draw(wds) Error in plot.default(x = x, y = zwr, main = main, sub = sub, xlab = xlab, : formal argument "type" matched by multiple actual arguments [[alternative HTML

HI everyone, I''m trying to assign colors to multiple lines in a graph. Problem is I don''t want to type in as many colors as there are lines....is there a way around this? In brief, I''m plotting the logratio for up to 60 samples and want a different color for each sample. Here is the code I''m using now.. Any help is greatly appreciated.. Best LT data <-

I am (desperately) trying to get used to using a Mac here at my new location. (Why *anyone* would ever use anything other than Linux, except under duress as I am, totally escapes me, but that's another story.) Fortunately much of the Mac OSX is actually Unix, so a civilized person can manage to carry on ... But there are some things. (Like this <expletive deleted> mailer ... But

Dear R-users, I have an outcome variable and I'm unsure about how to treat it. Any advice? I have spending data for each county in the state of California (N=58). Each county has been allocated money to spend on any one of the following four categories: A, B, C, and D. Each county may spend the money in any way they see fit. This also means that the county need not spend all the money that

Hi list I am just beginning to understand the complexities of linear mixed effects models. Maybe someone can give advise concerning the following problem: I have two groups of surgical patients in which repeated laboratory measurements were taken over time after surgery. I decided that lme would be the best model to fit the data. I already fitted the model lme(logratio ~ gr*I(pod-10) +

An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20091126/6366832d/attachment-0001.pl>

I am attempting to created colClasses for several tables, then read only specific columns. There are two different table layouts that I am working with. If I use exclusively one layout, the script works perfectly, but if I mix the layouts, it fails to grab the correct columns form layout that is read in second. It appears that colClasses fails to adopt the new structure after the first iteration.

Hi, I am new to R. I am trying to figure out how to graph expression data from the TCGA database. If I understand correctly the expression data I have downloaded is from a microarray using the AgilentG4502A. I've had trouble reading into R in the level I, level II, and the gene expression analysis data using >dat<-read.table("C:\\file.txt", header=T, row.names=1) for

Dear R-Help members, I am using Simon Wood`s mgcv package version1.7-22and R version 3.0.0 (2013-04-03) for fitting a GAM-Model to the LIDAR Data contained in the "SemiPar" package. Here is the code for fitting the model and for plotting the result: data("lidar") attach(lidar) ### # mgcv fitting ### gam_fit <- gam(logratio ~ s(range, k = 40, bs = "cr"), gamma

Is there a way of calling the method of the root class directly? Here is an example explaing my question. Everything relates to S3/UseMethod classes and methods. I have an object x of class MicroarrayMatrix, which inherits from class Matrix. The "core" object of x is a simple matrix. In other words, I basically have x <- matrix(logratios, ncol=nbrOfSpots, nrow=nbrOfSlides)

I'd just like to put in my two bob's worth: The ***only*** thing for which I prefer SAS to R/Splus, but it's a very important thing, is the analysis of mixed linear models. One can do the same analyses in both languages. The difference is that the SAS syntax for specifying such a model, and its output, is comprehensible to the human mind, which is what I'm equipped with. The

Hello to the list! Someone said I could use syslinux to get GRUB to boot SCSI and/or IDE. My problem is that I have Windows on IDE (XP) and Linux on SCSI (Ubuntu). When GRUB is on the SCSI drive it can't see the IDE drive, possibly because the BIOs doesn't support that. So I have to go into the BIOs and change the boot order on every O/S change. Will I just have the same problem as

Saw this command on slashdot (http://ask.slashdot.org/article.pl?sid=04/03/24/0434224&mode=flat&tid=126&tid=95 ) to basically make a box behave like a managed switch: "tc qdisc add dev eth1 root tbf rate 250kbit latency 20ms burst 2kb" Followup:

Hy all, I have written a small library that provide a new distribution function estimator based on IFS (that are essentially fractals). I would be pleased if any of you can build hte library for Unix-type machines and Windows implementations with the included makefiles as a Mac library is already working. As I'm not able to test these platforms any modification to makefiles is welcome.

While playing around with panther I discovered (with surprise) that tcltk seems to work (even if not smoothly) with RAqua without first calling tkStartGUI. I switched back to 10.2.6 and it works as well. I think this is due to the "last minute" fix in RAqua now using idle timers because of too much cpu usage. Some "mouse" trick is needed though. If you want to test it you

After I upgraded my installation of R for OS X (Tiger) to 1.12, it consistently crashes when calling the function acf() [stats lib]. Roberto Osorio