Displaying 20 results from an estimated 300 matches similar to: "Can boot return matrix?"
2003 May 23
1
variance components
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
2003 Feb 19
1
groups of boxplots
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
___
2003 Apr 24
1
bca ci's and NaN's in boot.out
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
2002 Oct 15
1
No subject
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
2003 Jan 22
1
lists to arrays
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
2003 Mar 12
1
apply with two matrixes
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
2003 Apr 04
1
Problems with Mac OS X Beta 3 display?
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.
2002 Nov 07
1
boot&table
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
2004 Oct 06
1
R-1.9.1 on OS 10.2
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]: ***
2001 May 19
2
calculations on diagonals of a matrix
Given an nxm matrix A I want to compute the nxm matrix B whose ij-th
element is the sum of the elements of A lying on the diagonal that ends
with element ij, i.e.,
b_ij = a_ij + a_(i-1)(j-1) + a_(i-2)(j-2) + ...
In APL (which I no longer use), I would use the 'rotate' operator to derive
an array whose columns are diagonals of the given array and then cumulate
down columns. Is
2006 Apr 01
1
Using vectorization instead of for loop for performing a calculation efficiently
I am trying to write an efficient function that will do the following:
Given an nxm matrix, 10 rows (observations) by 10 columns (samples)
for each row, test of all values in the row are greater than a value k
If all values are greater than k, then set all values to NA (or something),
Return an nxm matrix with the modified rows.
If I do this with a matrix of 20,000 rows, I will be waiting until
2008 Sep 17
1
Exact test in nxm contingency table
Hello,
I am trying to find a permutation test that works on a general nxm table. The data set is small enough to have cells with too small counts to make chi2-approximation invalid. If the table was a 2x2 contingency table I would like to use a Fsher exact test (fisher.test) but that wont work in this general table.
Does there exist a general function for this test.
Best regards,
Magnus
2003 Nov 03
3
A matrix is full rank is equal to having independent columns?
Dear R listers,
Just a simple question.
If we say an nxm matrix (n>m) is full rank of m,
does this mean that this matrix has linearly independent columns?
They are the same definition or needs some proof?
Thanks for your answer.
Fred
[[alternative HTML version deleted]]
2018 May 04
0
RFC: virtual-like methods via LLVM-style RTTI
On 3 May 2018, at 22:09, David Zarzycki via llvm-dev <llvm-dev at lists.llvm.org> wrote:
>
> Hello,
>
> In an effort to help LLVM-style projects save memory, I’ve been toying with some macros that provide an alternative to C++ vtables that use LLVM-style RTTI design patterns instead. Is this something that LLVM or sub-projects think is worth pursuing? Or are the macros below
2007 Jul 31
3
Nonlinear optimization with constraints
Hello R community,
I am using R for creating a model using optimization. I would like to ask if there is R-function/package for solving the problem below:
Minimize sum(abs(exp^(Ai1 x1 + Ai2 x2 + ... + Aim xm - bi) - 1)), for each i = 1, ..., n.
subject to Ai1 x1 + Ai2 x2 + ... + Ajm xm - bi <= c, where c is a scalar.
(x is a vector of variables, A is nxm matrix, b is a vector)
2007 Mar 05
2
Linear programming with sparse matrix input format?
Hi.
I am aware of three different R packages for linear programming: glpk,
linprog, lpSolve. From what I can tell, if there are N variables and M
constraints, all these solvers require the full NxM constraint matrix. Some
linear solvers I know of (not in R) have a sparse matrix input format. Are
there any linear solvers in R that have a sparse matrix input format?
(including the
2011 Jul 19
3
calculating the mean of a random matrix (by row) and some general questions
Hi everyone!
I'm trying to teach myself R in order to do some data analysis. I'm a
mathematics student and (only) familiar with matlab and latex. I'm working
trough the "official" introduction to R at the moment, while simultaneously
solving some exercises I found in the web. Before I post my (probably
stupid) question, I'd like to ask you for some general advice. How do
2004 Oct 18
1
installing package 'kinship'
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
2006 Sep 21
1
How to generating diagnal blocks ?
Hi,
I am trying to creat a matrix with diagnal blocks, say, I have a matrix X of any dimension (nxm) ,and would like to have:
X
X
X
2009 Mar 13
1
Help with Function!
Dear All,
I need to write 'n' functions on 'm' variables. The functions should be
constructed according to the values of an (nxm) matrix of '1/0' values as
follows. For example,
if row1 is equal to ,say [1 0 .......0 0] then f1 <- (1+x1)
if row 2 is equal to, say [1 1 1 0...0 1] then f2
<-(1+x1)*(1+x2)*(1+x3)*(1+xm)
if row n is equal to [0 1 0 1 1 0 ..... 0] then