Displaying 20 results from an estimated 190 matches similar to: "Correlation analysis"
2009 Jun 16
1
Error using getBM() to query BioMart archives
I'm trying to identify the positions of all genes within a specific
chromosomal region using biomart. When using the current biomart
database I'm able to do this without issue. However, I need to use
build 36 of the mouse genome which was last included in ensembl mart
46. I selected this mart and the mouse dataset as follows:
mart<-useMart(biomart="ensembl_mart_46",
2012 Sep 02
2
a newbie seeking for a simple problem
Dear Experienced R users,
I have a looks-like simple but complicated problem urgently needed to be solved. Below is the detail:
I have two dataframes, df1, df2. df1 contains two column and many thousands rows: column 1 is a "gene_name", column 2 is "value". df2 contains only one column which is "gene_name" with couple hundred rows. I want to change "value"
2008 Mar 25
2
help with rowsum/aggregate type functions
Hi--
This is a question with a trivial and obvious answer, I'm sure, but I can't seem to find it in the help files and books that I have handy. I have a dataframe consisting of two columns, "Gene_Name," a list of gene symbols, and "Number," a numeric measure of how frequently a tag representing that gene showed up in a SAGE library. Several of the genes are
2010 May 25
4
R eat my data
HI, Dear R community,
My original file has 1932 lines, but when I read into R, it changed to 1068
lines, how comes?
cdu@nuuk:~/operon$ wc -l id_name_gh5.txt
1932 id_name_gh5.txt
> gene_name<-read.table("/home/cdu/operon/id_name_gh5.txt", sep="\t",
skip=0, header=F, fill=T)
> dim(gene_name)
[1] 1068 3
--
Sincerely,
Changbin
--
Changbin Du
DOE Joint Genome
2011 Apr 07
1
Two questions about metacharacter in regexprs and function return
for the script, please kindly see the script below. At line 10 and line 13,
my problems occurs.
The first one is I try to retrieve the gene official name from a column of a
table. The pattern of official name is something starting with gene_name.
For detail problems, please see the according lines.
Any suggestions are appreciated
example of matching source (extract the Nnat, sometime it would
2008 Jun 16
1
Separator with " | " for read.table
Hi,
I have the following data file to be parsed and captured as a data frame:
__DATA__
#GDS_ID GENE_NAME GENE_DESCRIPTION GENE_FUNCTION
1007_s_at | DDR1 | discoidin domain receptor tyrosine kinase 1 | protein-coding
1053_at | RFC2 | replication factor C (activator 1) 2, 40kDa | protein-coding
117_at | HSPA6 | heat shock 70kDa protein 6 (HSP70B') | protein-coding
__END__
In particular it is
2010 Apr 07
2
help in attach function
Hi, r-community,
This morning, I MET the following problem several times when I try to attach
the data set.
When I closed the current console and reopen the R console, the problem
disappear. BUt with the time passed on, the problem occurs again.
Can anyone help me with this?
> attach(total)
The following object(s) are masked from total ( position 3 ) :
acid base cell_evalue
2008 Jun 16
2
Creating a Hash from Data.Frame
Hi,
I have the following data frame:
> print(mydatframe)
__DATAFRAME__
V1 V2 V3
1 1007_s_at DDR1 discoidin domain receptor tyrosine kinase 1
2 1053_at RFC2 replication factor C (activator 1) 2, 40kDa
3 117_at HSPA6 heat shock 70kDa protein 6 (HSP70B')
__END__
Is there a way to create a hash with
V2 as Key and V3 as its value?
- Gundala Viswanath
Jakarta - Indonesia
2006 Mar 21
3
ROWNAMES error message
I am getting an error message, which I do not know the source to.
I have a matrix SAMPLES that has preexisting rownames that I would like to
change.
GENE_NAMES contains these rownames.
> rownames(SAMPLES) = GENE_NAMES
Error in "dimnames<-.data.frame"(`*tmp*`, value = list(list(V1 = c(3843, :
invalid 'dimnames' given for data frame
> dim(SAMPLES)
[1] 12626
2005 Jul 07
1
Plotting Character Variable
Any ideas about the following problem:
I have a matrix (A) that looks like this:
gene_names values
hsa-mir-124 0.3
hsa-mir-234 0.1
hsa-mir-344 0.4
hsa-mir-333 0.7
..... .......
(This is a 2 by 22283 matrix: quite large)
I would like to plot the values, but output the gene_names as the plotting
symbol. I have tried
2016 Apr 05
0
Is that an efficient way to find the overlapped , upstream and downstream rangess for a bunch of rangess
I do have a bunch of genes ( nearly ~50000) from the whole genome, which
read in genomic ranges
A range(gene) can be seem as an observation has three columns chromosome,
start and end, like that
seqnames start end width strand
gene1 chr1 1 5 5 +
gene2 chr1 10 15 6 +
gene3 chr1 12 17 6 +
gene4 chr1 20 25 6 +
gene5
2009 Nov 20
1
how to specify the order of panels with xyplot
> chromosomes
id refseq name length
1 0 NC_000001.9 Homo sapiens chromosome 1 247249719
2 1 NC_000002.10 Homo sapiens chromosome 2 242951149
3 2 NC_000003.10 Homo sapiens chromosome 3 199501827
4 3 NC_000004.10 Homo sapiens chromosome 4 191273063
5 4 NC_000005.8 Homo sapiens chromosome 5 180857866
6 5 NC_000006.10 Homo sapiens chromosome 6
2004 Aug 06
1
questions related to ploting in R
Dear all.
I need to draw a scatter plot of 23 chromosome copy numbers (y axes) against
chromosome and physical location within each chromosome in one plot. The
data matrix looks as below:
chr location copy_num
1 118345 1.320118
1 3776202 1.133879
1 4798845 0.989997
1 5350951 1.100967
. more data here
.
.
2 118345 2.459119
2 157739 1.915919
2 1530065 1.924372
2
2009 Oct 21
0
drawing a line indicating extent of each factored data series in multipanel lattice xyplot
Hi,
Am am plotting aggregated frequency data (extracted from an RDBMS)
relating to DNA sequence features for each of the human chromosomes as
described in the table chromosomes below (the frequency data is in a
table 'hits' that has a value (or not) for each of a set of bins across
each chromosome).
I would like to mark the extent of the chromosome (according to the
length value in
2011 Apr 20
2
'Record' row values every time the binary value in a collumn changes
My question is twofold.
Part 1:
My data looks like this:
(example set, real data has 2*10^6 rows)
binary<-c(1,1,1,0,0,0,1,1,1,0,0)
Chromosome<-c(1,1,1,1,1,1,2,2,2,2,2)
start<-c(12,17,18,20,25,36,12,15,16,17,19)
Table<-cbind(Chromosome,start,binary)
Chromosome start binary
[1,] 1 12 1
[2,] 1 17 1
[3,] 1 18 1
[4,] 1
2009 Dec 04
1
Lattice : Help with changing the labels of x-axis in respective panels
Dear R-Helpers,
I am not very experienced in using lattice and I am still in the learning
stage
I have a data set which looks like this: (I have deleted a few lines in
order to save space)
Chromosome marker Marker.Name Distance
1 1 1 PeMm261 0.0000
2 1 2 Xtxp8 10.1013
..
20 1 20 EbMi148 210.3099
21 1 21 Xtxp25
2010 Jun 01
5
Help barplots
Dear All,
I am newbie to R, and I wanted to plot a barplots with R and in such a way
that It will also show me position which I can plot on the bar line.
Here is my code that I am using to plot,
> chromosome <- c(40.2, 35.6, 36.1, 29.6, 31, 29.6, 31, 29.4, 28.2, 23, 23,
28.2)
>barplot (chromosome, col="purple", xlab="Oryza sativa Chromosomes", border
= NA, space =
2007 Apr 02
3
Random number from density()
Hello,
I'm writing some genetic simulations in R where I would like to place
genes along a chromosome proportional to the density of markers in a
given region. For example, a chromosome can be presented as a list of
marker locations:
Chr1<-c(0, 6.5, 17.5, 26.2, 30.5, 36.4, 44.8, 45.7, 47.8, 48.7, 49.2,
50.9, 52.9, 54.5, 56.5, 58.9, 61.2, 64.1)
Where the numbers refer to the locations of
2008 Jul 09
2
replacing value in column of data frame
Dear all,
Probably a very basic question but I need some help.
I have a data frame (made by read.table from a text file) of microarray data, of which the first column is a factor and the rest of the columns are numeric.
The factor column contains chromosome names, so values 1 through 22 plus X, Y and XY. The numeric columns contain positions or intensity measurements.
What I need to do is
2012 Sep 10
3
plot: x and y chromosomes are missing
Hello,
I have a list with gene names, fold changes (=expression level) and
chromosomes.
Names fold change chromosome
hz 1.5 2
If I plot fold change versus chromosome (or vice versa):
plot (ch, fc)
I see only the chromosomes with numbers but not those with letter (x and
y). What can I do?
A second question:
How can I add a single line in that plot at a certain