mixed effects model in R

Book
http://lme4.r-forge.r-project.org/book/

Slides
http://www.stat.wisc.edu/~bates/PotsdamGLMM/LMMD.pdf

awk: verify same number of columns

## prints number of columns if it changed from previous line

> awk -F"\t" 'NF != oldnf { print NF; oldnf = NF }' filename |head

|head added just in case file is very long

## using a script

> awk -f script.awk filename | head

script.awk:

BEGIN {
 FS="\t"
 oldnf = 1
}
{
  if(NF != oldnf) 
  {
 print NF
 oldnf = NF
  }
}

Other examples in 
http://sparky.rice.edu/~hartigan/awk.html

merge eps files into one

convert *.eps all.pdf

running R64 in aquamacs

http://www.matthewckeller.com/html/aquamacs.html

USING 64-BIT R FROM WITHIN AQUAMACS

1) Aquamacs is terrific because it comes bundled with about everything you need to 
start using R with it. No need to separately download ESS (emacs speaks statistics) 
or mess with the .ess-lisp files. Download Aquamacs from http://aquamacs.org/ and 
install it.

2) Now any file with an “.R” extension that is opened with aquamacs will 
automatically have syntax highlighting AND will allow you to start 64-bit R. I 
highly recommend checking out this page, http://ess.r-project.org/, and downloading 
the ESS reference card therein if you are new. Emacs has a bit of a learning curve, 
but you’ll begin loving its capabilities as an aid to programming very soon.

3) To start 64-bit R from within aquamacs hit Control-u then Alt-x then “R” then 
. Aquamacs will ask you for “Starting Args”. Here’s where you tell it you 
want to run 64-bit R. Type “--arch=x86_64” then . A new window will open up 
that is the R session.

4) That’s almost it. Just one more issue: graphing will crash (something to do with 
the new R and X11). Aquamacs has X11 as the default (null) graphics device. To 
change this, type in R: “options(device=”quartz”)”. Now your null graphics device is 
quartz and it should all work out. If you really want to use X11 to write graphics, 
you can use x11(type="Xlib")' for each new graph. What I do is to set 
“options(device=”quartz”) within my Rprofile.site file.

trim white spaces

check ?gsub