similar to: apply with two matrixes

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

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 ___

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

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

Joyful. I'm adapting a FORTRAN 77 package for use with R. Pretty straightforward. Except for a glitch it took me some time to figure out. This existing package has subroutines which have parameters called "NA". So, I called subroutines like bnodes <- function(n, lst, lptr, lend, nodes, nb, na, nt) { ensure.all.numeric(list(n, lst, lptr, lend, nodes, nb, n.a, nt),

Hello everyone, Could you please help me find out if R supports matrixes inside matrixes? This is what I would like to do I have an area map of humidity per km. I would like at every cell to keep also information about the height of this area, the current temperature etc. Is something like that supported? I would like to thank you in advance for your help Best Regards Alex

Hello everyone. For reasons too long to explain I wanted to do plots similar to histograms with plot(type="h"). I ran into a problem - if I set line width too high, histogram isn't accurate anymore. For example: par(lend=2) plot(c(2,4,3,2),ylim=c(0,5), type="h") abline(h=3) Column 3 appears just as high as it should. But if I do par(lend=2)

Hi, I am trying to plot odds ratios and the corresponding confidence intervals in horizontal segments. It would be ideal if the confidence interval segment can be drawn with little vertical bars at both ends. I have tried very hard to change the type of ends by using 'lend' arguments, but cannot make it. I even tried 'arrows()', but still failed. Following is the code I use:

Dear R Users, close to the end of this I used wireframe to create a 3D plot from a matrix. The x and y axis tick labels (1-6) for each were created from the matrix being a 6X6 matrix. I need the axis tick labels to be the row and column headings which you can see in the output (mat.x). I have tried several work arounds, but they have not been successful. Thanks in advance. keith rm(list=ls())

Full_Name: Christophe Genolini Version: 2.8.1, but also 2.9 OS: Windows XP Submission from: (NULL) (82.225.59.146) I am using matplot with the option lend="butt", but only the first line (the black) is printed correctly : > matplot(matrix(1:9,3),type="c",lwd=10,lty=1,lend="butt") Gabor Grothendieck find the problem in matplot code: the ... is passed to plot

Dear all, I would like to multiply two matrixes with the different dimension column by column. Let make an example: If I have two matrixes "X" and "Y"as follow: X<- matrix(1:12, nrow=4, ncol=3, dimnames=list(c("A","B","C","D"), c("stage1","stage2","stage3"))) Y<- matrix(1:28, nrow=4, ncol=7,

Hello! I've a problem with the combn function and a set of vector. I would like to make a simple combination where, instead of scalars, i would like to combine vector, in order to form matrixes. In other words, i have nineteen 6-items vectors (for example coef1-coef19), that i would like to combine in n!/k!(n-k)! 6x6 matrixes. I tried with a code like this mma <-

The Sjobeck Company provides Asterisk Integration, Configuration, Support, and Training. We are a crack team of Unix, Windows and Apple system integrators with 10+ years of experience working with clients both large and small. The Sjobeck Company can provide turn-key solutions, or design, build, deliver, install, configure and deploy solutions any where in the world. We also do performance

Dear all I have two 10*10 matrixes and I would like to compare theirs contents. By the word content I mean to check visually (not with any mathematical formulation) how similar are the contents. I also know edit that prints my matrix in the scree but still one edit blocks the prompt to launch a second edit() screen. What is the best way to compare these two matrices? I would like to thank you