search for: weihenstephaner

Hi! In a current project, I am fitting loess models to subsets of data in order to use the loess predicitons for normalization (similar to what is done in many microarray analyses). While working on this I ran into a problem when I tried to predict from the loess models and the data contained NAs or NaNs. I tracked down the problem to the fact that predict.loess will not return a value at all

Dear R experts, Suppose I have a data frame of three variables: > foo <- data.frame(row=1:5, col=1:3, val=rnorm(15)) > foo row col val 1 1 1 -1.00631642 2 2 2 0.77715344 3 3 3 0.17358793 4 4 1 -1.67226988 5 5 2 1.08218836 6 1 3 1.32961329 7 2 1 -0.51186267 8 3 2 -1.20990127 9 4 3 -0.57786899 10 5 1 0.67102887 11 1 2

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:

Dear list, Can somebody tell me how to invoke user-defined functions from script files during run-time? Basically I have (almost) one function per script file. I am thinking of something like #include in C++ import in Java/Python use in Perl No, I don't need my functions every time R is started. Neither I am thinking about writing building my functions into packages. I just

Hi, I was wondering whether there was an easy way to change the attributes of all elements in a data.frame (rather than looping through elements)? Specifically, I would like to set the "dim" attributes to NULL Thanks for any help, Jarrod -- The University of Edinburgh is a charitable body, registered in Scotland, with registration number SC005336.

Hi, under the new version R 2.7.1 I can't use WinEdt. It seems that there is a bug in that version. To use older versions of the package RWinEdt doens't bring success. The Use of R 2.7.0 doesn't cause that Problem. Regards, Hans-J?rgen Auinger -- ---------------------------------------------- Hans-J?rgen Auinger Lehrstuhl f?r Pflanzenz?chtung Wissenschaftszentrum Weihenstephan

Hi all, I have a large data.frame, 1530 observation with 6 columns. I want to merge a 7th column, a transformation of the response variable (hospital admissions), namely trans<-sqrt(copd$admissions+0.25) trans<-data.frame(trans) And now when I do copd2<-merge(copd,trans) (copd being my original data.frame), R either crashes or is taking an extremely long time to do the computation. I

Dear All, I'm trying to read a text data file that contains several records separated by a blank line. Each record starts with a row that contains it's ID and the number of rows for the records (two columns), then the data table itself, e.g. 123 5 89.1791 1.1024 90.5735 1.1024 92.5666 1.1024 95.0725 1.1024 101.2070 1.1024 321 3 60.1601 1.1024 64.8023 1.1024 70.0593

I have Windows XP Professional Version 2002 and the R-Version 2.1.1. I want to do factor analysis with R. In Google Search I find the reference to the RScaLAPACK-Package, but my R-Version 2.1.1. told me that "RScaLAPACK is invalid package, bevor 2.0.0. installed" . Which package can I take for factor analysis? [[alternative HTML version deleted]]

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>>=

Dear ExpeRts, I am trying to plot a bunch of growth curves and would like to get some more control over groups and line colors than I seem to have. Example: # make some data dat <- Orange dat$group <- ifelse(dat$Tree%in%c('1','4','5'), 'A', 'B') # plot xyplot(circumference~age, dat, groups=group) # now use lines to make the growth curve more

Dear ExpeRts,t I am trying to read tab delimted data produced by somewhat brain dead software that seems to think it's a good idea to have an extra tab character after the last column - except for the header line. As explained in the help page, read.delim now assumes that the first column contains the row.names (which is not even wrong) but now and all col.names get shiftet by one column.

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, I have a string vector, say, x <- "a", "aab" and anther string vector, say, y <- "a", "aa", "aab", "aabc". Is there any R function to get the indices of y for the elements of x, that is, foo(x, y) will give me the index vector c(1, 3)? I know i can combine apply() and grep("^aab$", y) to do it. But is there any

Hello everyone, does anybody know how to increase the decimal places that R uses to calculate something. I think that in default R uses 6 decimal places but I need 12. Thank you very much! Best regards, Maja! --

Dear R-helpers I would like to do a fit of the form: y = a (x+c)**b, where a, b and c are unknown. Does anybody know how to do it? Thanks Thomas

All, I think this is pretty basic but I couldn't find the answer in any source. I have just built my own R for Linux (amd 64). It runs well, but R thinks that it only has 80 characters or so of screen width. I log on with Putty to the box and can stretch it as large as I want. Is there a simple way to tell R to use more screen width when displaying? From whatami: OS RELEASE :

Hi, * I am using read.table command as follow kegg<-read.table("c:/IDs.tab",header =TRUE,quote= "'", sep="\t") * Fragment of file is as follow: ID Pathway 04916 Melanogenesis 04920 Adipocytokine signaling pathway 04930 Type II diabetes mellitus 04940 Type I diabetes mellitus 04950 Maturity onset diabetes of the young 05010