similar to: R-help Digest, Vol 71, Issue 16

Hi, I'm trying to make a chromosomal map in R by using the locus. I have a list of genes and their locus, and I want to visualise that so you can see if there are multiple genes on a specific place on a chromosome. A example of what I more or less want is below: http://www.nabble.com/file/p21474206/untitled.JPG untitled.JPG The genes and locus are here:

Gracias Carlos. Yo también he visto el ejemplo que te pone chatGPT, pero la salida que te da no soy capaz de interpretarla. Os paso las ordenes y las respuestas de R de la propuesta de chatGPT Ejemplo.aov<- aov(P~TRAT+CORTE+REP) > summary (Ejemplo.aov) Df Sum Sq Mean Sq F value Pr(>F) TRAT 6 0.0028 0.00046 0.777 0.590 CORTE 2 0.5022 0.25110 424.542 <2e-16

Buenas tardes. Tengo un ensayo de dos factores (Variedad y dosis) Quisiera hacer anovas para cada una de las dos variedades que estoy evaluando. Lo puedo hacer cogiendo los datos primero de una variedad y luego de la otra. Pero creo que hay un comando "by variedad" que lo hace si necesidad de andar modificando el archivo de datos. Pero no sé dónde hay que colocarlo. Si alguno me puede

Buenas, En R, como en la mayoría del software estadístico, no se utiliza ningún nivel de confianza sino que lo que se calcula es el p-valor asociado al contraste. De forma que cuanto más cerca de 0 esté el p-valor "menos credibilidad le damos a la hipótesis nula". Dicho mejor, debemos rechazar la hipótesis nula si el p-valor está por debajo de nuestro nivel de confianza. Por ejemplo,

Un saludo: Estoy intentando con el comando "detach()" eliminar todos los archivos que durante la analíticas voy utilizando y "attached" al programa. El problema es que usando detach () no me borra todo lo adjuntado. No sé si estoy usando mal el comando o estoy usando un comando incorrecto que no es para esto. Gracias. Juan Bautista Relloso Barrio Técnico del Dpto

Hola a todos: Preguntar si alguien sabe si se puede conectar el Tinn-R con la versión 2.15 de "R" Yo lo he intentado y no lo consigo. si abro un script con el tinn-R no se me activan los botones para poder correr las lineas y que el resultado aparezca en la ventana de abajo de "R". Vamos como pasaba con la versión 2.14 Sin más un saludo a todos y gracias por

Hola, ¿qué tal? attach y detach son funciones cuyas páginas de ayuda están llenas de avisos y caveats además de sugerencias acerca de alternativas a su uso. Yo siempre recomiendo no utilizarlas y yo en particular, las uso, si acaso, de manera sumamente excepcional. Un saludo, Carlos J. Gil Bellosta http://www.datanalytics.com El día 21 de noviembre de 2013 13:18, Juan Bautista <jbautista

Hello, I am calculating this thing with vectors (b) and matrices (G,P): b'G/sqrt(b'Pb) where the denominator is a quadratic form and therefore always a scalar. In Scilab, it is quite simple: b'*G/sqrt(b'*P*b) However, in R, the denominator is an (1x1)matrix and R claims it is non conformable and I have to use drop() or as.numeric(). Like this: > b = 1:2 > G=diag(1,2) >

I've got two tables.... first one(table1): ID chrom start end Ex1 2 152 180 Ex2 10 2000 2220 Ex3 15 3000 4000 second one ( table2): chrom location name 2 160 Alv 2 190 GNN 2 100

Greetings, I am new to R, but trying to put in the time to learn. I have read the R manual and several other introductory texts; however, there is nothing like actually putting it into practice. So here is my problem, and its more of a learning exercise for myself than anything else, but I'm stuck and getting extremely frustrated that I can't figure it out. I'm trying to make a

Hi Jim I am trying to prepare a bed file to load as accustom track on the UCSC genome browser. I have a data frame that looks like the one below. > x V1 V2 V3 1 chr1 11255 55 2 chr1 11320 29 3 chr1 11400 45 4 chr2 21680 35 5 chr2 21750 84 6 chr2 21820 29 7 chr2 31890 46 8 chr3 32100 29 9 chr3 52380 29 10 chr3 66450 46 I would like to insert the following 4 lines at the beginning:

Hello, I am a total beginner with ?R? and found a package ?venn? to create a venn diagram. The problem is, I cannot create the vectors required for the diagram. The manual say: "R> venn(accession, libname, main = "All samples") where accession was a vector containing the codes identifying the RNA sequences, and libname was a vector containing the codes identifying the

Dear list, I wish to plot chromatin precipitation data: I would like to have a rectangles (x:end-start, y:peak) but I do not have an idea how to define x (in terms of qplot syntax) and to choose the correct geom. mydata is a subset of a larger file. > mydata chrom start end peak 1 chr11 5291000 5291926 8 2 chr11 10988025 10988526 7 3 chr11 11767950 11768676 8 4

Hi, I'm trying to match peaks between chromatographic runs. I'm able to match peaks when they are chromatographed with the same method, but not when there are different methods are used and spectra comes in to play. While searching I found the ALS package which should be usefull for my application, but I couldn't figure it out. I made some dummy chroms with R, which mimic my actual

Dear all, I have data.frame object in R. I want to export it in tab-delimited file with several lines of header initiated with comment sign (#). I do not know how to do that in R. Could you please give helps on this problem? Thanks in advance. Best, Jian-Feng, ################################################################## The lines I want to write in the header lines look like, with words

Hi, i have two files (file1.txt and file2.txt) which i would like to merge, based on certain criteria, i.e. it combines data based on matching geneID and exons. i have used the merge option, but it does not give me the desired outcome. merged.txt shows the result i would like. *File1. txt* ** AffyProbe ProbeType Flag GeneSymbol GeneID Exons Chrom Strand Affytart AffyEnd 1

Hi, i have two files (file1.txt and file2.txt) which i would like to merge, based on certain criteria, i.e. it combines data based on matching geneID and exons. i have used the merge option, but it does not give me the desired outcome. merged.txt shows the result i would like. *File1. txt* ** AffyProbe ProbeType Flag GeneSymbol GeneID Exons Chrom Strand Affytart AffyEnd 1

I have a data frame (narrow) with 431 rows and 6 columns containing information on chromosome, position, lod1, lod2, lod3, lod4, looking like this: > narrow chr pos lod1 lod2 lod3 lod4 1 1 3.456 -0.025 -0.003 -0.209 -0.057 2 1 5.697 -0.029 -0.005 -0.200 -0.058 3 1 8.434 -0.049 -0.012 -0.247 -0.092 4 1 9.466 -0.074 -0.025 -0.300 -0.136 5 1 9.706

Hello all, I have written a for loop to act on a dataframe with close to 3million rows and 6 columns and I would like to pass it to apply() to speed the process up (I let the loop run for 2 days before stopping it and it had only gone through 200,000 rows) but I am really struggling to find a way to pass the arguments. Below are the loop and the head of the dataframe I am working on. Any hints

Dear R users, I need to merge two data frames based on both equal and unequal comparisons. The "sqldf" package used to work well , but today, I cannot resolve the following error by reinstallation of the sqldf package. Can anyone suggest a different way to perform this kind of merge function? Thank you, Ding > DMRlog2pbde47DMS <- sqldf("select * from DMR_log2pbde47 as a