fork download
copy 
  1. regressionLine ps = let 
  2.     n = fromIntegral $ length ps
  3.     xs = map fst ps
  4.     ys = map snd ps
  5.     avg ds = (/n) $ sum ds
  6.     rval ds es 
  7.        = (avg $ zipWith (*) ds es) - (avg ds)*(avg es)
  8.     val ds = rval ds ds
  9.     avgx = avg xs
  10.     valx = val xs
  11.     avgy = avg ys
  12.     valy = val ys
  13.     rvalxy = rval xs ys
  14.     a = rvalxy/valx
  15.     b = avgy-a*avgx
  16.   in
  17.     (a,b)
  18.  
  19. dailyNewCasesInOsaka = []
  20.     ++ [74,52,55,88,40,85,54]
  21.     ++ [31,35,31,29,16,30,32]
  22.     ++ [44]
  23.  
  24. points = id
  25.     $ zip [0..] 
  26.     $ map log 
  27.     $ drop 0 dailyNewCasesInOsaka
  28.  
  29. patsAt a b x = exp $ a *x + b
  30.  
  31.  
  32. april15 = 0   + 15 - 15
  33. may06 = 30 + 06 - 15
  34.  
  35. main = do
  36.    let  (a,b) = regressionLine points
  37.    print $ "April 15th : " ++ (show $ patsAt a b april15)
  38.    print $ "The average ratio : " ++ (show $ exp a)
  39.    print $ "May 6th : " ++ (show $ patsAt a b may06)
  40.    print $ "Lists of the expecred number of new patients from May 1st"
  41.    print $ [round $ patsAt a b x | x<-[16..22]]
  42.    print $ [round $ patsAt a b x | x<-[23..29]]
  43.    print $ [round $ patsAt a b x | x<-[30..36]]
  44.    print $ [round $ patsAt a b x | x<-[37..43]]
Success #stdin #stdout 0s 4520KB
comments (?)
stdin
copy 
Standard input is empty
stdout
copy 
"April 15th : 68.97613696989895"
"The average ratio : 0.9316696264403739"
"May 6th : 15.60269248285008"
"Lists of the expecred number of new patients from May 1st"
[22,21,19,18,17,16,15]
[14,13,12,11,10,10,9]
[8,8,7,7,6,6,5]
[5,5,4,4,4,4,3]
https://ideone.com/Qhn1VZ
language:
Haskell (ghc 8.4.4)
created:
3 years ago
visibility:
public

Share or Embed source code

Discover > Sphere Engine API

The brand new service which powers Ideone!

Discover > IDE Widget

Widget for compiling and running the source code in a web browser!