ssvd: Sparse SVD
In ssvd: Sparse SVD

Description Usage Arguments Value Author(s) References See Also Examples

The function ssvd combines two functions ssvd.initial and ssvd.iter.thresh into one. Obtain sparse SVD using fast iterative thresholding method, together with a fast initialization algorithm

1
2
3

ssvd(x, method = c("theory", "method"), alpha.method = 0.05, alpha.theory = 1.5, 
huber.beta = 0.95, sigma = NA, r = 1, gamma.u = sqrt(2), gamma.v = sqrt(2), 
dothres = "hard", tol = 1e-08, n.iter = 100, n.boot = 100, non.orth = FALSE)

`x`	Input matrix, for which one would like to get a sparse SVD.
`method`	If method = "theory", then a theoretical procedure is adopted which is based on normal assumption on the noise. If method = "method", then the function bypass the normal assumption by some robust statistics. These two choices typically give similar solutions, but "theory" is much faster.
`alpha.method`	Alpha.method is the level of the hypothesis test when one performs Holm multiple hypothesis testing, which is used to select the candidate rows and columns in the initialization step. The value is usually set to be 0.05.
`alpha.theory`	Alpha.theory is a scaler that is used when normal assumption is true, method="theory", and a chisq tail bound is used to select the candidate rows and columns in the initialization step. Most of the time, users should keep it as it is.
`huber.beta`	Huber.beta is a scaler which is the cut-off point in the Huber function. The huberization is utilized to achieve robustness when normal assumption is violated in the initialization step.
`sigma`	Sigma is a scaler for the noise level. The user can set it to be NA, and the function will estimate it automatically.
`r`	A scaler, the number of components, i.e., the number of singular vectors to be computed.
`gamma.u`	When the method="theory", gamma.u=sqrt(2) corresponds to the sqrt(2 log(p)), which is the largest magnitude of p iid standard normals. If gamma.u is manually set to be smaller or larger than sqrt2, the left singular vectors will be denser or sparser respectively.
`gamma.v`	When the method="theory", gamma.u=sqrt(2) corresponds to the sqrt(2 log(p)), which is the largest magnitude of p iid standard normals. If gamma.u is manually set to be smaller or larger than sqrt2, the right singular vectors will be denser or sparser respectively.
`dothres`	Dothres has two choices, either "hard" or "soft", which means hard-thresholding or soft-thresholding.
`tol`	The tolerance level that determines when the algorithm stops.
`n.iter`	Maximum number of iterations allowed.
`n.boot`	Number of bootstrap to estimate the threshold level when method = "method"
`non.orth`	If non.orth=TRUE, then the last iteration of the algorithm will not involve orthoganolization, which should produce sparse solutions.

`u`	A matrix containing left singular vectors
`v`	A matrix containing right singular vectors
`d`	A vector containing singular values
`niter`	Number of iterations for the algorithm to converge
`sigma.hat`	An estimate of the noise level
`dist.u`	The distance between the left singular vectors of the last two iterations, can be used to see whether the algorithm indeed converges.
`dist.v`	The distance between the right singular vectors of the last two iterations, can be used to see whether the algorithm indeed converges.

Dan Yang

A Sparse SVD Method for High-dimensional Data

ssvd.initial and ssvd.iter.thresh

1	ssvd(matrix(rnorm(2^15),2^7,2^8), method = "method")

$u
       [,1]
  [1,]    0
  [2,]    0
  [3,]    0
  [4,]    0
  [5,]    0
  [6,]    0
  [7,]    0
  [8,]    0
  [9,]    0
 [10,]    0
 [11,]    0
 [12,]    0
 [13,]    0
 [14,]    0
 [15,]    0
 [16,]    0
 [17,]    0
 [18,]    0
 [19,]    0
 [20,]    0
 [21,]    0
 [22,]    0
 [23,]    0
 [24,]    0
 [25,]    0
 [26,]    0
 [27,]    0
 [28,]    0
 [29,]    0
 [30,]    0
 [31,]    0
 [32,]    0
 [33,]    0
 [34,]    0
 [35,]    0
 [36,]    0
 [37,]    0
 [38,]    0
 [39,]    0
 [40,]    0
 [41,]    0
 [42,]    0
 [43,]    0
 [44,]    0
 [45,]    0
 [46,]    0
 [47,]    0
 [48,]    0
 [49,]    0
 [50,]    0
 [51,]    0
 [52,]    0
 [53,]    0
 [54,]    0
 [55,]    0
 [56,]    0
 [57,]    0
 [58,]    0
 [59,]    0
 [60,]    0
 [61,]    0
 [62,]    0
 [63,]    0
 [64,]    0
 [65,]    0
 [66,]    1
 [67,]    0
 [68,]    0
 [69,]    0
 [70,]    0
 [71,]    0
 [72,]    0
 [73,]    0
 [74,]    0
 [75,]    0
 [76,]    0
 [77,]    0
 [78,]    0
 [79,]    0
 [80,]    0
 [81,]    0
 [82,]    0
 [83,]    0
 [84,]    0
 [85,]    0
 [86,]    0
 [87,]    0
 [88,]    0
 [89,]    0
 [90,]    0
 [91,]    0
 [92,]    0
 [93,]    0
 [94,]    0
 [95,]    0
 [96,]    0
 [97,]    0
 [98,]    0
 [99,]    0
