Bootstrap a pairwise distance matrix

Hi! I have to bootstrap a distance matrix to have the reliability of my phylogenetic tree. I found bootstrp as matlab function but I am not certain about what to do

%e.g. let us suppose that I have this matrix

M=magic(6);
% Just to have an idea
[bootstat,bootsam]=bootstrp(1000,@fun,M);
%I don't know what kind of fucntion can be correct 

I there anyone with some experience in boot strapping a distancematrix for phylogeny?

Thank you in advance

[Adam Danz]

Bootstrapping involves taking many randomized subsamples of your data, performing the same statistic or function on each subsample, and then analyzing the distribution of results. For example, if you want confidence intervals on the median value of some distribution, you should take randomized subsamples of the population, compute and store the median of the subpopulation, and repeat the process many times until you have a distribution of median values. In that example, @fun would be @median.

Let's say your distance matrix is computed by your homemade function,

output = computePhyloDistance(input)

then your @fun would be @computePhyloDistance.

Step through some of the examples provided by Matlab and the process will become clear.

https://www.mathworks.com/help/stats/bootstrp.html#brx09qf-190

Answer 1

[Shubh Sahu]

Hi,

Adding to Adam Danz, Bootstrapping is a random sampling conducted with replacement. It is straight forward way to derive estimates of standard error and confidence intervals for complex estimators of complex parameters of the distribution.

Here is the link you can refer to:

https://www.mathworks.com/help/stats/bootstrp.html#brx09qf-190.