I am getting a seg fault and the following error when trying to use an
embedded version of R in a multithreaded application.
Error: C stack usage is too close to the limit
When I do the same thing in a single threaded app, then there is no problem.
I've tried this with R version 2.3.0, 2.3.1
Dear R-projects users
I would like like to ask if there is any way to produce a multipanel plot
with the filled.contour function. In the help information of filled
contour it is said that this function is restricted to a full page
display.
With kind regards
Prodromos Zanis
--
Cal == Cal Stats [EMAIL PROTECTED]
on Sun, 25 Jun 2006 01:14:34 -0700 (PDT) writes:
Cal Hi.. I have to invert a 15000 x 15000 matrix
Cal (generalized inverse). I do run the process on a fairly
Cal powerful computer. but still complains indufficient
Cal memory.
Cal
Prodromos Zanis wrote:
Dear R-projects users
I would like like to ask if there is any way to produce a multipanel plot
with the filled.contour function. In the help information of filled
contour it is said that this function is restricted to a full page
display.
With kind regards
Just an idea: how about using the \usage for the formal syntax, and
\synopsis for the user syntax, i.e. x/y ?
Not sure it will work, but it might be worth a try... :-)
--
Bjørn-Helge Mevik
__
R-help@stat.math.ethz.ch mailing list
On Fri, Jun 23, 2006 at 01:43:54PM -0500, Marc Schwartz (via MN) wrote:
On Fri, 2006-06-23 at 14:02 -0400, Michael H. Prager wrote:
Previous posters have argued for EPS files as a desirable transfer
format for quality reasons. This is of course true when the output is
through a
Gregory == Gregory Gentlemen [EMAIL PROTECTED]
on Sat, 24 Jun 2006 14:41:37 -0400 (EDT) writes:
Gregory Can anyone tell me the trick for obtaining the
Gregory smoother matrix from smooth.spline when there are
Gregory non-unique values for x. I have the following code
On 6/26/06, Jan T. Kim [EMAIL PROTECTED] wrote:
On Fri, Jun 23, 2006 at 01:43:54PM -0500, Marc Schwartz (via MN) wrote:
On Fri, 2006-06-23 at 14:02 -0400, Michael H. Prager wrote:
Previous posters have argued for EPS files as a desirable transfer
format for quality reasons. This is of
Dear R Users,
I am looking for a function to use within optim s.t.
- f(x,a,b,c,d,e) where I am searching over parameter values a,b,c,d,e
- a is the lower bound (typically 1)
- b is the upper bound (typically 0)
- c,d,e give the function a shift, slope and curvature between a and b
Any good
Please do study the posting guide, as there is no `R 2.4.0'.
This seems to be related to the NEWS item
o [[ on a list does not duplicate the extracted element unless
necessary. (It did not duplicate in other cases, e.g. a
pairlist.)
Now in so far as I can follow it, in
I collected eggs laid by Springtails everyday over 28 days after swich to
isotopically enriched diet. The eggs were pooled at day 7, 14, and 28 (+ day 0
= initial value) and analyzed for isotopes. After the diet switch the isotopic
values of the adults and eggs change towards those of the new
Larsen, Thomas [EMAIL PROTECTED] writes:
I collected eggs laid by Springtails everyday over 28 days after swich to
isotopically enriched diet. The eggs were pooled at day 7, 14, and 28 (+ day
0 = initial value) and analyzed for isotopes. After the diet switch the
isotopic values of the
-BEGIN PGP SIGNED MESSAGE-
Hash: SHA1
Dear R-helpers,
my data.frame is of the form
x - data.frame( f=gl(4,3), X=rep(0:2,4), p=c(.1,.2,.3))
x
f X p
1 1 0 0.1
2 1 1 0.2
3 1 2 0.3
4 2 0 0.1
5 2 1 0.2
6 2 2 0.3
7 3 0 0.1
8 3 1 0.2
9 3 2 0.3
10 4 0 0.1
11 4 1 0.2
12 4 2 0.3
Hi folks,
I usually don't do version announcements, but several people have
requested things that are in this version (2.1). Prominent among these
are that soil.texture has inflated into yet another general triangle
plot function (triax.plot), the dreaded shadow effect has materialized,
pie3D
On Fri, 2006-06-23 at 21:38 -0700, Spencer Graves wrote:
Permit me to try to repeat what I said earlier a little more clearly:
When the outcomes are constant for each subject, either all 0's or all
1's, the maximum likelihood estimate of the between-subject variance in
Inf. Any
You need to specify the row/column name as character:
y
X1 X2 X3 X4
0 0.1 0.1 0.1 0.1
1 0.2 0.2 0.2 0.2
2 0.3 0.3 0.3 0.3
y[,'X3']
[1] 0.1 0.2 0.3
y['0','X3']
[1] 0.1
On 6/26/06, Matthias Braeunig [EMAIL PROTECTED] wrote:
-BEGIN PGP SIGNED MESSAGE-
Hash: SHA1
Dear
Hi,
I am trying to convert a dataset (dataframe) into time series object using ts
function in stats package. My dataset is as follows:
df
[1] 11.08 7.08 7.08 6.08 6.08 6.08 23.08 32.08 8.08 11.08 6.08 13.08
13.83 16.83 19.83 8.83 20.83 17.83
[19] 9.83 20.83 10.83 12.83
Hi,
I recently had a problem installing the rgl package on OS X and put together
a simple patch. The patched package is available here:
http://jinome.stanford.edu/files/rgl_0.66-patched_for_gcc4.tar.gz
It can be installed with R CMD INSTALL rgl_0.66-patched_for_gcc4.tar.gz as
normal at the
I think this has already been fixed. The rgl repository had a disk
crash and is being restored today, but you can see a recent build on my
web page,
http://www.stats.uwo.ca/faculty/murdoch/software/
We are planning a new release very soon; the disk crash is a bit of a
nuisance, though.
On Mon, 26 Jun 2006, Sachin J wrote:
I am trying to convert a dataset (dataframe) into time series object
using ts function in stats package. My dataset is as follows:
df
[1] 11.08 7.08 7.08 6.08 6.08 6.08 23.08 32.08 8.08 11.08 6.08 13.08
13.83 16.83 19.83 8.83 20.83 17.83
On Mon, 26 Jun 2006, Sachin J wrote:
Hi,
I am trying to convert a dataset (dataframe) into time series object
using ts function in stats package. My dataset is as follows:
df
[1] 11.08 7.08 7.08 6.08 6.08 6.08 23.08 32.08 8.08 11.08 6.08 13.08
13.83 16.83 19.83 8.83 20.83
Folks,
The contour() function wants x and y to be in increasing order. I have
a situation where I have a grid in x and y, and associated z values,
which looks like this:
x y z
[1,] 0.00 20 1.000
[2,] 0.00 30 1.000
[3,] 0.00 40 1.000
[4,] 0.00 50
Hi Achim,
I did the following:
df - read.csv(C:/data.csv, header=TRUE,sep=,,na.strings=NA, dec=,,
strip.white=TRUE)
Note: data.csv has 10 (V1...V10) columns.
df[1]
V1
111.08
2 7.08
3 7.08
4 6.08
5 6.08
6 6.08
7
We don't have data.csv so its still not ***reproducible*** by anyone
else. To be reproducible it means that anyone can copy the code
in your post, paste it into R and get the same answer.
Suggest you post the output of
dput(df)
and then post
dput - ...the output you got from dput(df)...
Now
You are right. The df is as follows:
df[1]
V1
111.08
2 7.08
3 7.08
4 6.08
5 6.08
6 6.08
723.08
832.08
9 8.08
10 11.08
116.08
12 13.08
13 13.83
14 16.83
15 19.83
16
Sorry I meant issue dput(df) and
post
df - ...the output your got from dput(df)...
...rest of your code...
Now its reproducible.
On 6/26/06, Gabor Grothendieck [EMAIL PROTECTED] wrote:
We don't have data.csv so its still not ***reproducible*** by anyone
else. To be reproducible it means
It seems I have problem in reading the data as dataframe. It is reading it as
factors. Here is the df
df - read.csv(C:/Data.csv,header=TRUE,sep=,,na.strings=NA, dec=,,
strip.white=TRUE)
dput(df)
df - structure(list(V1 = structure(c(2, 15, 15, 14, 14, 14, 12, 13,
+ 16, 2,
On 6/25/2006 9:33 AM, Ajay Narottam Shah wrote:
Folks,
The contour() function wants x and y to be in increasing order. I have
a situation where I have a grid in x and y, and associated z values,
which looks like this:
contour() wants vectors of x and y values, and a matrix of z values,
df[] - sapply(format(df), as.numeric)
will convert it to numeric but I think the real problem is the read.csv
statement. Do commas represent separators or decimals since
you have specified comma for both? Assuming it looks like:
A,B,C
1,2,3
4,5,6
just do:
DF - read.csv(Data.csv)
str(DF)
I think it would be helpful if this were added to the contour help file.
On 6/26/06, Duncan Murdoch [EMAIL PROTECTED] wrote:
On 6/25/2006 9:33 AM, Ajay Narottam Shah wrote:
Folks,
The contour() function wants x and y to be in increasing order. I have
a situation where I have a grid in x
On 6/26/2006 10:39 AM, Gabor Grothendieck wrote:
I think it would be helpful if this were added to the contour help file.
You mean an example of building up the z matrix from points, or just a
general discussion of the issue?
Duncan Murdoch
On 6/26/06, Duncan Murdoch [EMAIL PROTECTED]
Hi Gabor,
You are correct. The real problem is with read.csv. I am not sure why? My
data looks
V1,V2,V3
11.08,21.73,13.08
7.08,37.73,6.08
7.08,11.73,21.08
I never had this problem earlier. Anyway I did
df - read.csv(Data.csv)
tsdata - ts((df),frequency = 12, start =
I think the problem is that you specified the decimal character to be
a comma so when it saw the dot it figured its not a number.
On 6/26/06, Sachin J [EMAIL PROTECTED] wrote:
Hi Gabor,
You are correct. The real problem is with read.csv. I am not sure why? My
data looks
V1,V2,V3
I think its often the case that one has 3 tuples and does not know
how to use contour with that; so, it would be nice if the contour
help page gave advice and an example and a pointer to the relevant
functions if it cannot be done by contour.
On 6/26/06, Duncan Murdoch [EMAIL PROTECTED] wrote:
Sachin J [EMAIL PROTECTED] writes:
Hi Gabor,
You are correct. The real problem is with read.csv. I am not sure why? My
data looks
V1,V2,V3
11.08,21.73,13.08
7.08,37.73,6.08
7.08,11.73,21.08
read.csv(C:/Data.csv,header=TRUE,sep=,,na.strings=NA,
dec=,,
see inline
Rick Bilonick wrote:
On Fri, 2006-06-23 at 21:38 -0700, Spencer Graves wrote:
Permit me to try to repeat what I said earlier a little more clearly:
When the outcomes are constant for each subject, either all 0's or all
1's, the maximum likelihood estimate of the
Dear coleges,
How to use the genotype from microsatelite markers
from many different populations to construct a
tridimensional phylogenetic tree with R? Any
suggestions?
Thank you very much!
Baron, Erica
Dra. Erica Baron
Universidade dos Açores - Departamento de Ciências Agrárias
Grupo de
Rick Bilonick wrote:
I guess the moral is before you do any computations you have to make
sure the procedure makes sense for the data.
Is this a candidate for the fortunes package? (an oxymoronic profound, but
obvious comment).
-- Bert Gunter
Genentech Non-Clinical Statistics
South
The R2HTML package does have a driver for sweave (see help for
RweaveHTML). This allows you to write an HTML file and add the R
commands you want (use =, @ combinations to indicate R commands and
Sexpr r-code for inline replacement), then process it with sweave and
have a final HTML file (and
Hi Dear R users,
For a pair plotting, usaully we use par( ) function. Apparently it does not
work anywhere. I want to have 3 plots in a single figure, like this:
par(mfrow=c(3,1))
densityplot( a)
densityplot(b)
densityplot(c)
But it does not work. How is it possible to have such
Hi,
Is it possible to find corresponding color code in R
for the following RGB (R185, G35 B80)?
Thanks in advance.
Best regards,
Ezhil
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting
Amir Safari wrote:
Hi Dear R users,
For a pair plotting, usaully we use par( ) function. Apparently it does not
work anywhere. I want to have 3 plots in a single figure, like this:
par(mfrow=c(3,1))
densityplot( a)
densityplot(b)
densityplot(c)
But it does not work.
I'm getting
Error in strwidth(...) : X11 font at size 15 could not be loaded
kinds of errors when I try to change cex to something *other*
than 1 (works OK otherwise).
I can't get through to RSiteSearch() right now, but as I recall
my search of it only revealed people who had problems
A Ezhil wrote:
Hi,
Is it possible to find corresponding color code in R
for the following RGB (R185, G35 B80)?
Thanks in advance.
Best regards,
Ezhil
How about:
x - c(185, 35, 80)
class(x) - hexmode
paste(#, paste(format(x), collapse = ), sep = )
[1] #b92350
I found this
Well, rgb(185, 35, 80, max=255) seems like a good place
to start.
help.search(rgb)
turns up all kinds of color-related functions.
Sarah
On 6/26/06, A Ezhil [EMAIL PROTECTED] wrote:
Hi,
Is it possible to find corresponding color code in R
for the following RGB (R185, G35 B80)?
Thanks in
?rgb
You will need to divide your values by the maximum or specify it, as in
rgb(185, 35, 80, max=255)
on 6/26/2006 12:17 PM A Ezhil said the following:
Is it possible to find corresponding color code in R
for the following RGB (R185, G35 B80)?
--
Michael Prager, Ph.D.
Southeast
You can use the following:
rgb(185, 35, 80, max=255)
Which gives #B92350
Or if you want a color name, the closest I found is maroon which is
red: 176, green: 48, blue: 96
Hope this helps,
--
Gregory (Greg) L. Snow Ph.D.
Statistical Data Center
Intermountain Healthcare
[EMAIL PROTECTED]
On 6/26/06, Sundar Dorai-Raj [EMAIL PROTECTED] wrote:
Amir Safari wrote:
Hi Dear R users,
For a pair plotting, usaully we use par( ) function. Apparently it does
not work anywhere. I want to have 3 plots in a single figure, like this:
par(mfrow=c(3,1))
densityplot( a)
A Ezhil wrote:
Hi,
Is it possible to find corresponding color code in R
for the following RGB (R185, G35 B80)?
rgb(185, 35, 80, max=255)
[1] #B92350
?rgb
Thanks in advance.
Best regards,
Ezhil
__
R-help@stat.math.ethz.ch mailing list
Hi Greg,
Thank you very much for your response. It works.
Thanks also for Michael,Sarah Sundar.
Best regards,
Ezhil
--- Greg Snow [EMAIL PROTECTED] wrote:
You can use the following:
rgb(185, 35, 80, max=255)
Which gives #B92350
Or if you want a color name, the closest I found
Kuhn, Max wrote:
Jay,
You should use RCMD install --build pkgName to create the zip file on
Windows. The zip files you see on CRAN are Windows binaries. You could
also used RCMD build pkgName, but I remember seeing a post a while
back saying that using install instead of build was best
Hi there, is there any way to compare 2 odds ratios? I
have two tests that are supposed to detect a disease
presence. So for each test, I can compute an odds
ratio. My problem is how can I compare the 2 tests by
testing whether the 2 odds ratios are the same?
Appreciate
Afternoon folks:
I'm getting a program crash when I try to run demo(rgl). The following
error details result:
RGUI caused a stack fault in module NVOPENGL.DLL at 017f:695280f0.
Registers:
EAX=0002 CS=017f EIP=695280f0 EFLGS=00010246
EBX=0001 SS=0187 ESP=00572000 EBP=004d1208
Deepayan Sarkar wrote:
On 6/26/06, Sundar Dorai-Raj [EMAIL PROTECTED] wrote:
Amir Safari wrote:
Hi Dear R users,
For a pair plotting, usaully we use par( ) function. Apparently it
does not work anywhere. I want to have 3 plots in a single figure,
like this:
Brian Lunergan wrote:
Afternoon folks:
I'm getting a program crash when I try to run demo(rgl). The following
error details result:
RGUI caused a stack fault in module NVOPENGL.DLL at 017f:695280f0.
Registers:
EAX=0002 CS=017f EIP=695280f0 EFLGS=00010246
EBX=0001 SS=0187
Hi : I think I need to use sapply but I can't figure this out.
Suppose I have two vectors : tempa ( 4, 6,10 ) and tempb
( 11,23 ,39 )
I want a function that returns 4:11,6:23 and 10:39 as vectors.
I tried :
sapply(1:length(tempa) function (z) seq(tempa[z],tempb[z])
but i got 3 really
On Mon, 2006-06-26 at 12:40 -0500, [EMAIL PROTECTED] wrote:
Hi : I think I need to use sapply but I can't figure this out.
Suppose I have two vectors : tempa ( 4, 6,10 ) and tempb
( 11,23 ,39 )
I want a function that returns 4:11,6:23 and 10:39 as vectors.
I tried :
Michael Hopkins skreiv:
The kernel density plots don¹t have to be very sophisticated i.e.
default settings and greyscale are fine we can work on details
later. What has stumped us so far is how you ‘attach’ the kernel
density results to the scatterplot results and then overlay them.
Any
On Mon, 26 Jun 2006, [EMAIL PROTECTED] wrote:
Hi : I think I need to use sapply but I can't figure this out.
Suppose I have two vectors : tempa ( 4, 6,10 ) and tempb
( 11,23 ,39 )
I want a function that returns 4:11,6:23 and 10:39 as vectors.
I tried :
sapply(1:length(tempa) function
On 6/26/06, Sundar Dorai-Raj [EMAIL PROTECTED] wrote:
Deepayan Sarkar wrote:
On 6/26/06, Sundar Dorai-Raj [EMAIL PROTECTED] wrote:
Amir Safari wrote:
Hi Dear R users,
For a pair plotting, usaully we use par( ) function. Apparently it
does not work anywhere. I want to
Hi,
How can I produce the following output in .csv format using write.table
function.
for(i in seq(1:2))
{
df - rnorm(4, mean=0, sd=1)
write.table(df,C:/output.csv, append = TRUE, quote = FALSE, sep = ,,
row.names = FALSE, col.names = TRUE)
}
Current O/p:
x
Hello all!
I hope this is the right place to post this question.
The Oregon Chapter of ASA is working with local high school teachers as one
of its outreaching program.
We hope to use and test R as teaching tools.
So, we think that a menu system (like R commander) with a few packages
Do any of the R libraries have an implementation of the Inverse Error
Function (Inverse ERF)?
ref:
http://mathworld.wolfram.com/InverseErf.html
http://functions.wolfram.com/GammaBetaErf/InverseErf/
Thanks,
Nathan
[[alternative HTML version deleted]]
Nathan Dabney wrote:
Do any of the R libraries have an implementation of the Inverse Error
Function (Inverse ERF)?
ref:
http://mathworld.wolfram.com/InverseErf.html
http://functions.wolfram.com/GammaBetaErf/InverseErf/
Thanks,
Nathan
[[alternative HTML version deleted]]
Hi,
You can use the following relation between standard normal probability
distribution (\Phi) and error function:
Erf(z) = 2 * \Phi(\sqrt(2) z) - 1
to evaluate invErf(x) in R as follows:
invErf - function(x) {
# argument x must lie between -1 and 1
qnorm((1 + x) /2) / sqrt(2)
}
For example,
On 6/26/2006 3:14 PM, Dongseok Choi wrote:
Hello all!
I hope this is the right place to post this question.
The Oregon Chapter of ASA is working with local high school teachers as one
of its outreaching program.
We hope to use and test R as teaching tools.
So, we think that a
Dear r-helpers,
I wonder if you can figure out why the following is working:
sessionInfo()
Version 2.3.1 (2006-06-01)
powerpc-apple-darwin8.6.0
attached base packages:
[1] tcltk methods stats graphics grDevices
Dear r helpers,
In my exploration of the tcltk facilities of R I've had some success
but some failures, and wonder if someone could point me to a
solution. To begin:
**
sessionInfo()
Version 2.3.1 (2006-06-01)
powerpc-apple-darwin8.6.0
When I look through archives at
https://stat.ethz.ch/pipermail/r-help/2003-October/040525.html
I see this:
Liaw, Andy wrote:
In the `Detail' section of ?princomp:
princomp only handles so-called Q-mode PCA, that is feature extraction of
variables. If a data matrix is supplied (possibly via a
Hi there, is there any way to compare 2 odds ratios? I
have two tests that are supposed to detect a disease
presence. So for each test, I can compute an odds
ratio. My problem is how can I compare the 2 tests by
testing whether the 2 odds ratios are the same?
Appreciate
Hey everyone,
I have read the paper by Ritter and Tanner(1992) on Griddy-Gibbs sampler and
I am trying to implement it in R without much luck. I was wondering if anyone
had used this or could point me to any example code.
Thanks,
Liz
Michael Kubovy [EMAIL PROTECTED] writes:
Dear r helpers,
In my exploration of the tcltk facilities of R I've had some success
but some failures, and wonder if someone could point me to a
solution. To begin:
**
sessionInfo()
Hi all,
I have a function getSomeData() that is called from command line -
getSomeData(id='1240'). The function getSomeData() calls another function
getData that needs the exact same argument passed to getSomeData(). I am
using the function call getData(paste(names(args[1]), =,
On Mon, 2006-06-26 at 10:59 -0300, Erica Baron wrote:
Dear coleges,
How to use the genotype from microsatelite markers
from many different populations to construct a
tridimensional phylogenetic tree with R? Any
suggestions?
Thank you very much!
Baron, Erica
Hi Erica,
I know nothing about
Dear R-users
I have four simple linear models, which are all in the form of a*X+b
The estimated parameter vectors are
a - c(1,2,3,4)
b - c(4,3,2,1)
My goal is to draw a plot where x-axis is X (range from -100 to 100) and
y-axis is the sum of
all four linear models
X - seq(-100,100,
Hi Spencer,
what you're describing now seems to me to be similar to the idea of
intra-class correlation. I have previously used the relative and
absolute amounts of variance represented at different hierarchical
levels of random effects, and they seemed to me to be easy to
interpret and useful,
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
On Behalf Of Taka Matzmoto
Sent: Monday, June 26, 2006 5:21 PM
To: r-help@stat.math.ethz.ch
Subject: [R] looking for a more efficient R code.
Dear R-users
I have four simple linear models, which are all in
try this:
a - c(1,2,3,4)
b - c(4,3,2,1)
X - seq(-100,100, length=1)
z - colSums(sapply(X,function(x)a*x+b))
plot(X, z, type='l')
On 6/26/06, Taka Matzmoto [EMAIL PROTECTED] wrote:
Dear R-users
I have four simple linear models, which are all in the form of a*X+b
The estimated parameter
Little algebra will convince you that your 'summing' function is just
equal to:
sum(a) + sum(b)*X
So, it is very difficult to get it any faster than:
plot(X, sum(a)+sum(b)*X)
And, by the way, most of the execution time is spent plotting, not
computing the average:
system.time({y1 -
Hello,
with 4 different linear mixed models (continuous dependent) I find that my
residuals do not follow the normality assumption (significant Shapiro-Wilk
with values equal/higher than 0.976; sample sizes 750 or 1200). I find,
instead, that my residuals are really well fitted by a t
in my previous post in which i asked about creating sequences
from two vectors of numbers, all suggestions worked.
tradevectors-mapply(seq,from=tempa,to=tempb)
or
tradevectors-sapply(1:length(tempa),function(x) seq(tempa[x],tempb[x])
both return a list with 3 components.
the problem is that
See:
http://tolstoy.newcastle.edu.au/R/help/06/06/29301.html
On 6/26/06, Aarti Dahiya [EMAIL PROTECTED] wrote:
Hi all,
I have a function getSomeData() that is called from command line -
getSomeData(id='1240'). The function getSomeData() calls another function
getData that needs the exact
I used RSearchSite and found out about using
matrix(unlist(thelist),nrow=length(thelist),byrow=TRUE)
but then i realized that I can't use this because each
component of the list can be a different length
so the matrix would have to have variable column lengths which
makes no sense. thanks for
On Mon, 2006-06-26 at 20:46 -0500, [EMAIL PROTECTED] wrote:
in my previous post in which i asked about creating sequences
from two vectors of numbers, all suggestions worked.
tradevectors-mapply(seq,from=tempa,to=tempb)
or
tradevectors-sapply(1:length(tempa),function(x)
Thanks, this is not what what I meant. I need to reshape the original
dataframe that I can access p_f[X] for numerical X. Maybe I was not
clear enough.
The problem really is that X starts at 0. Note that in my example
changing the row names to 0:2 does not have the desired effect.
jim holtman
If I have two matrices:
a-matrix(1:5,ncol=1)
[,1]
[1,]1
[2,]2
[3,]3
[4,]4
[5,]5
b-matrix(1:50,ncol=10)
[,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10]
[1,]16 11 16 21 26 31 36 4146
[2,]27 12 17 22 27 32 37 42
86 matches
Mail list logo
