Example usage for org.apache.commons.math3.linear ArrayRealVector getDimension

Introduction

In this page you can find the example usage for org.apache.commons.math3.linear ArrayRealVector getDimension.

Prototype

@Override
public int getDimension()

Source Link

Usage

From source file:gamlss.distributions.GA.java

/** Second cross derivative.
 * @param y - vector of values of response variable
 * @return  a vector of Second cross derivative
 *//*from   w ww . j  ava  2s .c  o m*/
private ArrayRealVector d2ldmdd(final ArrayRealVector y) {
    //d2ldmdd = function(y)  rep(0,length(y)),
    // all elemnts are 0's
    return new ArrayRealVector(y.getDimension());
}

From source file:gamlss.distributions.GA.java

/**
 * Set tempV2 = sigma*sigma.// w w  w  .  j ava  2 s  .c  om
 * @param y - response variable
 */
private void setInterimArrays(final ArrayRealVector y) {

    muV = distributionParameters.get(DistributionSettings.MU);
    sigmaV = distributionParameters.get(DistributionSettings.SIGMA);
    size = y.getDimension();
    tempV2 = sigmaV.ebeMultiply(sigmaV);
}

From source file:gamlss.distributions.TF.java

/** Second cross derivative of likelihood function in 
    * respect to mu and sigma (d2ldmdd = d2l/dmu*dsigma).
    * @param y - vector of values of response variable
    * @return  a vector of Second cross derivative
    *//*w  ww.jav a 2 s  .  co  m*/
private ArrayRealVector d2ldmds(final ArrayRealVector y) {
    //d2ldmdd = function(y) rep(0,length(y)
    return new ArrayRealVector(y.getDimension());
}

From source file:gamlss.distributions.TF.java

/** Computes the global Deviance Increament.
 * @param y - vector of response variable values
 * @return vector of global Deviance Increament values 
 *//*from   w w  w  .  j  a  v a2s. com*/
public final ArrayRealVector globalDevianceIncreament(final ArrayRealVector y) {
    //G.dev.incr  = function(y,mu,sigma,nu,tau,...)  
    size = y.getDimension();
    double[] out = new double[size];

    double[] muArr = distributionParameters.get(DistributionSettings.MU).getDataRef();
    double[] sigmaArr = distributionParameters.get(DistributionSettings.SIGMA).getDataRef();
    double[] nuArr = distributionParameters.get(DistributionSettings.NU).getDataRef();

    for (int i = 0; i < size; i++) {

        out[i] = (-2) * dTF(y.getEntry(i), muArr[i], sigmaArr[i], nuArr[i], Controls.LOG_LIKELIHOOD);
    }
    return new ArrayRealVector(out, false);
}

From source file:gamlss.distributions.TF.java

/** Calculate and set initial value of nu.
 * @param y - vector of values of response variable
 * @return vector of initial values of nu
 *///from w ww.j  av  a 2s  .  c  o  m
private ArrayRealVector setNuInitial(final ArrayRealVector y) {
    //nu.initial = expression( nu <- rep(10, length(y)))
    tempV = new ArrayRealVector(y.getDimension());
    tempV.set(10.0);
    return tempV;
}

From source file:gamlss.distributions.TF2.java

/** Computes the global Deviance Increament.
 * @param y - vector of response variable values
 * @return vector of global Deviance Increament values 
 *///from  w w  w .ja va  2 s  .c o  m
public final ArrayRealVector globalDevianceIncreament(final ArrayRealVector y) {
    //G.dev.incr  = function(y,mu,sigma,nu,tau,...)  
    size = y.getDimension();
    double[] out = new double[size];

    double[] muArr = distributionParameters.get(DistributionSettings.MU).getDataRef();
    double[] sigmaArr = distributionParameters.get(DistributionSettings.SIGMA).getDataRef();
    double[] nuArr = distributionParameters.get(DistributionSettings.NU).getDataRef();

    for (int i = 0; i < size; i++) {

        out[i] = (-2) * dTF2(y.getEntry(i), muArr[i], sigmaArr[i], nuArr[i], Controls.LOG_LIKELIHOOD);
    }
    return new ArrayRealVector(out, false);
}

From source file:gamlss.distributions.PE.java

/** Second cross derivative of likelihood function in 
 * respect to mu and sigma (d2ldmdd = d2l/dmu*dsigma).
 * @param y - vector of values of response variable
 * @return  a vector of Second cross derivative
 *//*from   w  ww .  j  a v  a2s  .  c  o m*/
private ArrayRealVector d2ldmds(final ArrayRealVector y) {
    // d2ldmdd = function(y)  rep(0,length(y))
    return new ArrayRealVector(y.getDimension());
}

From source file:gamlss.distributions.PE.java

/** Second cross derivative of likelihood function
  * in respect to mu and nu (d2ldmdd = d2l/dmu*dnu).
 * @param y - vector of values of response variable
 * @return  a vector of Second cross derivative
 *///from w w w. ja v a 2s  .  c o  m
private ArrayRealVector d2ldmdn(final ArrayRealVector y) {
    //d2ldmdv = function(y)  rep(0,length(y)),
    return new ArrayRealVector(y.getDimension());
}

From source file:gamlss.distributions.TF2.java

/** Calculate and set initial value of nu.
 * @param y - vector of values of response variable
 * @return vector of initial values of nu
 *///from   w w  w  . j a  v  a 2s  . c  om
private ArrayRealVector setNuInitial(final ArrayRealVector y) {
    //nu.initial = expression( nu <- rep(4, length(y)))
    tempV = new ArrayRealVector(y.getDimension());
    tempV.set(4.0);
    return tempV;
}

From source file:gamlss.distributions.PE.java

/** Calculate and set initial value of nu.
 * @param y - vector of values of response variable
 * @return vector of initial values of nu
 *//*from  w w w. j av a  2s  . co  m*/
private ArrayRealVector setNuInitial(final ArrayRealVector y) {
    //nu.initial = expression( nu <- rep(1.8, length(y)))
    tempV = new ArrayRealVector(y.getDimension());
    tempV.set(1.8);
    return tempV;
}