Hi. I would like to fit data to a grafted normal-weibull distribution.
The red line is a normal distribution and the blue line is a weibull distribution.
The grafting point between the normal and weibulll distribution must share the same value and slope on both sides. This gives two equations. There are four parameters (a and b for the weibull distribution, and alpha and sigma for the normal distribution).
My current approach is to use fminunc to find values of mu and sigma which minizime the square error. The values of a and b are solved for with vpasolve using the two constraints, and the error is computed.
This works fairly well (it was how I created the figure); however it often hits an error when vpasolve is unable to find values for a and b.
Is there a better approach?