Create ProfiledDeviance.R

Cost function for a Linear Mixed Effects model. Requires explicit  input of X and Z (fixed and random effects) matrices and measurements vector y. 
This code is experimental and has known "shortcomings"; criticism will very appreciated.