projects
/
valse.git
/ blobdiff
commit
grep
author
committer
pickaxe
?
search:
re
summary
|
shortlog
|
log
|
commit
|
commitdiff
|
tree
raw
|
inline
| side by side
fix initialization and made some update
[valse.git]
/
pkg
/
R
/
computeGridLambda.R
diff --git
a/pkg/R/computeGridLambda.R
b/pkg/R/computeGridLambda.R
index
e2b6303
..
9d06aed
100644
(file)
--- a/
pkg/R/computeGridLambda.R
+++ b/
pkg/R/computeGridLambda.R
@@
-17,7
+17,7
@@
#'
#' @export
computeGridLambda = function(phiInit, rhoInit, piInit, gamInit, X, Y,
#'
#' @export
computeGridLambda = function(phiInit, rhoInit, piInit, gamInit, X, Y,
- gamma, mini, maxi, tau)
+ gamma, mini, maxi, tau
, fast=TRUE
)
{
n = nrow(X)
p = dim(phiInit)[1]
{
n = nrow(X)
p = dim(phiInit)[1]
@@
-25,7
+25,7
@@
computeGridLambda = function(phiInit, rhoInit, piInit, gamInit, X, Y,
k = dim(phiInit)[3]
list_EMG = EMGLLF(phiInit, rhoInit, piInit, gamInit, mini, maxi,
k = dim(phiInit)[3]
list_EMG = EMGLLF(phiInit, rhoInit, piInit, gamInit, mini, maxi,
- gamma, lambda=0, X, Y, tau)
+ gamma, lambda=0, X, Y, tau
, fast
)
grid = array(0, dim=c(p,m,k))
for (i in 1:p)
{
grid = array(0, dim=c(p,m,k))
for (i in 1:p)
{
@@
-33,5
+33,5
@@
computeGridLambda = function(phiInit, rhoInit, piInit, gamInit, X, Y,
grid[i,j,] = abs(list_EMG$S[i,j,]) / (n*list_EMG$pi^gamma)
}
grid = unique(grid)
grid[i,j,] = abs(list_EMG$S[i,j,]) / (n*list_EMG$pi^gamma)
}
grid = unique(grid)
- grid
+ sort(grid)
}
}