s-news
[Top] [All Lists]

Re: Matrix Problems

To: Lewis Jordan <lcj6444@owl.forestry.uga.edu>
Subject: Re: Matrix Problems
From: Tony Plate <tplate@blackmesacapital.com>
Date: Wed, 17 May 2006 11:52:08 -0600
Cc: s-news@lists.biostat.wustl.edu
In-reply-to: <s46b2328.055@redhorse.forestry.uga.edu>
References: <s46b2328.055@redhorse.forestry.uga.edu>
User-agent: Mozilla Thunderbird 1.0.5 (Windows/20050711)
It's unlikely that you have sufficient memory for your result, even if you are using a 64-bit version of S-PLUS: your matrix takes about 44Gb of storage.

The message you are seeing is due to integer overflow in the internal calculations, e.g.:

> 76365 * 76365 # Can't fit this result in a 32-bit integer
[1] NA
> # how many Gb to store a double matrix this size?
> 76365.0 * 76365 * 8.0 / 2^30
[1] 43.4489
>

-- Tony Plate

Lewis Jordan wrote:
I am "attempting" to perform some matrix multiplication with s-plus 6.0.
I have 2 matrices with the following dimensions: #> dim(C)
#[1] 76365  15
#> dim(delta)
#[1] 15 15
I am want to do the following:

a<-C%*%solve(delta)%*%t(C)

which should be a 76365 X 76365 matrix. However, I get the following error: "Problem in rep.int: No missing values allowed in "times": rep.int(1:n, times) " I have checked the data, and there are no "missing values". Is this what this message even means? Any suggestions? Thanks, Lewis Jordan ***************************************
Lewis Jordan
Research Coordinator:  Wood Quality Consortium
University of Georgia
Warnell School of Forest Resources
Athens, Georgia 30602
Office:  706-542-7587
Cell:     706-614-8884
Email:  lcj6444@forestry.uga.edu <mailto:lcj6444@forestry.uga.edu>
***************************************


<Prev in Thread] Current Thread [Next in Thread>