similar to: Duplicated genes

Hi all, is there a function (or variable) I can use in a custom init script that identifies the init script name? i.e. I'm porting some init scripts from gentoo, where the $SVCNAME variable identifies the init script name within the script itself... d /* Davide Cittaro Cogentech - Consortium for Genomic Technologies via adamello, 16 20139 Milano Italy tel.: +39(02)574303007 e-mail:

Hi all, I'm new to Centos, just moved here from Gentoo Linux. I have to install a server for bioinformatics purposes and I see that default yum repositories do not include any bioinformatics software (i.e. ncbi-toolkit, blat, and others). I'm googling a bit but I can't find a valuable solution: which is (or which are) the best repository I should add to have a satisfying list

Hey All, I have a dataset(dat1) like this: ObsNumber ID Weight 1 0001 12 2 0001 13 3 0001 14 4 0002 16 5 0002 17 And another dataset(dat2) like this: ID Height 0001 3.2 0001 2.6 0001

Hi all, I'm migrating from Gentoo to CentOS... I'm experiencing a rather low performance in NFS r/w (as client). NFS server is solaris (which exports zfs volumes via nfs). The very same exports were mounted with the same parameters (auto,nosuid,exec) on gentoo and centos server (bot x86_64)... It happens that centos is 5-10 times slower either in read and write operations... Ok,

Hi we have samba 3.0.14a on FreeBSD 5.4. We tried with different kind of locking and oplocks (both enabled and disabled). If we try to copy from a Windows XP client a file larger than 3g, we get these error: Cannot copy XXX. The specified network name is no longer available. We traced this problem in the logs (log level 10) and we got this error ------------ [2006/09/19 10:29:41, 5]

Hi all Say I have a function: myname=function(dat,x=5,y=6){ res<<-x+y-dat } for various input such as myname(dat1) myname(dat2) myname(dat3) myname(dat4) myname(dat5) how should I modify the 'res' line, to have new informative variable name correspondingly, such as dat1.res dat2.res dat3.res dat4.res dat5.res stored in the workspace. This is only an example of a complex

I was modifying a source to add OpenMP capabilities. I knew that gcc supports OpenMP since version 4.2, so I've compiled a recent (4.3.2) gcc and installed in my home directory, this because gcc shipped with CentOS 5.2 is 4.1. I've then noticed that my executable links /usr/lib64/libgomp.so, installed with gcc 4.1 from CentOS rpms... Is CentOS default compiler OpenMP ready? Does

Hi all, I have a data matrix likein "input.txt" 8 9 2 5 4 5 8 5 6 6 8 9 2 8 9 2 8 9 2 1 8 9 2 5 4 5 8 5 6 4 8 9 2 5 4 5 8 5 6 6 8 9 2 8 9 2 8 9 2 1 8 9 2 5 4 5 8 9 2 2 In this example will be an 6x10 matrix (or data frame) I want to detect how many times in a row appears this combination 8 follewd by 9 followed by 2, and create a new matrix with only this number of occurs then

Dear all, I'm calling Winbugs (1.4.3) through R2WinBugs (2.1-18 coda_0.14-7) to fit a switching random walk model, but come up with an instant trap with the log only displaying 'check('. I will paste the trap with session info below; I'd be very grateful for any ideas. Couple of leads: 1. I presume the problem relates to the r package itself or the way I call bugs(), because I

Hi, I have a pretty simple problem. Here is the code: dat1=complete(dat.mice,1) dat2=complete(dat.mice,2) dat3=complete(dat.mice,3) dat4=complete(dat.mice,4) dat5=complete(dat.mice,5) dat6=complete(dat.mice,6) dat7=complete(dat.mice,7) dat8=complete(dat.mice,8) dat9=complete(dat.mice,9) dat10=complete(dat.mice,10) dat11=complete(dat.mice,11) dat12=complete(dat.mice,12)

Dear R developers, I noticed a bug in the stats4 package, specifically in the confint method applied to ?mle? objects. In particular, when some ?fixed? parameters define the log likelihood, these parameters are stored within the mle object but they are not used by the ?confint" method, which retrieves their value from the global environment (whenever they still exist). Sample code: >

Hi Guys, in what kind of problem I will be if I try to transfer files from a server that has on it 3.0.9 on windows running under cygwin, and the client having 2.6.9 on Mac OS X ? I tried to do it and no problem occured, apparently.. but I'd like to be sure.. I used the option --vvrtW thx in advance -- Vito Mule' HelpDesk Technician - Area Ricerca Istituto Europeo di Oncologia Via

Hi all, Is there a way to dump the current packages installed on a machine and use it do install/uninstall packages on other machines? I've configured one at install time but to speed up other installations I would like to install default packages and then install/uninstall some starting from my first machine configuration... thanks d PS Oh, I'm going to RTFM too... but I'm

1. I'm attempting to test for Random Effects. I've grouped the data on subject (grid) but want to use lme to build the model without subject as a RE then add it and do anova between the 2 models. This is the result I get and it appears it's adding Random Effects. tmp.dat4 <- groupedData(Trials ~ 1 | grid, data = tmp.dat4) mod2a <- lme(Trials ~ factor(group_id) + reversal,

Hi everyone, I'm having a little trouble with working out what formula is better to use for a repeated measures two way anova. I have two factors, L (five levels) and T (two levels). L and T are both crossed factors (all participants do all combinations). So, I do: summary(aov(dat~L*T+Error(participant/(L*T)),data=dat4)) But get different results with:

Hi all: How can I calculate the cumulative variance (or variance for each component) in correspondence analysis? If were possible in ADE4 package Thank you -- Antonio Punzón Merino O__---- Instituto Español de Oceanografía c/ /'_ --- Centro Oceanográfico de Santander (*) \(*) -- Promontorio de San Martín S/N ~~~~~~~~~~ 39004-Santander; Spain PO BOX: 240 Tlf: +34 942 29 10 60 Fax: +34

Hello everybody: I get a strange behavior with seq, take a look at this: > msd <- seq(0.05,0.3, 0.01) > msd[13] [1] 0.17 > class(msd) [1] "numeric" > class(msd[13]) [1] "numeric" > typeof(msd[13]) [1] "double" now the problem: > msd[13] == 0.17 [1] FALSE It is strange only to me? Consider that: > 0.17 == 0.17 [1] TRUE and also > a

Dear R users, I would like to combine lattice plot (xyplot) and standard R plot (plot and plotCI) in an unique figure. I use the function "par()" to combine plot and plotCI and I use the function "print()" to combine xyplot. I tried to use these functions to combine xyplot and plotCI and plots but they do not work. Does anybody know how I can do this? Thank you very much in

1. This relates to the package *rhdf5* and its implementation of the HDF5 dataspace interface. I am asking for an example of how other people who use this package make use of the HDF5 data space interface exposed by the library. Longer answer: As per my understanding, the dataspace interface exposes data locations within a dataspace, but even while retrieving data from an hdf5 file using methods

Dear all, I have a question with respect to counting process formulation of the coxph(survival) model. I have two groups of observations for which I have partitioned each observation into two distinct time intervals, namely, entry day till day 13, and day 13 till death or censorship day (of course the latter only for the observations that survived the first 13 day interval), and added a