Exemplo n.º 1
0
void tabulaterun(dataobject & o)
  {

  // datamatrix X;

  unsigned i;

  vector<ST::string> vnames = o.m.getModelVarnamesAsVector();

  // testing, wether all variables specified are already existing
  vector<ST::string> notex;
  bool failure=false;
  if ((o.allexisting(vnames,notex)) == false)
    {
    for (i=0;i<notex.size();i++)
      o.outerror("ERROR: variable " + notex[i] + " is not existing\n");

    failure = true;

    } // end: if ((datap->allexisting(modelvarnamesv,notex)) == false)

  ST::string ifexpression = o.methods[9].getexpression();


  if (failure==false)
    {

    // o.makematrix(vnames,X,ifexpression);

    // statmatrix<int> index(X.rows(),1);

    for(i=0;i<vnames.size();i++)
        {
//        o.m.makeModelMatrix_j(o.d,X,i);
        datamatrix Xi;
        o.makematrix(vnames[i],Xi,ifexpression);
        statmatrix<int> index(Xi.rows(),1);

        index.indexinit();

        Xi.indexsort(index,0,Xi.rows()-1,0,0);

        o.out("\n");
        o.out("Variable: " + vnames[i] + "\n",true, false, 14);

        o.out("\n");

        unsigned k = 0;
        double cum = 0;
        unsigned zaehler = 0;

        while(k<index.rows() && zaehler<101)
            {
            double anz=0;

            int* p = index.getV() + k;
            int* q = index.getV() + k;
            for(unsigned j=k;j<index.rows();j++,p++)
                 {
                 if (Xi.get(*p,0) == Xi.get(*q,0))
                    anz = anz+1;
                 }

                 k = k + anz;
                 zaehler = zaehler + 1;

            }

        if(zaehler>=101)
           o.outerror("ERROR: too many values");
        else
           {
           ST::string value = "Value";
           ST::string obs = "Obs";
           ST::string freqst = "Freq";
           ST::string cumst = "Cum";
           o.out(value.helpfill(15) + obs.helpfill(10) +
                 freqst.helpfill(15) + cumst.helpfill(15) + "\n",true);
           k = 0;

           while(k<index.rows())
                 {
                 double anz=0;

                 int* p = index.getV() + k;
                 int* q = index.getV() + k;
                 for(unsigned j=k;j<index.rows();j++,p++)
                    {
                    if (Xi.get(*p,0) == Xi.get(*q,0))
                       anz = anz+1;
                    }

                 double wert;
                 wert = Xi.get(*q,0);
                 double freq = anz / Xi.rows();
                 cum = cum + freq;

                 ST::string str_wert = ST::doubletostring(wert,8);
                 ST::string str_anz = ST::inttostring(anz);
                 ST::string str_freq = ST::doubletostring(freq,4);
                 ST::string str_cum = ST::doubletostring(cum,4);

                 o.out(str_wert.helpfill(15) + str_anz.helpfill(10) +
                       str_freq.helpfill(15) + str_cum.helpfill(15) + "\n");

                 k = k + anz;
                 }
           }

        }
    o.out("\n");
    }
  }
