Skip to content

bug in calculation of R2s #47

@JasperVanDenBosch

Description

@JasperVanDenBosch
~/Projects/pyGLMdenoise/glmdenoise/pyGlmdenoise.py in fit(self)
    385         self.results['R2s'] = calccodStack(
    386             whitened_data,
--> 387             modelfits)
    388         """ TO DO
    389         self.results['R2runs'] = [calccod(whitened_data[c_run], modelfits[c_run], 0)

~/Projects/pyGLMdenoise/glmdenoise/utils/optimiseHRF.py in calccodStack(y, yhat)
    198         with np.errstate(divide="ignore", invalid="ignore"):
    199             yrun = np.array(y[run])
--> 200             yhatrun = np.array(yhat[run])
    201             nom.append(np.sum((yrun - yhatrun) ** 2, axis=0))
    202             denom.append(np.sum(yrun ** 2, axis=0))  # Kendricks denominator

~/Projects/pyGLMdenoise/env/lib/python3.6/site-packages/numpy/matrixlib/defmatrix.py in __getitem__(self, index)
    193 
    194         try:
--> 195             out = N.ndarray.__getitem__(self, index)
    196         finally:
    197             self._getitem = False

IndexError: index 8 is out of bounds for axis 0 with size 8
> /home/adf/vandejjf/Projects/pyGLMdenoise/env/lib/python3.6/site-packages/numpy/matrixlib/defmatrix.py(195)__getitem__()
    193 
    194         try:
--> 195             out = N.ndarray.__getitem__(self, index)
    196         finally:
    197             self._getitem = False

ipdb> y   

Metadata

Metadata

Assignees

No one assigned

    Labels

    No labels
    No labels

    Type

    No type

    Projects

    No projects

    Milestone

    No milestone

    Relationships

    None yet

    Development

    No branches or pull requests

    Issue actions