求一組數(shù)的偏差和累加

q1208c

回復(fù) 8# yestreenstars

連最小二乘法都上來(lái)了, 你覺得shell能做這么復(fù)雜的計(jì)算么?!

我確實(shí)是來(lái)抬杠的, 因?yàn)槲矣X得shell的專長(zhǎng)并不在此.

殺雞雖然不需要用牛刀, 但殺牛, 用雞刀可是個(gè)體力活.
   

Herowinter

@wpd0508
先是平均數(shù)，然后是方差，
接下來(lái)是最小2乘法和曲線擬合，
恭喜星辰大大入坑。

Herowinter

@yestreenstars@q1208c
先是平均數(shù)，然后是方差，
接下來(lái)是最小2乘法和曲線擬合，
恭喜星辰大大入坑。

yestreenstars

回復(fù) 12# Herowinter
咱數(shù)學(xué)沒學(xué)得那么深入~{:2_169:}

   

用戶名注冊(cè)后不能更改

回復(fù) 12# Herowinter

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

ly5066113

貌似是這個(gè)意思：

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

復(fù)制代碼

wpd0508

看著很靠譜，能順便加上最后一步 E = sqrt( sum (Yi - (a1*xi +a0) )^2)的步驟嗎回復(fù) 16# ly5066113


   

wpd0508

我剛試了一下大大的代碼，為什么輸出的是y=-nan+-nanx回復(fù) 16# ly5066113



   

damcool

這一點(diǎn)點(diǎn)都不會(huì)改，那就沒人能幫你了！

wpd0508 發(fā)表于 2014-06-24 14:45
我剛試了一下大大的代碼，為什么輸出的是y=-nan+-nanx回復(fù) 16# ly5066113

wpd0508

提示一下唄，
回復(fù) 19# damcool