similar to: Inverting a large Matrix (14000 x 14000)

Hi.., is there a function in R to fit bivariate splines ? I came across 'polymars' (POLSPLINE) and 'mars' (mda) packages. Are these the one to use or are there other specific commands? Thanks. Harsh

Hello R-Users, I have to invert a matrix 3000X3000 and the solve method doesn't work or it is too slow. Are there any methods to invert a big matrix? Regards Stefan --------------------------------- [[alternative HTML version deleted]]

How to invert match parameter? Like this: tc filter add dev eth0 protocol ip parent 1: prio 1 u32 \ match ip src 10.20.30.40 \ match ip dst !10.30.0/24 \ match ip dst !10.40.0/24 \ flowid 1:20 _______________________________________________ LARTC mailing list / LARTC@mailman.ds9a.nl http://mailman.ds9a.nl/mailman/listinfo/lartc HOWTO: http://lartc.org/

Hello, I am trying to generate random survival times by inverting the function,  S(t)= exp(b*F(t)), where b is constant and F(t) is some cumulative distribution function, let say that F(t) is cdf of normal distribution or any others distributions.   as we know that S(t) has uniform distribution on  (0,1) so we can write that U= exp(b*F(t)), where U is uniform (0,1). Now to generat the time t, we

Full_Name: Ravi Varadhan Version: 2.8.1 OS: Windows Submission from: (NULL) (162.129.251.19) Inverting a matrix with solve(), but using LAPACK=TRUE, gives erroneous results: Here is an example: hilbert <- function(n) { i <- 1:n; 1 / outer(i - 1, i, "+") } h5 <- hilbert(5) hinv1 <- solve(qr(h5)) hinv2 <- solve(qr(h5, LAPACK=TRUE)) all.equal(hinv1, hinv2) #

Is there an easy way to invert a table? (not to solve for the inverted matrix, just swap rows for columns & vice versa). I've gone through my data manipulation bible (Phil Spector's book), but to no avail. [[alternative HTML version deleted]]

Hi.. is there a limit on the number of explanatory variables in nls ? i have a dataframe with the columns names x1,x2..,x300 when i run nls it gives the error: " x181 not found" thought it does run when i have x1,x2,...,x170 variables. Thanks Harsh --------------------------------- [[alternative HTML version deleted]]

Refiling this. The actual fix was slightly more complicated. Will soon be committed to R-Patched (aka 2.9.1 beta). -p rvaradhan at jhmi.edu wrote: > Full_Name: Ravi Varadhan > Version: 2.8.1 > OS: Windows > Submission from: (NULL) (162.129.251.19) >=20 >=20 > Inverting a matrix with solve(), but using LAPACK=3DTRUE, gives erroneo= us > results: Thanks, but there seems

Hi, is there a lm which will implement the moore-pensrose generalized inverse. Thanks. Harsh --------------------------------- [[alternative HTML version deleted]]

Hi.. i was trying to integrate the indicator funtion but had problems when limits where negative or equal to the indicator condition my function is ________________________ fun1<-function(x){ as.numeric(x>=2) } _________________________ which should be Ind(x>=2)*x seems to work for the following two cases

Hi, I have a large number of parameters to estimate in nls say 100: beta1--beta100 lets say i have 100 values in a vector is there a way where i can create the start vector for nls using a loop instead of individually filling the 100 values. Thanks Harsh --------------------------------- [[alternative HTML version deleted]]

Hi All, when i run cv.lars(x3,y3) it runs fine. but when i run cv.lars(x3,y3,fraction=seq(0,0.1,100)) I get the following error. Error in apply((y[omit] - fit)^2, 2, mean) : dim(X) must have a positive length Any help/suggestions will be appreciated. Thanks. Harsh --------------------------------- [[alternative HTML version deleted]]

Hi.. I a have a symmetric matrix to plot . I would like to plot only the Upper triangle but with the diagonal as the Base of the rectangle. Is there an easy way to do it. Thanks. Harsh --------------------------------- [[alternative HTML version deleted]]

Full_Name: David Pleydell Version: 2.2.1 OS: Debian Etch Submission from: (NULL) (193.55.70.206) There appears to be a conflict between the chol functions from the Matrix and the SparseM packages. chol() can only be applied to a matrix of class dspMatrix if SparseM is not in the path. with gratitude David > library(Matrix) > sm <- as(as(Matrix(diag(5) + 1), "dsyMatrix"),

rvaradhan at jhmi.edu wrote: > Full_Name: Ravi Varadhan > Version: 2.8.1 > OS: Windows > Submission from: (NULL) (162.129.251.19) > > > Inverting a matrix with solve(), but using LAPACK=TRUE, gives erroneous > results: Thanks, but there seems to be a much easier fix. Inside coef.qr, we have coef[qr$pivot, ] <- .Call("qr_coef_real", qr, y, PACKAGE =

Hi.. i have an expression of the form: model1<-nls(y~beta1*(x1+(k1*x2)+(k1*k1*x3)+(k2*x4)+(k2*k1*x5)+(k2*k2*x6)+(k3*x7)+(k3*k4*x8)+(k3*k2*x9)+(k3*k3*x10)+ (k4*x11)+(k4*k1*x12)+(k4*k2*x13)+(k4*k3*x14)+(k4*k4*x15)+(k5*x16)+(k5*k1*x17)+(k5*k2*x18)+(k5*k3*x19)+

I have been running into difficulties with dispatching on an S4 class defined in the SparseM package, when the method calls are inside a function passed as the f= argument to optimize() in functions in the spdep package. The S4 methods are typically defined as: setMethod("det","matrix.csr", function(x, ...) det(chol(x))^2) that is within setMethod() rather than by name before

Full_Name: Julien Gagneur Version: 2.2.1 OS: Linux (Suse 9.3) Submission from: (NULL) (194.94.44.4) Using the SparseM library (SparseM_0.66) and the e1071 library (e1071_1.5-12) I fail using svm method with a sparse matrix. Here is a sample example. I experienced the same problem under Windows. > library(SparseM) [1] "SparseM library loaded" > library("e1071")

I'm having trouble loading the package SparseM in R 1.8.1, OS = Windows 2000. Installing appeared to go well; I saw no error messages, html documentation was installed, and "installed.packages()" lists SparseM among the installed packages. When I try to load the library, however, I get the following: > library(SparseM) Error in slot(mlist, "argument") : Can't get

Can R invert a 160000x160000 matrix with all positive numbers? Thanks a lot!