ExponentialRegression returns an exponential regression on data series
(s Series)
| 47 | |
| 48 | // ExponentialRegression returns an exponential regression on data series |
| 49 | func ExponentialRegression(s Series) (regressions Series, err error) { |
| 50 | |
| 51 | if len(s) == 0 { |
| 52 | return nil, EmptyInputErr |
| 53 | } |
| 54 | |
| 55 | var sum [6]float64 |
| 56 | |
| 57 | for i := 0; i < len(s); i++ { |
| 58 | if s[i].Y < 0 { |
| 59 | return nil, YCoordErr |
| 60 | } |
| 61 | sum[0] += s[i].X |
| 62 | sum[1] += s[i].Y |
| 63 | sum[2] += s[i].X * s[i].X * s[i].Y |
| 64 | sum[3] += s[i].Y * math.Log(s[i].Y) |
| 65 | sum[4] += s[i].X * s[i].Y * math.Log(s[i].Y) |
| 66 | sum[5] += s[i].X * s[i].Y |
| 67 | } |
| 68 | |
| 69 | denominator := (sum[1]*sum[2] - sum[5]*sum[5]) |
| 70 | a := math.Pow(math.E, (sum[2]*sum[3]-sum[5]*sum[4])/denominator) |
| 71 | b := (sum[1]*sum[4] - sum[5]*sum[3]) / denominator |
| 72 | |
| 73 | for j := 0; j < len(s); j++ { |
| 74 | regressions = append(regressions, Coordinate{ |
| 75 | X: s[j].X, |
| 76 | Y: a * math.Exp(b*s[j].X), |
| 77 | }) |
| 78 | } |
| 79 | |
| 80 | return regressions, nil |
| 81 | } |
| 82 | |
| 83 | // LogarithmicRegression returns an logarithmic regression on data series |
| 84 | func LogarithmicRegression(s Series) (regressions Series, err error) { |
no outgoing calls
searching dependent graphs…