fork download
copy 
  1. library(SuppDists)
  2. jdist<-list(gamma=-0.7, delta=0.8, xi=26, lambda=1.4, type="SU")
  3. xvals<-seq(10,60,by=1)
  4. yvals<-dJohnson(xvals,jdist)
  5. print(jdist)
  6. #plot(xvals,yvals,col="blue",type="l")
  7.  
  8. #this is some trial observation - say we see an observation of "40" and give it a weighting of one-third
  9. observation<-40
  10. observationweight<-1/3
  11.  
  12. #your code here - modifies jdist for example
  13. observationnumber<-5000000
  14. randomobservations<-rJohnson(observationnumber, jdist)
  15. newobservations<-rep(observation,times=round(observationweight*observationnumber/100))
  16. randomobservations<-append(randomobservations,newobservations)
  17. jdist<-JohnsonFit(randomobservations)
  18. #
  19.  
  20. xvals<-seq(10,60,by=1)
  21. yvals<-dJohnson(xvals,jdist)
  22. print(jdist)
  23. lines(xvals,yvals,col="red")
Runtime error #stdin #stdout 3.58s 122112KB
comments (?)
stdin
copy 
Standard input is empty
stdout
copy 
$gamma
[1] -0.7

$delta
[1] 0.8

$xi
[1] 26

$lambda
[1] 1.4

$type
[1] "SU"

$gamma
[1] -0.6893282

$delta
[1] 0.7737824

$xi
[1] 26.03562

$lambda
[1] 1.334896

$type
[1] "SU"
https://ideone.com/4ySNv
language:
R (R 3.5.2)
created:
14 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!