similar to: setting pairwise comparisons of columns

Hi, I have a number of genes (columns) for which I want to examine pairwise associations of genotypes (each row is an individual)...For example (see data below), I would like to compare M1 to M2, M2 to M3, and M1 to M3 (i.e. does ac from M1 tend to be found with bc from M2 more often than expected.) Down stream I will be performing chi square tests for each pair. But I am looking for a way to

Dear R experts, For quite some time I have been trying to solve a mistery of reading a seemingly trouble-free text file. The data is temperature reconstruction arranged as a huge grid, preceded by seven "header lines" (which you see better if file is opened in Firefox or Chrome). This is the data (gridded temperature reconstruction)

Hardware: blade HS20 on RHEL4 U1 base kernel. # hg tip changeset: 6811:a8edb2655c5d220b1bbf2b8d8415aae0ac136da3 tag: tip user: cl349@firebug.cl.cam.ac.uk date: Wed Sep 14 10:47:48 2005 summary: Add methods to read/write Vm/Domain store entries. I created domU and got the messages: # xm create -c vm2.cfg Using config file "vm2.cfg". Started domain

Noob alert! Hi, I am having trouble figuring out how to filter a select helper. The following shows all of the rows in the table but I would like to filter the collection from the list_name column (:list_name => "Project Status"): @project_status = Valuelist.all.collect {|s| [s.list_value,s.id]} Thanks in advance, James -- Posted via http://www.ruby-forum.com/. -- You received

Dear all, apologies for this (perhaps recurrent) question but I did not found a question when searching mailing lists. How to write a list of a simple kind, e.g.: abc <- list(one=(1:2), two=(1:5)) # to a file? I understand that write() & co. cannot work but when I try sink("aa.txt", append=T, split=T) abc sink() # the output is indeed "split by rows" in the

Hello, I was wondering who to contact to see about somehow interconnecting or integrating the programs R and Epi Info. I have been involved in the development of Epi Info which is a software program developed by the Centers for Disease Control and Prevention (CDC) in Atlanta (www.cdc.gov/epiinfo<http://www.cdc.gov/epiinfo>). Epi Info allows for creating data entry systems and performs

?? Hola. Estoy intentando ordenar el gráfico tal que las columnas se apilen en el siguiente orden: debajo, los EPI; los SemiEPI en el medio y los NoEpi arriba de todo. Estos son mis datos: > tabla.barras Rangos Clase Asistentes 1 De 45d EPI 405 2 De 1a EPI 812 3 De 2a EPI 1639 4 De 3a EPI 2202 5 De 4a EPI 1994 6 De 5a EPI

Hola, Mira la versión de ggplot2 que estás usando por tenerla actualizada a la última... A mi me funciona introduciendo algún cambio salvo la última línea cuando cambias el tema... #----------- Lines <- " id Rangos Clase Asistentes 1 45d EPI 405 2 1a EPI 812 3 2a EPI 1639 4 3a EPI 2202 5 4a EPI 1994 6 5a EPI

The name conflicts in Epi and ROC packages (2 'ROC' functions are the problem) cause the following code to work once, but not twice: library(MASS); data(cats); x = cats[,2] y = ifelse(cats[,1]=='F',0,1) library(Epi); ROC(x,y,grid=0)$AUC library(ROC); AUC(rocdemo.sca(y, x, dxrule.sca)) What is the standard way of resolving name conflicts? Ask maintainers to resolve

Dear all, I would like to calculate the optimal cut off (threshold) of a test using the Epi package. Here I am presenting some data based on the output of two tests. I am interested in identifying the optimal cut off an its 95% CI. Running the ROC() function with the Epi package I obtain a nice picture that returns what I interpret as the optimal cut off with lr.eta=0.431. would be this the

This is a really simple question, I'm sure,but I can't make EpiR work! I keep getting the following: > epi.2by2(47, 263483, 282, 935028, method="cohort.time", conf.level=0.95) Error in epi.2by2(47, 263483, 282, 935028, method = "cohort.time", conf.level = 0.95) : unused argument(s) (935028) and I really don't know why!. Any ideas very very welcome thank

Hello: Does anybody knows if Epi Info, from CDC (www.cdc.gov/epiinfo) runs on Wine? I did try AppDB first but there are no entries for Epi Info. Epi Info is a data base for epidemiology and public-heath developed by Centers for Disease Control and Prevention (USA) and heavilly used in the world. It is based on Access 97 or 2000 witch I found in Frank's Corner AppDB. Thank you, Gerson

Hello all! I'm a newcomer to R, so please bear with me. My problem is the following: I'm trying to use R to do some geostatistics stuff (mostly kriging), and am having trouble printing a multiple figure. When I create a single figure (one image/contour command), it works fine. The instant I try to create a multiple figure file, I start getting postscript errors. I'm using the

I have a question, I generated some data and tried survreg to analyze them, but here, I want to get the covariance matrix of the coefficients. I know how to get corr from survreg, say, summary(f,corr=T), but how to get covariance matrix of the coefficient ? can I put them in a matrix and output ? thanks below is what I did <mailto:r-help@stat.math.ethz.ch> n=100; beta=1.0

There is now an Epi-package on CRAN. It is intended for epidemiological analysis in R. It has its own homepage, http://www.pubhealth.ku.dk/~bxc/Epi The package has been used at the course "Statistical practise in Epidemiology with R", see http://www.pubhealth.ku.dk/~bxc/SPE. A mailing list for people using R for epidemiology, R-sig-Epi, has been set up, see

Dear all, I would like to calculate the optimal cut off (threshold) of a test using the Epi package. Here I am presenting some data based on the output of two tests. I am interested in identifying the optimal cut off and its 95% CI. Running the ROC() function with the Epi package I obtain a nice picture that returns what I interpret as the optimal cut off with lr.eta=0.491. would be this the

There is now an Epi-package on CRAN. It is intended for epidemiological analysis in R. It has its own homepage, http://www.pubhealth.ku.dk/~bxc/Epi The package has been used at the course "Statistical practise in Epidemiology with R", see http://www.pubhealth.ku.dk/~bxc/SPE. A mailing list for people using R for epidemiology, R-sig-Epi, has been set up, see

The stat.table function in the Epi package won't do standard deviations. It didn't seem that it would be difficult to add an "sd" function to the stat.table function. Following the example for the mean, I set up a similar function for the sd (and included it as an options) but it just won't work. (I tried sending messages to the Epi mailing list after subscribing but my mail

Hi all, I have a question and I do not know If I am at the right place to ask this question. But is there someone that has the formula of CKD-Epi in code in R? I have tried a lot of loops, but none of the approaches give me the right answer. Is there someone who has this formula coded? Thank you! [[alternative HTML version deleted]]

A new major upgrade of the Epi package for Epidemiological data analysis has been put on CRAN, it is now at version 1.0.7. It contains an entirely new way of representing follow-up data on multiple timescales and multiple states. See the function Lexis(). Plus a lot of other useful stuff for epidemiological analysis. See more on the package homepage, www.biostat.ku.dk/~bxc/Epi Note also that