similar to: variance components

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

Does anyone know how to coerce a list to an array? I have a list with equal sized matrixes and I want to coerce it to a 3D array. Why is that that I cannot feed in equal sized matrixes as layers to an array even if I set dim 3? Thanks! Katalin ___ Katalin Csillery Division of Biological Sciences University of Montana, Missoula MT 59801 Phone: 406 243 6106, E-mail: csillery at selway.umt.edu

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

Any Rails Developers in Montana? Especially Helena or Bozeman? I don''t expect there''s enough of us to form a user group, but it might be nice to talk shop ''offline'' every once and a while. Eric (in Helena) --~--~---------~--~----~------------~-------~--~----~ You received this message because you are subscribed to the Google Groups "Ruby on Rails:

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

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, 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:

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

I want all the file-names to be stored on a samba share from my XP clients to be automatically converted to 8.3 syntax, so that i do not have any further problem while transporting those files over to other networks/servers like netware. As for now i am using these below mentioned lines to support file-name mangling over my samba. --------- #I do not want my filenames to be case sensitive,

Newb here; I feel like im cheating when i actually write to the list, but i''ve been stuck on this for a while now and cant seem to find the answer in tutorials... heres the simplified version of the problem controller says: def show @exhibition = Exhibition.find(params[:id]) end content says: <% for exhibition in @exhibitions %> <%= link_to

I am trying to use varcomp to obtain the variance partitioning across different nested levels of random effects (say x,y and z). I get the three variance components (for each of my along with an additional one called 'within' from varcomp output. I am using the 'scale total variance to 1' option and though the within component is small, it does form a part of what explains the

Hello, I am trying to use varcomp to decompose the variance across multiple nested levels on a lme object. I am able to successfully do this and when I view the varcomp object I can see the individual values / estimates for the variance at different levels. However, I want to be able to extract each of them separately, as I need to build a confidence interval using bootstrapping on the sample