Feng,
I had the same question as you, how to read a subset of data, and the same
reaction as Wensui when I discovered that read.table could not. Even if my
computer's memory were up to it, I am troubled by the idea of reading in 1.8
GB of data (in my case) to get just 4,000 numbers, for
A BIG thanks to Luke Tierney and Jim Holtman (and perhaps others? I haven't
seen the latest digest) for pointing out two fixes to some code that was
troubling me. I since found a third problem, and with that, have good code.
Their comments are summarized below this fixed version of the code.
Hello, The code below is supposed to be a wrapper for matplot to
do columnwise visible comparison of several matrices, but I'm
doing something wrong because I can't access an argument called
'colnum'. I'd be most grateful for some insight.
Thanks,
John Thaden
Little Rock, AR
Apparently you're right that colnum doesn't exist when it needs to
be evaluated, but why? Why is 'paste' being evaluated so early? It is,
after all, the value of an argument ('main') of my mmatplot function
with colnum being another argument. I thought arguments were lazy-loaded.
Does using
I'd like to use someattributes(), as described
in documentation for R version 2.4.1 (windows build)
help(attributes)
however, someattributes() does not seem to exist.
someattributes()
Error: could not find function someattributes
Is this true or am I doing something wrong?
-John
I had written
...someattributes() does not seem to exist.
And Haris Skiadas replied
My help shows it as moreattributes, not
someattributes. (MacOSX), though doesn't
sound like it should be platform-specific).
Thanks for correcting me. Actually, my windows
R documentation says
Haris Skiadas replied
Thanks for correcting me. Actually, my windows
R documentation says mostattributes(), but it
makes no difference -- none of the three show
up as function names or R objects.
That's because there is no mostattributes function,
it only works as an assignment:
I am pretty pleased with baselines I fit to chromatograms using the
runquantile() function in caTools(v1.6) when its probs parameter is
set to 0.2 and its k parameter to ~1/20th of n (e.g., k ~ 225 for n ~
4500, where n is time series length). This ignores occasional low-
side outliers, and,
Is it possible there is a typo in your attached code? You wrote...
loglik - function(theta, param) {
sum(d1*log(hazard1(x1, param)) - chazard1(x1,param) +
d2*log(hazard2(x2, theta, param)) - chazard2(x2, theta, param))
} / 1
Confidentiality Notice: This e-mail message, including any
Ive been converting elements of matrices and arrays, e.g., from
Integers to double-precision, by vectorizing the matrix and then
remaking it. Alternatively, I can redefine one element as double
which then redefines them all. Both methods are quick, so I guess
I shouldn't complain, but I would
## While in R v. 2.3.1 (the mid-July patch for Windows)
## on a Windows XP machine, this call to order() works fine...
order(1:10,decreasing = TRUE)
## [1] 10 9 8 7 6 5 4 3 2 1
## ...however, the argument name 'decreasing'
## must be typed in toto (note the missing 'g'
## in the
On Sunday, July 09, 2006 12:31 PM, Roger Koenker = RK
[EMAIL PROTECTED] wrote
RK On 7/8/06, Thaden, John J [EMAIL PROTECTED] wrote:
JT As there is nothing inherent in either compressed, sparse,
JT format that would prevent recognition and handling of
JT duplicated index pairs, I'm
has written. He
has a book that is scheduled for publication this September
DB Tim Davis (2006), Direct Methods for Sparse Linear Systems, SIAM,
Philadelphia, PA
DB I hope we will be able to refer to that book for details of the
representation and algorithms.
On 7/8/06, Thaden, John J [EMAIL
With thanks to Matrix package co-author Martin
Maechler, I'm happy to report satisfactory closure
of two recent threads I initiated about that package:
- Warning while subsetting with Matrix
- R Crash with 'library(Matrix);as(x,dgCMatrix)
In the first, I reported seeing a warning
In the Matrix package v. 0.995-11 I see that the dgTMatrix
Class for compressed, sparse, triplet-form matrices handles
Identically indexed data instances by summing their values,
e.g.,
library(Matrix)
(Mt - new(dgTMatrix,
i = as.integer(c(0,0,1,1,4)),
j = as.integer(c(0,1,2,2,4)),
x =
Martin Maechler replied to my query Warning while subsetting...:
MartinM JohnT == Thaden, John J [EMAIL PROTECTED]
MartinM on Thu, 6 Jul 2006 00:02:10 -0500 writes:
JohnT ...
JohnT # While subsetting x, I was surprised to get this warning:
JohnT y-x[1:300,]
JohnT Warning
Hello,
Could someone please explain the following warning while subsetting in
the Matrix package?
Thanks,
John Thaden, PhD
U. Arkansas for Med. Sci.
Little Rock AR USA
# In the Matrix package...
library(Matrix)
# ...I had previously created a sparse matrix in triplet form:
str(x)
Formal
Sorry, I omitted background information:
R version: 2.3.0
OS: Windows XP
CPU: Pentium III,
RAM: 768 MB
Also, what command-line memory settings might prevent R from crashing
while using the Matrix package to convert my 600 X 4482 dgTMatrix to the
dgCMatrix class or to an expanded
:
JJTh == Thaden, John J [EMAIL PROTECTED]
on Sun, 2 Jul 2006 07:38:21 -0500 writes:
JJTh Dear R-Help list: I'm using the Matrix library to
JJTh operate on 600 X ~5000 element unsymmetrical sparse
JJTh arrays. So far, so good, but if I find I need more
JJTh speed
Dear R-Help list:
I'm using the Matrix library to operate on 600 X ~5000 element
unsymmetrical sparse arrays. So far, so good, but if I find I need more
speed or functionality, how hard would it be to utilize other sparse
matrix toolsets from within R, say MUMPS, PARDISO or UMFPACK, that do
20 matches
Mail list logo