hi,

I ran a bootstrap using the `Maxboot` function (from the Maximum Likliehood application package) and save the parameters into a file let say bootparameter.

My question is , how can i see the parameters distribution or read that bootparameter file?

If I try to read it with the `loadd` function, I get the error:

error G0125:Read past end of file

## 4 Answers

0

You are wanting to view the bootstrapped parameter estimates that the `maxboot` function saves to a GAUSS dataset, correct? If so, they should be in a file with a `.dat` extension, such as `bootparameter.dat`. You can load the contents of a GAUSS dataset into a file with the `loadd` function like this:

```
p_boot = loadd("bootparameter");
```

If `bootparameter.dat` is not in your GAUSS current working directory, make sure to add the path, i.e.:

```
p_boot = loadd("C:\\gauss13\\myproj\\bootparameter");
```

0

You can produce a plot showing the distribution using the MaxDensity procedure that comes with Maxlik. The arguments are the dataset containing the bootstrapped coefficients, and a list of columns that you want to be displayed.

**> MAXdensity

**

** Purpose: To compute kernel density estimate and plot.

**

**

** Format: ( px,py,sm } = MAXdensity(dataset,pars);

**

** Input: dataset string, name of GAUSS dataset

** containing data.

**

** pars Kx1 vector, selected columns for

** estimation and display.

**

**

** Output: px _max_NumPointsx1 vector, abscissae.

** py _max_NumPointsx1 vector, ordinates.

** sm Kx1, or Nxk, or Nx1 smoothing coefficients.

**

** Remarks:

**

** kernel density plots of the selected parameters are

** generated.

**

** Globals:

**

** _max_Kernel Kx1 character vector, type of kernel:

**

** NORMAL - normal kernel

** EPAN - Epanechnikov kernel

** BIWGT - biweight kernel

** TRIANG - triangular kernel

** RECTANG - rectangular kernel

** TNORMAL - truncated normal kernel

**

** If _max_Kernel is scalar, the kernel is the same

** for all parameters. Default = { NORMAL };

**

** _max_NumPoints scalar, number of points to be computed for plots

**

** _max_EndPoints Kx2 matrix, lower (in first column) and upper

** (in second column) endpoints of density. Default is

** minimum and maximum, respectively, of the parameter

** values. If 1x2 matrix, endpoints will be the same

** for all parameters.

**

** _max_Smoothing Kx1 vector or Nx1 vector or NxK matrix, smoothing

** coefficients for each plot. If scalar, smoothing

** coefficient will be the same for each plot. If zero,

** smoothing coefficient will be computed by CMLdensity.

** If matrix, smoothing coefficient will be different for

** each observation.

** Default = 0;

**

** _max_Truncate Kx2 matrix, lower (in first column) and upper (in

** second column) truncation limits for truncated normal

** kernel. If 1x2 matrix, truncations limits will be the

** same for all plots. Default is minimum and maximum,

** respectively.

0

thank all of you. I have tried your suggestions but that the error message i keep receiving:

C:\gauss12\src\saveload.src(54) : error G0125 : Read past end of file

Currently active call: loadd [54] C:\gauss12\src\saveload.src

Stack trace:

loadd

0

I opened the dataset file that you sent in and checked the dimensions of the dataset like this:

open fh = eltcoeff for read; r = rowsf(fh); c = colsf(fh);

It reported that the dataset had 0 rows and 36 columns. I next checked to see what variables were in the dataset like this:

//using same 'fh' file handle from above vnames = getnamef(fh); print vnames;

and saw that the variable names:

PAR_1 PAR_2 PAR_3 . . . PAR_36

had been added to the dataset, but no values were assigned to them. It seems that a problem is occurring such that the data is not being written to the dataset.

## Your Answer

## 4 Answers

You are wanting to view the bootstrapped parameter estimates that the `maxboot` function saves to a GAUSS dataset, correct? If so, they should be in a file with a `.dat` extension, such as `bootparameter.dat`. You can load the contents of a GAUSS dataset into a file with the `loadd` function like this:

```
p_boot = loadd("bootparameter");
```

If `bootparameter.dat` is not in your GAUSS current working directory, make sure to add the path, i.e.:

```
p_boot = loadd("C:\\gauss13\\myproj\\bootparameter");
```

You can produce a plot showing the distribution using the MaxDensity procedure that comes with Maxlik. The arguments are the dataset containing the bootstrapped coefficients, and a list of columns that you want to be displayed.

**> MAXdensity

**

** Purpose: To compute kernel density estimate and plot.

**

**

** Format: ( px,py,sm } = MAXdensity(dataset,pars);

**

** Input: dataset string, name of GAUSS dataset

** containing data.

**

** pars Kx1 vector, selected columns for

** estimation and display.

**

**

** Output: px _max_NumPointsx1 vector, abscissae.

** py _max_NumPointsx1 vector, ordinates.

** sm Kx1, or Nxk, or Nx1 smoothing coefficients.

**

** Remarks:

**

** kernel density plots of the selected parameters are

** generated.

**

** Globals:

**

** _max_Kernel Kx1 character vector, type of kernel:

**

** NORMAL - normal kernel

** EPAN - Epanechnikov kernel

** BIWGT - biweight kernel

** TRIANG - triangular kernel

** RECTANG - rectangular kernel

** TNORMAL - truncated normal kernel

**

** If _max_Kernel is scalar, the kernel is the same

** for all parameters. Default = { NORMAL };

**

** _max_NumPoints scalar, number of points to be computed for plots

**

** _max_EndPoints Kx2 matrix, lower (in first column) and upper

** (in second column) endpoints of density. Default is

** minimum and maximum, respectively, of the parameter

** values. If 1x2 matrix, endpoints will be the same

** for all parameters.

**

** _max_Smoothing Kx1 vector or Nx1 vector or NxK matrix, smoothing

** coefficients for each plot. If scalar, smoothing

** coefficient will be the same for each plot. If zero,

** smoothing coefficient will be computed by CMLdensity.

** If matrix, smoothing coefficient will be different for

** each observation.

** Default = 0;

**

** _max_Truncate Kx2 matrix, lower (in first column) and upper (in

** second column) truncation limits for truncated normal

** kernel. If 1x2 matrix, truncations limits will be the

** same for all plots. Default is minimum and maximum,

** respectively.

thank all of you. I have tried your suggestions but that the error message i keep receiving:

C:\gauss12\src\saveload.src(54) : error G0125 : Read past end of file

Currently active call: loadd [54] C:\gauss12\src\saveload.src

Stack trace:

loadd

I opened the dataset file that you sent in and checked the dimensions of the dataset like this:

open fh = eltcoeff for read; r = rowsf(fh); c = colsf(fh);

It reported that the dataset had 0 rows and 36 columns. I next checked to see what variables were in the dataset like this:

//using same 'fh' file handle from above vnames = getnamef(fh); print vnames;

and saw that the variable names:

PAR_1 PAR_2 PAR_3 . . . PAR_36

had been added to the dataset, but no values were assigned to them. It seems that a problem is occurring such that the data is not being written to the dataset.