Learn how easy it is to sync an existing GitHub or Google Code repo to a SourceForge project! See Demo

Close

Diff of /inst/unscale_parameters.m [000000] .. [afe5a1] Maximize Restore

  Switch to unified view

a b/inst/unscale_parameters.m
1
# Copyright (C) 2003,2004  Michael Creel <michael.creel@uab.es>
2
# 
3
# This program is free software; you can redistribute it and/or modify
4
# it under the terms of the GNU General Public License as published by
5
# the Free Software Foundation; either version 2 of the License, or
6
# (at your option) any later version.
7
# 
8
# This program is distributed in the hope that it will be useful,
9
# but WITHOUT ANY WARRANTY; without even the implied warranty of
10
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
11
# GNU General Public License for more details.
12
# 
13
# You should have received a copy of the GNU General Public License
14
# along with this program; if not, write to the Free Software
15
# Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA  02111-1307  USA
16
17
# Unscales parameters that were estimated using scaled data
18
# primarily for use by BFGS
19
function [theta_us, vartheta_us] = unscale_parameters(theta, vartheta, scalecoefs);
20
  k = rows(theta);
21
  A = nth(scalecoefs, 1);
22
  b = nth(scalecoefs, 2);
23
  
24
  kk = rows(b);
25
  B = zeros(kk-1,kk);
26
  B = [b'; B];
27
  
28
  C = A + B;
29
  
30
  # allow for parameters that aren't associated with x
31
  if (k > kk)
32
      D = zeros(kk, k - kk);
33
      C = [C, D; D', eye(k - kk)];
34
  endif;  
35
  
36
  
37
  
38
  theta_us = C*theta;
39
  vartheta_us = C * vartheta * C';
40
endfunction