On Mon, 29 Aug 2005, Ferran Carrascosa wrote:
Hi,
I have a matrix with 700.000 x 10.000 cells with floating point data.
I would like to work with the entire table but I have a lot of memory
problems. I have read the ?memory
I work with Win 2000 with R2.1.0
The only solution that I have
Thanks Prof Brian for your answers,
I have read about 'ref' package to work with more efficient memory
work. Anybody know if this package could help me to work with a
700.000 x 10.000 matrix?
I will have problems with ref package on:
- Limit of 2 Gb in R for Windows.
-The maximum cells in one
Re: [R] memoryFrankly:
you probably CAN NOT unles you have a very, very powerful machine.
5.6e10 bytes = (if I got it right) 35 Gb of memory.
Do you have that much?
If you have sparse matrices (those with many, many zeros), perhaps you can work
with sparse matrices methods... but if you had them
Ferran,
What are you trying to do with such a large matrix? with 7e9 cells
and a linear algorithm which is quite unlikely, your problem solution
is likely to take a very long time(tm)... just quickly... at one
micro-second per operation (very optimistic?) and 7e9 operations,
thats
7e9/1e6/60
Sorry, your message was not sent out to 'sap-r3-apo'
because it appears not to have a Subject.
If you believe you are receiving this message in error, resend
your message with a Subject in it.
For your information, the message is quoted below:
Return-Path: r-help@stat.math.ethz.ch
Received:
b wrote:
Hi
Tried to install version 2011 of R onto my XP
machine. It installs
fine, and I verified the checksums, but I click and
nothing happens?
i.e the icon pointing to RGUI.exe does nothing??
any clues??
What happens if you run it from a command window? Change to the
Hi,
here is some info about the first part of my homework, for those, who want
to break down their signal (heart beat or whatever) into a collection of pure
sin waves to analyse main frequency magnitudes and phases.
First some very un-mathematical applied theory:
If you sample a waveform
On Mon, 29 Aug 2005, Peter Dalgaard wrote:
Torsten Hothorn [EMAIL PROTECTED] writes:
Dear All,
is there a stratified version of the Wilcoxon test (also known as van
Elteren test) available in R?
you can plug it together using the `coin' infrastructure (see the
examples in
Hello,
I use the function qda (package MASS) to obtain a qda object like below.
x.qda = qda(x, group)
the group is a factor of two levels
and use this object to do the prediction below.
y.pred = predict(x.qda, y)
after that, I set different prediction priors like below, but the
results of
Hi
library(magic)
?minmax
[
the basic idea is min(x) == max(x)
]
best wishes
Robin
On 29 Aug 2005, at 20:35, Vincent Goulet wrote:
Is there a canonical way to check if all elements of a vector or
matrix are
the same? Solutions below work, but look hackish to me.
x - rep(1, 10)
Hi R-masters!
I have a problem and need your help.
I have 9 discrete variables with 2 levels each.
In exploratory analisys I generate one matrix with chi-square for tables
with 2 ariables each with this script
setwd(F:/)
dados-read.csv(log.csv)[,2:10]
dados.x-matrix(NA,ncol=9,nrow=9)
for(i in
Dear list,
I cannot make the latex command to output a ftable objet the way I
want it. Is it posible?
I found a post in the archives saying that one should use the rgroup
and n.rgroup arguments to supply the row names, but so far I have been
unsuccessful.
This is what I have:
Dear R-list readers:
Royston described the effect of sample size on the p-value obtained from
the Shapiro-Francia test (Estimating departure from normality. Stat Med
1991;10:1283-93). He developed two indices from the Shapiro-Francia
test (i) V' - an index of departure from normality
I have own function wrote. I used an algorithm, which was written in Matlab.
in matlab:
...
gamma = inv(v)*g;
...
#v = matrix of variable size, v=vv(k) = k=2 = dimension of v 2x2
#g = a line vector with 4 elements e.g. g=[1,0,2,0];
my rewritten r-file:
...
gamma = solve(v)*g;
...
which is my
Torsten Hothorn [EMAIL PROTECTED] writes:
On Mon, 29 Aug 2005, Peter Dalgaard wrote:
Torsten Hothorn [EMAIL PROTECTED] writes:
Dear All,
is there a stratified version of the Wilcoxon test (also known as van
Elteren test) available in R?
you can plug it together using
Weiwei:
Job searches are difficult! One obvious answer is Amstat news and the ASA
job site, but there may be many not posted in these places. Most large
Pharmas (Pfizer, GSK, Merck, etc.) have (relatively small) pre-/non-
clinical research groups, so you might check on their websites for open
Hi Fredrik,
What you need to do is to massage your table a little to get it into the
appropriate structure and then use the rgroup and n.rgroup options.
Here's an example using the Titanic data that come with R. Note that
this is not necessarily (or even remotely likely?) the best way to get
your
Christian Hinz wrote:
I have own function wrote. I used an algorithm, which was written in Matlab.
in matlab:
...
gamma = inv(v)*g;
...
#v = matrix of variable size, v=vv(k) = k=2 = dimension of v 2x2
#g = a line vector with 4 elements e.g. g=[1,0,2,0];
my rewritten r-file:
...
Please provide an example that doesn't work, I made one up quickly just to
see what would happend and it worked without errors.
g=c(1,0,2,0)
g
[1] 1 0 2 0
v-c(5,10,15,20)
dim(v) - c(2,2)
v
[,1] [,2]
[1,] 5 15
[2,] 10 20
solve(v)*g
[,1] [,2]
[1,] -0.4 0.6
[2,] 0.0 0.0
On 8/30/05,
On Tue, 30 Aug 2005, Peter Dalgaard wrote:
Torsten Hothorn [EMAIL PROTECTED] writes:
On Mon, 29 Aug 2005, Peter Dalgaard wrote:
Torsten Hothorn [EMAIL PROTECTED] writes:
Dear All,
is there a stratified version of the Wilcoxon test (also known as van
Elteren test)
I have two matrices (see example below) and I want the differences for the
matching row numbers, but the row numbers are not identical in the two
matrices. There are probably many ways to do this. Anyone know of any easy
way to do this? I could loop over them, but you know what they say about
Hi again Fredrik,
Here's a slightly better version (sex is no longer the first column, it
is used by the rowname option in latex instead).
library(Hmisc)
x - ftable(Titanic, row.vars = 1:2)
x
Age Child Adult
SurvivedNo YesNo Yes
Class Sex
1st Male
roger bos wrote:
I have two matrices (see example below) and I want the differences for the
matching row numbers, but the row numbers are not identical in the two
matrices. There are probably many ways to do this. Anyone know of any easy
way to do this? I could loop over them, but you know
Uwe,
Thanks. I have used merge a lot, but I didn't realize you could merge on
row names. I got it working now:
test - merge(out.r[1:5, 1:3], out.p[1:5, 1:3], by=row.names, all.x=TRUE)
Thanks,
Roger
On 8/30/05, Uwe Ligges [EMAIL PROTECTED] wrote:
roger bos wrote:
I have two matrices
Thank you so much David.
That example provided me with enough information for me to convert the
table to latex, with
subdivisions intact.
/Fredrik
On 8/30/05, David Whiting [EMAIL PROTECTED] wrote:
Hi Fredrik,
What you need to do is to massage your table a little to get it into the
[EMAIL PROTECTED] writes:
Nice to see that the old code made sense.
Nice to see that the new code is working correctly :-)
A bit surprising that it
gives _exactly_ the same result as the blockwise ranking in coin...
why? Without ties, the conditional and unconditional versions of
A quick reply,
are you looking for the matrix multiplication operator? In R it's %*%
Henrik
Christian Hinz wrote:
I have own function wrote. I used an algorithm, which was written in Matlab.
in matlab:
...
gamma = inv(v)*g;
...
#v = matrix of variable size, v=vv(k) = k=2 = dimension
Henrik Bengtsson [EMAIL PROTECTED] writes:
A quick reply,
are you looking for the matrix multiplication operator? In R it's %*%
... and solve(v,g) is probably preferable anyway.
my rewritten r-file:
...
gamma = solve(v)*g;
...
--
O__ Peter Dalgaard Ă˜ster
Hi,
I have two doubts about the nested aov output.
1) I have this:
anova.ratos - aov(Glicogenio~Tratamento+Error(Tratamento/Rato/Figado))
summary(anova.ratos)
Error: Tratamento
Df Sum Sq Mean Sq
Tratamento 2 1557.56 778.78
Error: Tratamento:Rato
Df Sum Sq Mean Sq F
How can I aggregate this data.frame to list the min and max date for
each unique id?
From this :
r = data.frame(id=rep(seq(1:3), 3), date= as.Date(c(rep(2005-08-25,3),
rep(2005-08-26,3), rep(2005-08-29, 3)), %Y-%m-%d))
r
id date
1 2005-08-25
2 2005-08-25
3 2005-08-25
1
Hello,
I guess a have a very simple problem though up to now couldn't solve it:
I want to plot two datasets wihtin one plot like plot(x) provides it for
one dataset(type=b that is: points connected by lines).
Example data 'x':
Befragung1 Befragung2 Befragung3 Geschlecht
2.25 2.34 1.78
maybe you could use something like this:
dat - data.frame(id = rep(1:3, 3), date = as.Date(rep(c(2005-08-25,
2005-08-26, 2005-08-29), each = 3)))
do.call(rbind, lapply(split(dat, dat$id), function(x) data.frame(id
= x$id[1], start = min(x$date), end = max(x$date
I
Dear All,
Can anyone give me some hints about how to set starting values for a lmer
model? For complicated models like this, good starting values can help the
numerical computation and make the model converge faster. Thanks!
Shige
[[alternative HTML version deleted]]
Dear Henrik, dear Peter
many thanks for your reply. I've now done tests on other computers and yes,
this is a local problem on my computer. Unfortunately, my second HD produces
the same flaw so it must be the disk controller or whatever. When correcting
the files with vedit, the non-standard
Since I have not seen a reply to this post, I will attempt a brief
comment: I won't comment on the specifics, but your general approach
seems appropriate. You already seem to know that a factor with k levels
is converted into (k-1) separate numerical variables, and a separate
I have installed package xtable with
su -c 'R CMD INSTALL xtable'
and got this promising feedback:
* Installing *source* package 'xtable' ...
** R
** data
** help
Building/Updating help pages for package 'xtable'
Formats: text html latex example
* DONE (xtable)
Despite that, R returns:
I was puzzled as to what the question actually was. If you set
options(contrasts=c(contr.sum, contr.poly))
the interaction contrasts are precisely those created manually. Nothing
fancier is required. But I am not sure what you want to do with them once
you have them.
On Tue, 30 Aug 2005,
Have you load the package with library(xtable) before attempting?
Andy
From: Mag. Ferri Leberl
I have installed package xtable with
su -c 'R CMD INSTALL xtable'
and got this promising feedback:
* Installing *source* package 'xtable' ...
** R
** data
** help
Building/Updating
Dear Prof. Ripley:
Agreed. I thought that coding the factor levels manually and / or
looking at model.matrix(y~xma) after further simplifying the problem
might help Scot clarify his question.
Thanks,
Spencer Graves
Prof Brian Ripley wrote:
I was puzzled as to
On 8/30/05, Shige Song [EMAIL PROTECTED] wrote:
Dear All,
Can anyone give me some hints about how to set starting values for a lmer
model? For complicated models like this, good starting values can help the
numerical computation and make the model converge faster. Thanks!
Shige
I agree
Wolfgang == Wolfgang Waser [EMAIL PROTECTED] writes:
I would be most obliged for any comments and help.
Wolfgang,
I've used R's fft to filter ECG signals and will comment on your
commentary based on my experience. First, as an easily
accessible reference, I suggest The Scientist and
XLSolutions Corporation (www.xlsolutions-corp.com) is proud to
announce 2-day R/S-plus Fundamentals and Programming
Techniques in New York City.
www.xlsolutions-corp.com/Rfund.htm
New York City September 22nd-23rd, 2005
Reserve your seat now at the early bird
On 8/30/05, Karsten Rincke [EMAIL PROTECTED] wrote:
Hello,
I guess a have a very simple problem though up to now couldn't solve it:
I want to plot two datasets wihtin one plot like plot(x) provides it for
one dataset(type=b that is: points connected by lines).
Example data 'x':
Befragung1
1) Do you really need to turn rows 3, 6, 9 of mat into NA ? I would
suggest the use of usediff() to remove them after complete.cases
2) The trick is to create a matrix of the same dimension as 'mat' but
initialised with NAs followed by replacing the required rows.
Example :
w - which(
You can shorten the notation by=0 .
On Tue, 2005-08-30 at 09:30 -0400, roger bos wrote:
Uwe,
Thanks. I have used merge a lot, but I didn't realize you could merge on
row names. I got it working now:
test - merge(out.r[1:5, 1:3], out.p[1:5, 1:3], by=row.names, all.x=TRUE)
Thanks,
Roger
Dear R wizards:
Has anyone gotten the TeXtext font encoding to work? If I execute:
if (is.null(postscriptFonts()$lucida)) {
luafmfiles - c(/usr/share/texmf/fonts/afm/yandy/lubright/lbr.afm,
/usr/share/texmf/fonts/afm/yandy/lubright/lbd.afm,
Dear all,
I am trying to use the
rstable(n, alpha, beta, gamma = 1, delta = 0, pm = c(0, 1, 2)))
function to generate positive stable random numbers. For positive stable
distribution, beta==1 and alpha is in (0,1), which defines random variables
with support (0, infinity). So, I used
Dear list.
New to R, I'm looking for a way of using crosstab to output low-dimensional
(higher than 2) contingency tables (frequencies, per-cents by rows, % by
columns, mean, quantiles) I'm looking for something of the following sort
dataframe: singers,
categorical variates: voice
On Tue, 2005-08-30 at 12:28 -0700, Isotta Felli wrote:
Dear list.
New to R, I'm looking for a way of using crosstab to output
low-dimensional (higher than 2) contingency tables (frequencies,
per-cents by rows, % by columns, mean, quantiles) I'm looking for
something of the following
X. Cong [EMAIL PROTECTED] writes:
Dear all,
I am trying to use the
rstable(n, alpha, beta, gamma = 1, delta = 0, pm = c(0, 1, 2)))
function to generate positive stable random numbers. For positive stable
distribution, beta==1 and alpha is in (0,1), which defines random variables
with
From where did you get rstable()? It isn't part of R itself, and
poking around shows at least three different sources for it
(CircStats, fBasics, stable (Lambert/Lindsey, not on CRAN)). I seem to
recall that the third one also has a distribution function which is
almost, but not quite
X. Cong wrote:
Dear all,
I am trying to use the
rstable(n, alpha, beta, gamma = 1, delta = 0, pm = c(0, 1, 2)))
function to generate positive stable random numbers. For positive stable
distribution, beta==1 and alpha is in (0,1), which defines random variables
with support (0, infinity).
Hi,
How can I get the p-value, S.E., and Z from the result of lrm fit? Thank you.
Donglei
__
R-help@stat.math.ethz.ch mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide!
[EMAIL PROTECTED] wrote:
Hi,
How can I get the p-value, S.E., and Z from the result of lrm fit? Thank you.
Donglei
That's automatic. The print method for lrm in the Design package does
that. Better though is the anova method which gives you the right
pooled tests when you have
Dear R wizards: I believe some more font encoding info. some of the
font encodings work, others do not: IsoLatin1, MacRoman, WinAnsi, and
PDFDoc seem fine. AdobeStd, AdobeSym, ISOLatin2, ISOLatin9, and
TeXtext seem broken, in that the resulting output file is silently
corrupt. The font
Hi Vincent,
Thanks for your good work!
Best,
--
Jose Claudio Faria
Brasil/Bahia/UESC/DCET
Estatistica Experimental/Prof. Adjunto
mails:
[EMAIL PROTECTED]
[EMAIL PROTECTED]
[EMAIL PROTECTED]
tel: 73-3634.2779
__
R-help@stat.math.ethz.ch mailing
Solved the problem with the 'R' installation!
One of the processes that loads up on startup was
conflicting with it. It was a program called
Gameutil.exe which seems to enhance games by
altering display modes. For some reason a conflict
existed. Anyhow I dont really need the Gameutil.exe,
so I'm
Hi Isotta,
You can do this with the reshape package (available from CRAN). eg
install.packages(reshape)
library(reshape)
data(singer, package=lattice)
singer$type - c(drammatic, spinto, lirico-spinto, lirico,
leggero)[sample(1:5, 235, replace=T)]
singer$school - c(german, italian, french,
Since I have not seen a reply to this, I will offer a brief comment.
I first tried, 'RSiteSearch(panel regression)', which produced many
irrelevant hits. Then I tried, 'RSiteSearch(analysis of panel data)'.
The fourth item on that list was
From 'RSiteSearch(polynomial distributed lag)', I learned that
there was a discussion of this issue not quite two years ago, which may
interest you. In brief, Thomas Lumley provded code with some noted
limitations.
spencer graves
[EMAIL PROTECTED] wrote:
Dear r-help
Dear Rexperts,
I intend to burn some R CDs to colleagues in Vietnam. I want to put all
binary files for base as well as contributed packages (for both Windows
and Linux). It is very time consuming if I download files by files. Is
there a place a can buy a CD with all the mentioned files?
G'day Nam-Ky,
NKN == Nam-Ky Nguyen [EMAIL PROTECTED] writes:
NKN I intend to burn some R CDs to colleagues in Vietnam. I want
NKN to put all binary files for base as well as contributed
NKN packages (for both Windows and Linux).
I don't believe that precompiled binary files for
62 matches
Mail list logo