similar to: lists to arrays

Dear All, I need to calculate the variance components in a mixed effect model (one fixed and one random effect) with REML (maximizing the proportion of the likelihood that does not depend on the fixed effects). In S+ there is the varcomp function, but I would like to do it in R. Is there a way to do that? Thanks! Katalin ___ Katalin Csillery Division of Biological Sciences University of

Hi, does anyone know how to make pairs of boxplots? E.g. 5 pairs of boxes. I figured a not very elegant way, by inserting a dummy NA field to the list where the data is given. So, where I have the dummy I got more distance between boxes, as a result it looks like groups. But I am having problems with putting the names under the groups, and the thick marks. Any better idea? Thanks, Kati ___

Hi, I have a coding problem. I have data pairs, and I want to calculate a parameter for the data pairs as a whole data set (so not using the information that thely are pairs) I am looking for bias in that certain parameter so I want to bootstrap these data pairs(!), and than calculate the parameter. I run into a coding problem. How should I define theta in order to resample the pairs, but than

Dear All, I have a function which takes a n x m matrix as an argument and returns an n x n matrix. I want to take bootstrap samples form the input matrix in the way as each row represent a multivariate observation, so each bootstrap sample would be an n x m matrix, and on each sample I want to calculate the n x n matrix. This task can be done with the sample function, but I would like to use

Dear All, I am trying to use the bca.ci function on a boot.out object which consists a few NaN's and I want to ignore those NaN's, and get a ci only for the "normal" values. boot.out$t has R number of values for 3000 different statistics, so when I use boot.ci(boot.out, index=i) and i happens to be a column in boot.out$t with some NaN's in there I get an error message. I

Hi, I have a function which does a certain task with two vectors, f1 <- function(a,b){body} I also have a list of matrixes (all with the same dim's), and for each column of each matrix in the list I want to use "f1", in such way that it gives the vector "a" in the first argument of "f1". The second argument of the function "b" also comes form a

Hello, I want to bootstarp a sample, and I need table() for each sample taken, like if R=100, I need all 100 table(). It seems to mee that boot() works anly if you assing statistics, that gives you back one number, like mean. Although the ?boot/statistics says that it can be a vector containing statistics. table() gives an array, so it might work. Does anybody have an idea how? I konw it can be

Hi, I've been using R 1.7.1 but now I need to upgrade to R-1.9.1. Since it is not on the Fink tree yet I tried to install it from source to my sw library. The configuration runs fine, but make doesn't. Here is the error message I get: /sw/lib/libg2c.a(err.o) definition of common f(short, void, int, char) (size 4) /usr/bin/libtool: internal link edit command failed make[4]: ***

Beyond a general problem with X11, that it sometimes doesn't display -especially when it is inactive for a while- new windows (like new terminal, R device, or emacs) and only the restart helps, I found a problem which might also be especially the X11's problem. >From an apple Terminal window, the following code (which is the first boot.ci example) works just fine.

Dear All, I have problem installing the 'kinship' package. I used the 'install.packages' from R command line and had no problem installing other packages before. Here is what I get: * Installing *source* package 'kinship' ... ** libs gcc -no-cpp-precomp -I/Library/Frameworks/R.framework/Resources/include -I/usr/local/include -fno-common -g -O2 -c agfit6b.c -o agfit6b.o

It''s a new month, and that means another meeting of the Missoula Ruby and Ruby on Rails Users Group. As usual, we''ll gather at a table in the back half of Sean Kelly''s, tomorrow June 5th, at 6:00pm. http://groups.google.com/group/ms-rug --~--~---------~--~----~------------~-------~--~----~ You received this message because you are subscribed to the Google Groups

The first Ruby Tuesday gathering for 2007 of the Missoula Ruby & Ruby on Rails Group will be January 2nd. (Note that''s the evening of the *first* workday after the long weekend. It''ll feel very much like Monday, but it''s still Tuesday.) Nothing very formal or organized. Just a friendly gathering of enthusiasts to chat about Ruby & RoR and related topics over

Hello All, I am looking for an R library/function that allows the specification of a custom link function in a linear mixed model. I've been using glm() in library MASS to fit fixed-effect models with a custom link but my study design demands mixed models. Any suggestions on the best R library/function to achieve this would be greatly appreciated. I have tried, to no avail, to

Hello, I am interested in fitting a generalized nonlinear regression (gnlr) model with negative binomial errors. I have found Jim Lindsay's package that will do gnlr, but I have having trouble with the particular model I am interested in fitting. It is a threshhold model, where below a certain value of one of the parameters being fitted, the model changes. Here is a sample: Cones:

Dear All, I am using xyplot() with many groups like this: statselect <- levels(dat$stat) xyplot(relmse~T|lambda, groups=stat, data=dat, panel = panel.superpose, key=simpleKey(statselect, lines=T)) However, I want lines not scatterplots and if I set panel.groups=panel.linejoin that connects the lines according to the relmse and not according the consecutive values of T. Is there

[Arrow] [Idea] [Exclamation] Subject: ~RED ALERT**URGENT 2009 [CDC::"girl's suicide rate 'skyrockets' 76%" [Shocked] .....What is causing this [Question] and, Help stop this [Exclamation] ]~ ******* ------>>> Which is more important, $$$ or your daughters [Question] <<<------ ******* ~~~RED ALERT~~~ {~See CDC quote---76% 'skyrocket' rise is real

Hi again, I was contemplating the solution using base 3: set.seed(3) mat2 <- matrix(sample(0:2, 15, replace=T), 5, 3) Extracting the line numbers is simple: bases <- c(3, 3, 3)^(2:0) # or just 3^(2:0) colSums(apply(mat2, 1, function(x) x*bases)) + 1 [1] 7 23 25 8 1 The problem is sometimes the columns have different number of levels, as in: mat1 <- expand.grid(0:2, 0:2,

Hello, I am trying to import a list of "structure" objects into R (data for BUGS, actually), but am having little success. In the help file for scan there dies not seem to be a list of valid values for the "what" keyword. Can anyone provide me with one? I have pasted in a sample input file that I am trying to import below this message. I have tried

Hello all, I have a .csv file like below. Tool,Step_Number,Data1,Data2... etc up to 100 columns. A,1,0,1 A,2,3,1 A,3,2,1 . . B,1,3,2 B,2,1,2 B,3,3,2 . . ...... so on upto 50 rows where the column "*Tool*" has distinct steps in second column "*Step_Number*",but both have same entries in Step_Number column. I want the output like below.

Hi I have the following array: test <- array(c(1:16), dim = c(3,4,3)) test ## I call some enries using an index array test.ind <- array(rbind(c(1,2,1), c(3,3,2)), dim = c(2,3)) test[test.ind] ## suppose I want all values in the 2nd row and 4th col over ## all three 3rd dimensions test[2,4,] how to specify a test.ind array with the last index left with ',' i.e test.ind should be