search for: bioinformatiker

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:

Hi all, I have to calculate an expression using acos -function. A strange result of acos appears *1. case with error* ss <- sin(10.74*pi/180)**2 +(cos(10.74*pi/180)*cos(10.74*pi/180)*cos(0*pi/180)) ss acos(ss) [1] NaN Warning message: NaNs produced in: acos(ss) *2. case without error* ss <- sin(10.7*pi/180)**2 +(cos(10.7*pi/180)*cos(10.7*pi/180)*cos(0*pi/180)) ss acos(ss)

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

Hi, How can I overlap two density plots? A <- c(8,10,10,11,11,11,12,12,12,12,11,11,11,10,10,7) B <- c(11,13,13,14,14,14,15,15,15,15,14,14,14,13,13,10) plot(density(A)) plot(density(B)) Regards, Emre P.S. There's a overlap.density function in package DAAG Is there another way to do this? -- --- Emre ÜNAL http://www.geocities.com/dusemre [[alternative HTML version

Greetings helpRs -- I would like to use plot() to plot two cumulative distribution curves so that a user of the plot can compare the distributions of the two variables. The following code draws two distributions separately, but I cannot find the instruction necessary to add a second cumulative distribution to the first one. Any suggestion would be very welcome. x1 <-

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.

Hello, I have difficulties in reading in a text file with multiple (=2) header lines. My table looks like this: id code gr grcode AA AB AC AD time 30 40 50 60 1 1234 1 c 0 1 0 0 2 1346 1 c 0 0 0 1 3 1456 1 c 0 0 1 1 4 1893 1 c 0 0 1 1 Can someone help me out? Thanks Deborah --

Hi Rowland, > Remove 'map untrusted to domain = Yes', it has been removed. > > Add 'domain logons = Yes' > > This gets it back to being a PDC: thanks for the hints! I did that, but it doesn't help. I guess the problem is not on the PDC server but on SERVER2. That's the one that got upgraded and stopped working (even with the non-pdc config of SERVER1). I

Hi list, I have a question on how to identify the overlapped items from two vectors: >x<-c(0,1,2) >y<-c(1,2,2,3) and plot the number of the overlapped as well as non-overlapped in a diagram? Thanks much, Allen [[alternative HTML version deleted]]

Dear R Community, A simple problem (for some of you): I wish to index a data.frame by all elements NOT in my index E.g.: > a<-as.data.frame(matrix(rnorm(100),nrow=10,ncol=10)) > b<-which(a$V1>0.8) > b [1] 1 4 6 10 > a_indexb<-a[b,] > a_notIndexB<-a[!b,] > nrow(a_notIndexB) [1] 0 Indexing a on b is not a problem (a_indexb), but how can do get only the

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

Hi all, I have a data management question. I am using an panel dataset read into R as a dataframe, call it "ex". The variables in "ex" are: id year x id: a character string which identifies the unit year: identifies the time period x: the variable of interest (which might contain NAs). Here is an example: > id <- rep(c("A","B","C"),2)

Rowland penny via samba wrote: > So your server doesn't seem to be able to find winbindd, are you sure it > is running ? > > What does this show: > > ps ax | grep '[w]inbind' > > What OS is this ? It's SuSE Linux Enterprise 15sp1. winbindd is definitely running, I showed that in the first mail in the output of "rcwinbind status", there you can

Dear I try to add three different lines (solid, short dash, long dash) in to current barplot. I saw there are types of lines ("p","b", ect). However, they are not what I request. Has nayone have this experience? Many Thanks! Xin [[alternative HTML version deleted]]

Hallo, I was trying some code, but couldn't make one step of the code properly. Can anybody please help me? I have one matrix like this > values [,1] [,2] [,3] [,4] [,5] [1,] 0.7777778 0.36111111 0.22222222 0.1388889 0.0000000 [2,] 1.0000000 0.00000000 0.53846154 0.0000000 0.5384615 [3,] 0.5200000 0.48000000 0.64000000 0.0000000 0.8800000 [4,] 0.8928571