Learning R

Beginning at the beginning: R inferno - nice pdf, probably helps if you are already a programmer. -------------------------------------------------------------------------------------- Youtube! apply, effectively removes the need for a double for loop. tapply expects vectors. ----------------------- This is wrong!: mean(2, -100, -4, 3, -230, 5) This is correct!: mean(c(2, -100, -4, 3, -230, 5))

11/12/2012

If using SQLForge, need to have installed (but not necessarily loaded) human.db0 --------------------------------------------------------------------------------- plot(.., col=ifelse(foo>2, "red", "black") -------------------------------------------

The aim of the moderated t-test is to test whether the log-ratios are consistently greater or lesser than zero, not to test whether the log-ratios are equal.  The moderated t-test does not penalize a gene for having a large variance as much as the ordinary t-test would do.  Compared to the ordinary t-test, the moderated t-test gives more weight to the treatment fold change.  It can therefore be viewed as a compromise between the ordinary t-test and ranking by fold-change.
--------------------------------------------------------------------------------------------------------------
The basic difference between linear models and glms is that with linear models, the F-test and all the individual t-tests that make it up can be conducted as one computation, but for glms each individual 1df test and overall F-test have to all be separate computations.
--------------------------------------------------------------------------------------------------------------
Delta beta refers to change in Beta value relative to a group.
----------------------------------------------------------------------------------------------------------------
Pearson correlations are a marginal measure of association and thus sensitive to confounding factors.
------------------------------------------------------------------------------------------------------
Agilent arrays - can't use robust-spline normalization, only global loess normalization is meaningful.

22/11/2012

Because I always forget:
2^number
###################
(HAO) AND "Homo sapiens"[porgn:__txid9606]

19/11/2012

scvitic
spartak 25
[Note the space!]
##############

Tuesday, 13th November 2012

Bpipe
hs_er_pid4.log - is a Java crash dump (a 'hotspot error log')
See if something is still running:
bpipe jobs
go to the directory of the pipeline and type
bpipe log
bpipe run -n 8
########
48 cores, 264G of memory
#####################
pathognomonic - characteristic or symptomatic of a disease or condition.
-------------------------------------------------------------------------
CNV arrays:
For any future work, ideally batches of > 50.
Less than 10 is basically ineffective.
-----------------------------------
ClinVar aggregates information about sequence variation and its relationship to human health.  Because the resource is still under active development, our preliminary website release is limited to our preview site.

Thursday, 8th November 2012

grep -E "7$" hasoverlaps.bed.complete - get all lines ending in 7.
###################################################
Making a coloured .bed?
http://cloford.com/resources/colours/namedcol.htm
#######################################