search

easycure1 / vnpc.avg

0 stars 2 forks source link

Getting an error running `Rscript test_example.R` #2

Open avivajpeyi opened 1 week ago

avivajpeyi commented 1 week ago 
❯ Rscript test_example.R 

Attaching package: ‘vnpc.avg’

The following object is masked from ‘package:beyondWhittle’:

    psd_varma

Error in hasEigenValueSmallerZero(A, tol) : 
  could not find function "hasEigenValueSmallerZero"
Calls: <Anonymous> -> gibbs_m_avg_nuisance -> stopifnot -> is_hpd
In addition: Warning messages:
1: In vnpc.avg::gibbs_vnpc_avg(data = data, seg_n = 200, Ntotal = Ntotal,  :
  Data has been mean centered
2: In gibbs_m_avg_nuisance(data = data, mcmc_params = mcmc_params,  :
  Using conjugate sampling for theta. TODO: Get that consistent with MH!!
Execution halted
easycure1 commented 1 week ago 
❯ Rscript test_example.R 

Attaching package: ‘vnpc.avg’

The following object is masked from ‘package:beyondWhittle’:

    psd_varma

Error in hasEigenValueSmallerZero(A, tol) : 
  could not find function "hasEigenValueSmallerZero"
Calls: <Anonymous> -> gibbs_m_avg_nuisance -> stopifnot -> is_hpd
In addition: Warning messages:
1: In vnpc.avg::gibbs_vnpc_avg(data = data, seg_n = 200, Ntotal = Ntotal,  :
  Data has been mean centered
2: In gibbs_m_avg_nuisance(data = data, mcmc_params = mcmc_params,  :
  Using conjugate sampling for theta. TODO: Get that consistent with MH!!
Execution halted

Hi Avi, sorry for the late reply. I am having a fever... Please try devtools::document() before you install the package. One Rd file was not updated before.

easycure1 commented 1 week ago

In addition, if you only want to check fft, you can directly call vnpc.avg:::mdft()

avivajpeyi commented 1 week ago

oh no! hope you feel better!

Ill try your install suggestion and see if it works! thanks :)

avivajpeyi commented 1 week ago

Great! was able to get the example to work:

psd_matrix_plot

Added the plotting script for the example in #1

avivajpeyi commented 1 week ago

(we can close this issue now if youd' like :) )