search for: csilleri

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 ___

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

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

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

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