求一组数的偏差和累加

q1208c

回复 8# yestreenstars

连最小二乘法都上来了, 你觉得shell能做这么复杂的计算么?!

我确实是来抬杠的, 因为我觉得shell的专长并不在此.

杀鸡虽然不需要用牛刀, 但杀牛, 用鸡刀可是个体力活.
   

Herowinter

@wpd0508
先是平均数，然后是方差，
接下来是最小2乘法和曲线拟合，
恭喜星辰大大入坑。

Herowinter

@yestreenstars@q1208c
先是平均数，然后是方差，
接下来是最小2乘法和曲线拟合，
恭喜星辰大大入坑。

yestreenstars

回复 12# Herowinter
咱数学没学得那么深入~{:2_169:}

   

用户名注册后不能更改

回复 12# Herowinter

我知道最小是minimum，二是silly，乘法是multiplication，但最小二乘法是虾米意思？

ly5066113

貌似是这个意思：

1. $ cat file
   
2. 1
   
3. 3
   
4. 3
   
5. 5
   
6. 4
   
7. 2
   
8. $ awk '{
   
9.         a[NR] = $1
   
10.         s += $1
    
11. }
    
12. END{
    
13.         n = NR
    
14.         avg = s / n
    
15.         for(x=1; x<=n; x++)
    
16.         {
    
17.                 y = a[x] - avg
    
18.                 sx += x
    
19.                 sy += y
    
20.                 sxy += x * y
    
21.                 sx2 += x * x
    
22.         }
    
23.         a0 = (sx2 * sy - sx * sxy) / (n * sx2 - sx * sx)
    
24.         a1 = (n * sxy - sx * sy) / (n * sx2 - sx * sx)
    
25.         printf "y = %s + %sx\n" ,a0 ,a1
    
26. }' file
    
27. 
    
28. y = -1 + 0.285714x

复制代码

wpd0508

看着很靠谱，能顺便加上最后一步 E = sqrt( sum (Yi - (a1*xi +a0) )^2)的步骤吗回复 16# ly5066113


   

wpd0508

我刚试了一下大大的代码，为什么输出的是y=-nan+-nanx回复 16# ly5066113



   

damcool

这一点点都不会改，那就没人能帮你了！

wpd0508 发表于 2014-06-24 14:45
我刚试了一下大大的代码，为什么输出的是y=-nan+-nanx回复 16# ly5066113

wpd0508

提示一下呗，
回复 19# damcool