sdp: Semi-Definite Programming Algorithm

Description Usage Arguments Details Value See Also Examples

View source: R/sdp.R

Description

sdp returns a completed Euclidean Distance Matrix D, with dimension d, from a partial Euclidean Distance Matrix using the methods of Alfakih et. al. (1999)

Usage

1
sdp(D, A, toler = 1e-08)

Arguments

D

An nxn partial-distance matrix to be completed. D must satisfy a list of conditions (see details), with unkown entries set to NA.

A

a weight matrix, with h_{ij} = 0 implying a_{ij} is unknown. Generally, if a_{ij} is known, h_{ij} = 1, although any non-negative weight is allowed.

toler

convergence tolerance for the algorithm

Details

This is an implementation of the Semi-Definite Programming Algorithm (sdp) for Euclidean Distance Matrix Completion, as proposed in 'Solving Euclidean Distance Matrix Completion Problems via Semidefinite Programming' (Alfakih et. al., 1999).

The method seeks to minimize the following:

||A \cdot (D - psd2edm(S))||_{F}^{2}

where the function psd2edm() is that described in psd2edm(), and the norm is Frobenius. Minimization is over S, a positive semidefinite matrix.

The matrix D is a partial-distance matrix, meaning some of its entries are unknown. It must satisfy the following conditions in order to be completed:

Value

D

an nxn matrix of the completed Euclidean distances

optval

the minimum value achieved of the target function during minimization

See Also

psd2edm

Examples

 1
 2
 3
 4
 5
 6
 7
 8
 9
10
11
12
13
14
D <- matrix(c(0,3,4,3,4,3,
             3,0,1,NA,5,NA,
             4,1,0,5,NA,5,
             3,NA,5,0,1,NA,
             4,5,NA,1,0,5,
             3,NA,5,NA,5,0), byrow=TRUE, nrow=6)
A <- matrix(c(1,1,1,1,1,1,
             1,1,1,0,1,0,
             1,1,1,1,0,1,
             1,0,1,1,1,0,
             1,1,0,1,1,1,
             1,0,1,0,1,1), byrow=TRUE, nrow=6)
             
edmc(D, method="sdp", A=A, toler=1e-2)

great-northern-diver/edmcr documentation built on Dec. 20, 2021, 12:52 p.m.