search for: imhof

Hi all. I'm trying to numerically invert the characteristic function of a quadratic form following Imhof's (1961, Biometrika 48) procedure. The parameters are: lambda=c(.6,.3,.1) h=c(2,2,2) sigma=c(0,0,0) q=3 I've implemented Imhof's procedure two ways that, for me, should give the same answer: #more legible integral1 = function(u) { o=(1/2)*sum(h*atan(lambda*u)+sigma^2*lambda*u/(1+l...

...asking the community: Is there a way to fix this or do I have to code around it? This is in R 2.11.1 (Linux), by the way. Thanks in advance Philipp -- Dr. Philipp Pagel Lehrstuhl f?r Genomorientierte Bioinformatik Technische Universit?t M?nchen Wissenschaftszentrum Weihenstephan Maximus-von-Imhof-Forum 3 85354 Freising, Germany http://webclu.bio.wzw.tum.de/~pagel/

Dear All: I do see I have weeks of codes in my console when I check with my arrow up keys. I have been clearing them with Control L but it seems to clear it clear the screen temporally. I do see the previous codes again when I open R the next day, after quitting the session! Q: How do I clear this? Thanks; YA -- View this message in context:

This may be asking too much, but I'm wondering if anyone has a solution (even a hack) for creating multiple (overlay) plots in an Sweave file and post-processing the overlays in beamer appropriately. For example, suppose I have a series of figure blocks in my .Rnw file: <<plot1,fig=TRUE>>= [stuff] @ <<plot2,fig=TRUE>>= [stuff] @ <<plot3,fig=TRUE>>=

...ts own color Of course I can now use the col argument to manually assign the colors by group but is there a more elegant way that I missed? cu Philipp -- Dr. Philipp Pagel Lehrstuhl f?r Genomorientierte Bioinformatik Technische Universit?t M?nchen Wissenschaftszentrum Weihenstephan Maximus-von-Imhof-Forum 3 85354 Freising, Germany http://webclu.bio.wzw.tum.de/~pagel/

...st. Did I just not find it or does it not exist? And if it doesn't - does anyone else think it would be a nice item for the wishlist? cu Philipp -- Dr. Philipp Pagel Lehrstuhl f?r Genomorientierte Bioinformatik Technische Universit?t M?nchen Wissenschaftszentrum Weihenstephan Maximus-von-Imhof-Forum 3 85354 Freising, Germany http://webclu.bio.wzw.tum.de/~pagel/

I have a 64-bit windows box - Intel Xeon CPU E7340 @ 2.4GHz 31.9GB of RAM I have R 2.11.1 (64bit) running on it. My csv data is 3.6 GB (with about 15 million obs, 120 variables.) ------------------------------------------------ I have successfully imported the data above into R. No problem. Now I am trying to run 'rpart' on my data. But I got the following error : Error: cannot

Dear all I have written a function to perform a very simple but useful task which I do regularly. It is designed to show how many values are missing from each variable in a data.frame. In its current form it works but is slow because I have used several loops to achieve this simple task. Can anyone see a more efficient way to get the same results? Or is there existing function which does

Hi, I have been working with R for the past couple of years; analyzing data and producing some graphics. I was just wondering if people use graphics from R straight into articles or are they always edited in some way; fonts, headers, axis, color etc? Using photoshop or some other programs? I would like to think it is possible, better and more profession to do it all in R. I tried google and

Hi everybody, I have a strange behavior of the as.Date() function. For example: as.Date(36525, origin="1900-01-01'") I would expect to get 2000-01-01. But R gives me [1] "2070-01-01" I tried with several origins but the results are still the same. It seems that the origin does not apply and remains 1970-01-01. Does anyone can give me a hint how to

Hello All, How do I save the output of the R console and the graphic output to the same PDF file and append these to each other? I need to have a frequency table and a corresponding graph, one below the other in a file. I have tried with sending the cross table to the graph window using 'textplot' and then saving the graphic output. However, the table does not look nice in the graph

Dear Sir/Madam; I'm not sure whether this is the correct contact for help. I've been recently working with R on my project, unfortunately It sudenly crashes! It gives me the following message: "FATAL ERROR: unable to restore saved data in .RDATA" I decided to uninstall the copy (a R2.11.0) and installed a new version (2.11.1) but I'm still receiving the same message. When I

Hello, I have seen a few posts about changing a vector of numbers to days of the week where the numbers are 1-7. I have a similar problem with a twist. I have a vector of numbers ranging form 20 to 500. The numbers represent days since a starting point. The list is not consecutive, some numbers skipped and some numbers duplicated. I know day 1 was a Monday. I want to use this vector in a lm

Hi guys, I am trying to shift down the "ylab" of my plot but can't find how to do it. I tried to tune "mar" but it enable more room for the labels to be displayed but it does not move to "ylab" as I would like. Is there a way with "par" to shift down my "ylab" ?? Cheers -- View this message in context:

Hi all, Is there a way to save the currently displayed plot to an image format just after we view it? I think this would be more intuitive as a user if I wish to save it just after I visualize the plot. I am aware that we need to do some thing like this jpeg('somefilename.jpg') ... plot... commands... dev.off() KM [[alternative HTML version deleted]]

Hello, I'd like to fit a logit function to my data. The data is distributed like a logit (like in this plot on wikipedia http://en.wikipedia.org/wiki/File:Logit.png) but the values on the x-axis are not between 0 and 1. I don't think using a glm is the solution because I simply want to infer the parameters of the logit function (offset, compression, slope...), so I can apply it to all

Hi Group, I am having a strange problem for this simple task. I am using read.table to read a plan 3 column CSV file . the file is getting read . But the first column has datetime in the csv file in the following format: 20110221.114041 But this is being read as 20110221 only . the time portion (decimal is missing) in the data frame Any idea why? Thanks, Santosh

Hi all, When I attempt to run a script, I keep getting the error message shown below. I have seen from searching that the "Use of save versions prior to 2 is deprecated" could be an error resulting from a permissions problem, but I have been unable to find anything describing the "file has magnic number '# Coh'". The file I am attempting to read into the script is a

Dear R helpers Suppose I have a vector as vect1 = as.character(c("ABC", "XYZ", "LMN", "DEF")) > vect1 [1] "ABC" "XYZ" "LMN" "DEF" I want to reverse the order of this vector as vect2 = c("DEF", "LMN", "XYZ", "ABC") Kindly guide Regards Vincy [[alternative HTML

dear lady and gentalmen: i am gaoshan from kansas university. i used such coding to deal with gel data data <- ReadAffy() Warning messages: 1: In file(out, "wt") : cannot open file 'C:\Users\gaoshan\AppData\Local\Temp\RtmpvsyXOV\Rhttpd3f0b2e85': No such file or directory 2: In file(out, "wt") : cannot open file