[100,]    0
[101,]    0
[102,]    0
[103,]    0
[104,]    0
[105,]    0
[106,]    0
[107,]    0
[108,]    0
[109,]    0
[110,]    0
[111,]    0
[112,]    0
[113,]    0
[114,]    0
[115,]    0
[116,]    0
[117,]    0
[118,]    0
[119,]    0
[120,]    0
[121,]    0
[122,]    0
[123,]    0
[124,]    0
[125,]    0
[126,]    0
[127,]    0
[128,]    0

$v
       [,1]
  [1,]    0
  [2,]    0
  [3,]    0
  [4,]    0
  [5,]    0
  [6,]    0
  [7,]    0
  [8,]    0
  [9,]   -1
 [10,]    0
 [11,]    0
 [12,]    0
 [13,]    0
 [14,]    0
 [15,]    0
 [16,]    0
 [17,]    0
 [18,]    0
 [19,]    0
 [20,]    0
 [21,]    0
 [22,]    0
 [23,]    0
 [24,]    0
 [25,]    0
 [26,]    0
 [27,]    0
 [28,]    0
 [29,]    0
 [30,]    0
 [31,]    0
 [32,]    0
 [33,]    0
 [34,]    0
 [35,]    0
 [36,]    0
 [37,]    0
 [38,]    0
 [39,]    0
 [40,]    0
 [41,]    0
 [42,]    0
 [43,]    0
 [44,]    0
 [45,]    0
 [46,]    0
 [47,]    0
 [48,]    0
 [49,]    0
 [50,]    0
 [51,]    0
 [52,]    0
 [53,]    0
 [54,]    0
 [55,]    0
 [56,]    0
 [57,]    0
 [58,]    0
 [59,]    0
 [60,]    0
 [61,]    0
 [62,]    0
 [63,]    0
 [64,]    0
 [65,]    0
 [66,]    0
 [67,]    0
 [68,]    0
 [69,]    0
 [70,]    0
 [71,]    0
 [72,]    0
 [73,]    0
 [74,]    0
 [75,]    0
 [76,]    0
 [77,]    0
 [78,]    0
 [79,]    0
 [80,]    0
 [81,]    0
 [82,]    0
 [83,]    0
 [84,]    0
 [85,]    0
 [86,]    0
 [87,]    0
 [88,]    0
 [89,]    0
 [90,]    0
 [91,]    0
 [92,]    0
 [93,]    0
 [94,]    0
 [95,]    0
 [96,]    0
 [97,]    0
 [98,]    0
 [99,]    0
[100,]    0
[101,]    0
[102,]    0
[103,]    0
[104,]    0
[105,]    0
[106,]    0
[107,]    0
[108,]    0
[109,]    0
[110,]    0
[111,]    0
[112,]    0
[113,]    0
[114,]    0
[115,]    0
[116,]    0
[117,]    0
[118,]    0
[119,]    0
[120,]    0
[121,]    0
[122,]    0
[123,]    0
[124,]    0
[125,]    0
[126,]    0
[127,]    0
[128,]    0
[129,]    0
[130,]    0
[131,]    0
[132,]    0
[133,]    0
[134,]    0
[135,]    0
[136,]    0
[137,]    0
[138,]    0
[139,]    0
[140,]    0
[141,]    0
[142,]    0
[143,]    0
[144,]    0
[145,]    0
[146,]    0
[147,]    0
[148,]    0
[149,]    0
[150,]    0
[151,]    0
[152,]    0
[153,]    0
[154,]    0
[155,]    0
[156,]    0
[157,]    0
[158,]    0
[159,]    0
[160,]    0
[161,]    0
[162,]    0
[163,]    0
[164,]    0
[165,]    0
[166,]    0
[167,]    0
[168,]    0
[169,]    0
[170,]    0
[171,]    0
[172,]    0
[173,]    0
[174,]    0
[175,]    0
[176,]    0
[177,]    0
[178,]    0
[179,]    0
[180,]    0
[181,]    0
[182,]    0
[183,]    0
[184,]    0
[185,]    0
[186,]    0
[187,]    0
[188,]    0
[189,]    0
[190,]    0
[191,]    0
[192,]    0
[193,]    0
[194,]    0
[195,]    0
[196,]    0
[197,]    0
[198,]    0
[199,]    0
[200,]    0
[201,]    0
[202,]    0
[203,]    0
[204,]    0
[205,]    0
[206,]    0
[207,]    0
[208,]    0
[209,]    0
[210,]    0
[211,]    0
[212,]    0
[213,]    0
[214,]    0
[215,]    0
[216,]    0
[217,]    0
[218,]    0
[219,]    0
[220,]    0
[221,]    0
[222,]    0
[223,]    0
[224,]    0
[225,]    0
[226,]    0
[227,]    0
[228,]    0
[229,]    0
[230,]    0
[231,]    0
[232,]    0
[233,]    0
[234,]    0
[235,]    0
[236,]    0
[237,]    0
[238,]    0
[239,]    0
[240,]    0
[241,]    0
[242,]    0
[243,]    0
[244,]    0
[245,]    0
[246,]    0
[247,]    0
[248,]    0
[249,]    0
[250,]    0
[251,]    0
[252,]    0
[253,]    0
[254,]    0
[255,]    0
[256,]    0

$d
[1] 2.990177

$niter
[1] 1

$sigma.hat
[1] 0.9985177

$dist.u
[1] 0

$dist.v
[1] 0

Warning message:
In ssvd.iter.thresh(x, method = method, u.old = ans.initial$u, v.old = ans.initial$v,  :
  SSVD.iter.thresh: Overthresh: all zero!