function st = prepareDoubleExpResult(st,cfun,gof) if ~isempty(cfun) cf1 = confint(cfun); st.tau1 = cfun.tau1; st.tau1_confint_low = cf1(1,1); st.tau1_confint_high= cf1(2,1); st.tau2 = cfun.tau2; st.tau2_confint_low = cf1(1,2); st.tau2_confint_high= cf1(2,2); st.A1 = cfun.A1; st.A1_confint_low = cf1(1,3); st.A1_confint_high = cf1(2,3); st.A2 = cfun.A2; st.A2_confint_low = cf1(1,4); st.A2_confint_high = cf1(2,4); st.Rsquare_double = gof.rsquare; else st.tau1 = NaN; st.tau1_confint_low = NaN; st.tau1_confint_high = NaN; st.tau2 = NaN; st.tau2_confint_low = NaN; st.tau2_confint_high = NaN; st.A1 = NaN; st.A1_confint_low = NaN; st.A1_confint_high = NaN; st.A2 = NaN; st.A2_confint_low = NaN; st.A2_confint_high = NaN; st.Rsquare_double = NaN; end end