/* * CalibratedYuleInitialTree.java * * Copyright (C) 2002-2006 Alexei Drummond and Andrew Rambaut * * This file is part of BEAST. * See the NOTICE file distributed with this work for additional * information regarding copyright ownership and licensing. * * BEAST is free software; you can redistribute it and/or modify * it under the terms of the GNU Lesser General Public License as * published by the Free Software Foundation; either version 2 * of the License, or (at your option) any later version. * * BEAST is distributed in the hope that it will be useful, * but WITHOUT ANY WARRANTY; without even the implied warranty of * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the * GNU Lesser General Public License for more details. * * You should have received a copy of the GNU Lesser General Public * License along with BEAST; if not, write to the * Free Software Foundation, Inc., 51 Franklin St, Fifth Floor, * Boston, MA 02110-1301 USA */ package beast.evolution.speciation; import java.util.ArrayList; import java.util.List; import org.apache.commons.math.MathException; import beast.core.Description; import beast.core.Input; import beast.core.StateNode; import beast.core.StateNodeInitialiser; import beast.evolution.tree.Tree; /** * @author Joseph Heled */ @Description("A tree compatible with a set of monophyletic clades and hard limits on the clade root.") public class CalibratedYuleInitialTree extends Tree implements StateNodeInitialiser { final public Input<List<CalibrationPoint>> calibrationsInput = new Input<>("calibrations", "Set of calibrated nodes", new ArrayList<>(), Input.Validate.REQUIRED); @Override public void initAndValidate() { super.initAndValidate(); initStateNodes(); } @Override public void initStateNodes() { // Would have been nice to use the MCMC CalibratedYuleModel beastObject directly, but at this point // it does not exist since the tree being initialized is one of its arguments. So, build a temporary // one using the initializer tree. final List<CalibrationPoint> cals = calibrationsInput.get(); final CalibratedYuleModel cym = new CalibratedYuleModel(); for( final CalibrationPoint cal : cals ) { cym.setInputValue("calibrations", cal); } cym.setInputValue("tree", this); cym.setInputValue("type", CalibratedYuleModel.Type.NONE); cym.initAndValidate(); Tree t; try { t = cym.compatibleInitialTree(); } catch (MathException e) { throw new IllegalArgumentException(e.getMessage()); } m_initial.get().assignFromWithoutID(t); } @Override public void getInitialisedStateNodes(final List<StateNode> stateNodes) { stateNodes.add(m_initial.get()); } }