EVOLUTION-MANAGER
Edit File: dim_select.html
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Strict//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-strict.dtd"><html xmlns="http://www.w3.org/1999/xhtml"><head><title>R: Dimensionality selection for singular values using profile...</title> <meta http-equiv="Content-Type" content="text/html; charset=utf-8" /> <link rel="stylesheet" type="text/css" href="R.css" /> </head><body> <table width="100%" summary="page for dim_select {igraph}"><tr><td>dim_select {igraph}</td><td style="text-align: right;">R Documentation</td></tr></table> <h2>Dimensionality selection for singular values using profile likelihood.</h2> <h3>Description</h3> <p>Select the number of significant singular values, by finding the ‘elbow’ of the scree plot, in a principled way. </p> <h3>Usage</h3> <pre> dim_select(sv) </pre> <h3>Arguments</h3> <table summary="R argblock"> <tr valign="top"><td><code>sv</code></td> <td> <p>A numeric vector, the ordered singular values.</p> </td></tr> </table> <h3>Details</h3> <p>The input of the function is a numeric vector which contains the measure of ‘importance’ for each dimension. </p> <p>For spectral embedding, these are the singular values of the adjacency matrix. The singular values are assumed to be generated from a Gaussian mixture distribution with two components that have different means and same variance. The dimensionality <i>d</i> is chosen to maximize the likelihood when the <i>d</i> largest singular values are assigned to one component of the mixture and the rest of the singular values assigned to the other component. </p> <p>This function can also be used for the general separation problem, where we assume that the left and the right of the vector are coming from two Normal distributions, with different means, and we want to know their border. See examples below. </p> <h3>Value</h3> <p>A numeric scalar, the estimate of <i>d</i>. </p> <h3>Author(s)</h3> <p>Gabor Csardi <a href="mailto:csardi.gabor@gmail.com">csardi.gabor@gmail.com</a> </p> <h3>References</h3> <p>M. Zhu, and A. Ghodsi (2006). Automatic dimensionality selection from the scree plot via the use of profile likelihood. <em>Computational Statistics and Data Analysis</em>, Vol. 51, 918–930. </p> <h3>See Also</h3> <p><code><a href="embed_adjacency_matrix.html">embed_adjacency_matrix</a></code> </p> <h3>Examples</h3> <pre> # Generate the two groups of singular values with # Gaussian mixture of two components that have different means sing.vals <- c( rnorm (10, mean=1, sd=1), rnorm(10, mean=3, sd=1) ) dim.chosen <- dim_select(sing.vals) dim.chosen # Sample random vectors with multivariate normal distribution # and normalize to unit length lpvs <- matrix(rnorm(200), 10, 20) lpvs <- apply(lpvs, 2, function(x) { (abs(x) / sqrt(sum(x^2))) }) RDP.graph <- sample_dot_product(lpvs) dim_select( embed_adjacency_matrix(RDP.graph, 10)$D ) # Sample random vectors with the Dirichlet distribution lpvs.dir <- sample_dirichlet(n=20, rep(1, 10)) RDP.graph.2 <- sample_dot_product(lpvs.dir) dim_select( embed_adjacency_matrix(RDP.graph.2, 10)$D ) # Sample random vectors from hypersphere with radius 1. lpvs.sph <- sample_sphere_surface(dim=10, n=20, radius=1) RDP.graph.3 <- sample_dot_product(lpvs.sph) dim_select( embed_adjacency_matrix(RDP.graph.3, 10)$D ) </pre> <hr /><div style="text-align: center;">[Package <em>igraph</em> version 1.3.5 <a href="00Index.html">Index</a>]</div> </body></html>