Displaying 2 results from an estimated 2 matches for "cg10000".
2012 Jul 18
3
Mean of matched data
...value (of gene expression, in
the column called "value") for that row's gene&gender combination, and in
the other column the standard deviation for the gene&gender combination.
Any suggestions?
Rob
Example of the top of the data frame:
gene variable value gender line rep
1 CG10000 X208.F1.30456 4.758010 Female 208 1
2 CG10000 X365.F2.30478 4.915395 Female 365 2
3 CG10000 X799.F2.30509 4.641636 Female 799 2
4 CG10000 X306.M2.32650 4.550676 Male 306 2
5 CG10000 X712.M2.30830 4.633811 Male 712 2
6 CG10000 X732.M2.30504 4.857564 Male 732 2
7 CG10000 X707.F1.31120 5.104165 Female...
2012 May 14
1
Extract Variance Components
...ne,beta$line,beta$residual)
rownames(beta)<-gene.list[[1]] #gene.list is just a list of all my focal
genes/
Produces this:
row.names # sexline #line
#residual #varcomp #gene
1 (Intercept) #0.005087154 #0.007013756 #0.01839387 #0.6031809
#CG10000
I have thought about creating a loop to then produce my big full data frame
but I can't get it to work (nor can my friend at work who is really good
with R) we have tried numerous ways and are now at our wits end...
/## try to loop for all genes...
newobj <- matrix(nrow=0, ncol=3)...