Exemplo n.º 2
0
void pctilerun(dataobject & o)
  {

  datamatrix X;

  unsigned i;

  vector<ST::string> vnames = o.m.getModelVarnamesAsVector();

  ST::string ifexpression = o.methods[10].getexpression();

  // testing, wether all variables specified are already existing
  vector<ST::string> notex;
  bool failure=false;
  if ((o.allexisting(vnames,notex)) == false)
    {
    for (i=0;i<notex.size();i++)
      o.outerror("ERROR: variable " + notex[i] + " is not existing\n");

    failure = true;

    } // end: if ((datap->allexisting(modelvarnamesv,notex)) == false)

  if (failure==false)
    {

    o.makematrix(vnames,X,ifexpression);

    o.out("\n");

    //Ausgabe:
    //ST::string varst = "Variable";
    //ST::string

    for(i=0;i<vnames.size();i++)
       {
 //       o.m.makeModelMatrix_j(o.d,X,i);
        o.out("Variable: " + vnames[i] + "\n",true, false, 14);
        o.out("\n");

        for(unsigned j=1; j<6; j++)
           {
           double quant = X.quantile(j,i);
           ST::string quantst = ST::doubletostring(quant,8);
           o.out(" " + ST::inttostring(j) + "%" + quantst.helpfill(15) + "\n");
           }
        unsigned j = 25;
        while(j<=75)
           {
           double quant = X.quantile(j,i);
           ST::string quantst = ST::doubletostring(quant,8);
           o.out(ST::inttostring(j) + "%" + quantst.helpfill(15) + "\n");
           j = j + 25;
           }
        for(unsigned j=95; j<100; j++)
           {
           double quant = X.quantile(j,i);
           ST::string quantst = ST::doubletostring(quant,8);
           o.out(ST::inttostring(j) + "%" + quantst.helpfill(15) + "\n");
           }
           o.out("\n");
        }
   }

 }
Exemplo n.º 3
0
void descriptiverun(dataobject & o)
  {

  // datamatrix X;

  unsigned i;

  vector<ST::string> vnames = o.m.getModelVarnamesAsVector();

  // testing, wether all variables specified are already existing
  vector<ST::string> notex;
  bool failure=false;
  if ((o.allexisting(vnames,notex)) == false)
    {
    for (i=0;i<notex.size();i++)
      o.outerror("ERROR: variable " + notex[i] + " is not existing\n");

    failure = true;

    } // end: if ((datap->allexisting(modelvarnamesv,notex)) == false)


  ST::string ifexpression = o.methods[8].getexpression();

  if (failure == false)
    {

    // o.makematrix(vnames,X,ifexpression);

    o.out("\n");
    ST::string varst = "Variable";
    ST::string obsst = "Obs";
    ST::string meanst = "Mean";
    ST::string medst = "Median";
    ST::string stdst = "Std";
    ST::string minst = "Min";
    ST::string maxst ="Max";
    o.out(varst.helpfill(11) + obsst.helpfill(10) + meanst.helpfill(15) + medst.helpfill(15)
           + stdst.helpfill(15) + minst.helpfill(15) + maxst.helpfill(15) + "\n",true);

    for(i=0;i<vnames.size();i++)
      {
//      o.m.makeModelMatrix_j(o.d,X,i);
      datamatrix Xi;
      o.makematrix(vnames[i],Xi,ifexpression);      // für jede Variable eigene Matrix erzeugen?

      unsigned int obs = Xi.rows();
      double mean = Xi.mean(0);
      double var = 0;
      if(obs > 1)
        var = Xi.var(0) * obs / (obs-1);

      ST::string str_std;
      if (var>0)
        {
        double std = sqrt(var);
        str_std = ST::doubletostring(std,8);
        }
      else
        {
        int std = 0;
        str_std = ST::inttostring(std);
        }
      double min = Xi.min(0);
      double max = Xi.max(0);
      double med = Xi.quantile(50,0);

      int l_vn = vnames[i].length();
      if (l_vn > 9)
         {
         vnames[i] = vnames[i].substr(0, 9) + "~";
         }
      ST::string str_obs = ST::inttostring(obs);
      ST::string str_mean = ST::doubletostring(mean,8);
      ST::string str_med = ST::doubletostring(med,8);
      ST::string str_min = ST::doubletostring(min,8);
      ST::string str_max = ST::doubletostring(max,8);

      o.out(vnames[i].helpfill(11) + str_obs.helpfill(10) + str_mean.helpfill(15) +
            str_med.helpfill(15) + str_std.helpfill(15) + str_min.helpfill(15) +
            str_max.helpfill(15) + "\n");

      }
    o.out("\n");
    }

  }