similar to: Error in using matest function

Hi, I am using "*Maanova* package" to do anova. I have created *datafile* with probeID as the first column, which is a tab limited text file and also created *designfile*. I have created *readma object* which is named as abf1. >From that readma object, i have to create data object by using *createData*function and also i hav to create model object by using *makemodel* function,

Hi, I am using "*Maanova* package" to do anova. I have created *datafile* with probeID as the first column, which is a tab limited text file and also created *designfile*. I have created *readma object* which is named as abf1. >From that readma object, i have to create data object by using *createData*function and also i hav to create model object by using *makemodel* function,

I have just installed R 2.4.0 (Windows XP OS) and have tried to run code that had worked with previous versions of R. However, I now get an error message: Error in La.svd(iCmat, method = "dgesvd") : unused argument(s) (method = "dgesvd") The R NEWS page about the release and issue 6/4 of the newsletter states: La.svd(method = "dgesvd") is defunct. As a

Generator Microsoft Word 11 (filtered medium) Hi, I'm the mannova package maintainer. We used La.svd(X, method="dgesvd") in maanova package before. After R-2.3.0, the old La.svd() method was deprecated for option method="dgesvd". I changed maanova code correspondingly, which will call method="dgesdd" instead. But after that, we keep getting below error message

Dear All, Under Debian GNU/Linux La.svd (with method = "dgesdd") sometimes gives the error "Error in La.svd(data, nu = 0, nv = min(nrow, ncol), method = "dgesdd") : LAPACK routine DGESDD gave error code -12" It seems not to depend on the data per se, but on the relationship between numbers of rows and columns. For example, if the number of columns is 100,

Dear R helpers, I am working with R 2.4.1 GUI 1.18 (4038) for MacOSX. I have a matrix of 10 000 genes and try to run the following commands: > model.mix<-makeModel (data=data, formula=~Dye+Array+Sample+Time, random=~Array+Sample) > anova.mix<-fitmaanova (data, model.mix) > test.mix<-matest (data, model=model.mix, term="Time", n.perm=100, test.method=c(1,0,1,1))

Hi, Am using maanova package for doing anova.But am getting error like this..plz, help me regarding this.. > TGR=read.madata("rmaexpr.dat",designfile="design.dat") Reading one color array. Otherwise change arrayType='twoColor' then read the data again Warning messages: 1: In read.madata("rmaexpr.dat", designfile = "design.dat") : Assume that

Hi, I'm just wondering if it is planned to include the Lapack routine DGESDD (and friends) in R-1.4.0? This is faster (supposedly by a factor of ~6 for large matrices) than DGESVD which is currently (R-1.3.1) called by La.svd. And if it is not in the plans yet, is there a chance it could be? I've added it to my local version of R-1.3.1 and so far see a factor of 4 improvement over

Full_Name: Chanseok Park Version: 1.8.1 OS: linux-gnu Submission from: (NULL) (130.127.112.183) digamma with any negative value does not give a right answer. It gives -1.797693e+308 for any negative arguments. For example, digamma(-1.1) gives -1.797693e+308. The right answer should be 10.15416 This bug can be easily fixed by using the following digamma identity. digamma(x) = digamma(1-x) -

I have two matrix s1 and s2, each of them is 1000*1. and I have two equations: digamma(p)-digamma(p+q)=s1, digamma(q)-digamma(p+q)=s2, and I want to sovle these two equations to get the value of x and y, which are also two 1000*1 matrices. I write a program like this: f <- function(x) { p<- x[1]; q <- x[2]; ((digamma(p)-digamma(p+q)-s1[2,]) )^2 +((digamma(q)-digamma(p+q)-s2[2,]) )^2

Full_Name: Henrik Bengtsson Version: 1.5.1 OS: WinMe Submission from: (NULL) (217.210.0.243) In the Perl script $R_HOME/bin/check there is a bug under the section "Check R code for syntax errors" where the 'Rfiles <- c(...)' is build up. If there are too many files in @Rfiles the source code line generated will be too long and weird things will happen, e.g. strange

Hey, I'm trying to build something like this http://flowingdata.com/2011/05/11/how-to-map-connections-with-great-circles/ but with my own data in csv files. The code runs well if I use the same csv files as the author, but with mine , this is what I get *Code* library(maps) library(geosphere) map("world") xlim <- c(-180.00, 180.00) ylim <- c(-90.00, 90.00)

On Wed, 14 Aug 2002, Henrik Bengtsson wrote: > Sorry, but it was indeed the redirection of the standard output in > Cygwin/bash that cause the first problem, not R (I should stop doing > troubleshooting at 1:00 AM). So please forget about the problems reported in > R_CMD_check.out. However, it would still be nice if you still update R CMD > check to do join with "\n".

Is there a way to make a shading interpolation on an image plot? Something similar to matlab 'shading interp', I think it is called Gouraud shading. What I want is to make a image plot look nicer. with image() it looks very facetted, and I would like to make it look smoother. I've tried with interp.surface() in fields package but it (obviously) makes nan values at the borders and

This message is in MIME format. The first part should be readable text, while the remaining parts are likely unreadable without MIME-aware tools. Send mail to mime@docserver.cac.washington.edu for more info. ------=_NextPart_000_0006_01C2432C.7ACD6BA0 Content-Type: TEXT/PLAIN; CHARSET=iso-8859-1 Content-Transfer-Encoding: 8BIT Content-ID: <Pine.GSO.4.44.0208140920312.15226@auk.stats>

Given a time series of length N, I am trying to figure out its maximum DWT decomposition level. With reference to "dwt" function of R package "wavelets", running the provided example I get the following and wonder how the maximum decomposition level (which is not an integer number) is calculated. I tried myself the formula in the on-line documentation but could not get the

Hi there, Could somebody help me disect this legacy R script I inherited at work, I have two questions: 1. I've tried to upgrade our R version from 1.6.2 (yeah, I know), to R 2.0, but some of the lines in this script are not compatible with R 2.0, could someone help me figure out where the problem is? 2. the jpeg generated (attached) seems to be off on some of the data, is there a better way

Bueno, ya lo he solucionado gracias a Carlos me ha enviado un correo en privado (supongo que se ha despistado, si no lo querías hacer público ya es tarde) con esta info: -- Creo que lo tienes (en formato básico) aquí: https://stat.ethz.ch/pipermail/r-help/2007-November/146285.html Aunque pensaba que si tu objetivo último es el de que usuarios "potencialmente tontos" sean capaces de

Hello everyone I'm a beginner in Stats and R, I'm using R 2.10.1. I need to create a multivariate qq plot, there is 8 variable group with each has 55 number of input. An example of what I did so far, just to get my point out: > data=read.csv(file.choose(),header=T) > data country village group av_expen P2ary_ed no_fisher 1 Cook Islands Aitutaki D

I have a classification tree model similar to the following (slightly simplified here): > treemod<-tree(y~x) where y is a factor and x is a matrix of numeric predictors. They have dimensions: > length(y) [1] 1163 > dim(x) [1] 1163 75 I?ve evaluated the tree model and am happy with the fit. I also have a matrix of cases that I want to use the tree model to classify. Call it