Advertisement
Guest User

Untitled

a guest
Aug 20th, 2019
108
0
Never
Not a member of Pastebin yet? Sign Up, it unlocks many cool features!
text 1.27 KB | None | 0 0
  1. var1,var2,var3
  2. 0.942856823,0.568425866,0.325885379
  3. 1.227681099,1.335672206,0.925331054
  4. 1.952671045,1.829479996,1.512280854
  5. 2.45428731,1.990174152,1.534456808
  6. 2.987783477,2.78975186,1.725095748
  7. 3.651682331,2.966399127,1.972274564
  8. 3.768010479,3.211381506,1.993080807
  9. 4.509429614,3.642983433,2.541071547
  10. 4.81498729,3.888415006,3.218031802
  11.  
  12. open System.IO
  13. open MathNet.Numerics.LinearAlgebra
  14. let rows = [|for line in File.ReadAllLines("Z:\mypath.csv")
  15. |> Seq.skip 1 do yield line.Split(',') |> Array.map float|]
  16. let data = DenseMatrix.ofRowArrays rows
  17.  
  18. let data_logdiff =
  19. DenseMatrix.init (data.RowCount-1) (data.ColumnCount)
  20. (fun j i -> if j = 0 then 0. else data.At(j, i) / data.At(j-1, i) |> log)
  21.  
  22. let alpha = vector [for i in data_logdiff.EnumerateColumns() -> i |> Statistics.Mean]
  23.  
  24. let sigsq (values:Vector<float>) (avg: float) =
  25. let sqr x = x * x
  26. let result = values |> (fun i -> sqr (i - avg))
  27. result
  28.  
  29. sigsq (data_logdiff.Column(i), alpha.[0]) |> printfn "%A"
  30.  
  31. let sigsq (values:Vector<float>) (avg: float) =
  32.  
  33. sigsq (data_logdiff.Column(i), alpha.[0]) |> printfn "%A"
  34.  
  35. sigsq (data_logdiff.Column(0)) (alpha.[0])
  36.  
  37. sigsq <| data_logdiff.Column(0) <| alpha.[0]
  38.  
  39. data_logdiff.Column(0) |> sigsq <| alpha.[0]
Advertisement
Add Comment
Please, Sign In to add comment
Advertisement