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<-1000/observationweight
  14. randomobservations<-rJohnson(observationnumber, jdist)
  15. jdist<-JohnsonFit(randomobservations)
  16. #
  17.  
  18. xvals<-seq(10,60,by=1)
  19. yvals<-dJohnson(xvals,jdist)
  20. print(jdist)
  21. #lines(xvals,yvals,col="red")
Success #stdin #stdout 0.36s 24104KB
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.6796815

$delta
[1] 0.7643596

$xi
[1] 26.07033

$lambda
[1] 1.328601

$type
[1] "SU"
https://ideone.com/mIsQ4
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!