diff options
| author | Suraj Yerramilli | 2016-02-14 11:35:52 +0530 |
|---|---|---|
| committer | Suraj Yerramilli | 2016-02-14 11:35:52 +0530 |
| commit | f0242cac8f81ad471069b0e40226cea585e85a73 (patch) | |
| tree | 6fb74dfb9e9b3e45580551cf9f88524d29014f59 /R | |
| parent | b2ced3712b06118a15d32c586bbc1deeac2c304f (diff) | |
| download | SysID-R-code-f0242cac8f81ad471069b0e40226cea585e85a73.tar.gz SysID-R-code-f0242cac8f81ad471069b0e40226cea585e85a73.tar.bz2 SysID-R-code-f0242cac8f81ad471069b0e40226cea585e85a73.zip | |
minor bug fix
Diffstat (limited to 'R')
| -rw-r--r-- | R/estUtil.R | 3 |
1 files changed, 1 insertions, 2 deletions
diff --git a/R/estUtil.R b/R/estUtil.R index 2d1a76f..474aa5d 100644 --- a/R/estUtil.R +++ b/R/estUtil.R @@ -13,7 +13,6 @@ levbmqdt <- function(...,obj,theta0,N,opt){ l <- obj(theta=theta0,e=NULL,dots) e <- l$Y-l$X%*%theta0 sumsq0 <- sum(e^2)/N - theta <- theta0 # variable to count the number of times objective function is called countObj <- 0 @@ -70,7 +69,7 @@ levbmqdt <- function(...,obj,theta0,N,opt){ break } } - + theta <- theta0 e <- e[1:N,] sigma2 <- sum(e^2)/df vcov <- 1/N*Hinv*sigma2 |