library(SuppDists)
jdist<-list(gamma=-0.7, delta=0.8, xi=26, lambda=1.4, type="SU")
xvals<-seq(10,60,by=1)
yvals<-dJohnson(xvals,jdist)
print(jdist)
#plot(xvals,yvals,col="blue",type="l")
 
#this is some trial observation - say we see an observation of "40" and give it a weighting of one-third
observation<-40
observationweight<-1/3
 
#your code here - modifies jdist for example
observationnumber<-10000/observationweight
randomobservations<-rJohnson(observationnumber, jdist)
jdist<-JohnsonFit(randomobservations)
#
 
xvals<-seq(10,60,by=1)
yvals<-dJohnson(xvals,jdist)
print(jdist)
#lines(xvals,yvals,col="red")

bGlicmFyeShTdXBwRGlzdHMpCmpkaXN0PC1saXN0KGdhbW1hPS0wLjcsIGRlbHRhPTAuOCwgeGk9MjYsIGxhbWJkYT0xLjQsIHR5cGU9IlNVIikKeHZhbHM8LXNlcSgxMCw2MCxieT0xKQp5dmFsczwtZEpvaG5zb24oeHZhbHMsamRpc3QpCnByaW50KGpkaXN0KQojcGxvdCh4dmFscyx5dmFscyxjb2w9ImJsdWUiLHR5cGU9ImwiKQoKI3RoaXMgaXMgc29tZSB0cmlhbCBvYnNlcnZhdGlvbiAtIHNheSB3ZSBzZWUgYW4gb2JzZXJ2YXRpb24gb2YgIjQwIiBhbmQgZ2l2ZSBpdCBhIHdlaWdodGluZyBvZiBvbmUtdGhpcmQKb2JzZXJ2YXRpb248LTQwCm9ic2VydmF0aW9ud2VpZ2h0PC0xLzMKCiN5b3VyIGNvZGUgaGVyZSAtIG1vZGlmaWVzIGpkaXN0IGZvciBleGFtcGxlCm9ic2VydmF0aW9ubnVtYmVyPC0xMDAwMC9vYnNlcnZhdGlvbndlaWdodApyYW5kb21vYnNlcnZhdGlvbnM8LXJKb2huc29uKG9ic2VydmF0aW9ubnVtYmVyLCBqZGlzdCkKamRpc3Q8LUpvaG5zb25GaXQocmFuZG9tb2JzZXJ2YXRpb25zKQojCgp4dmFsczwtc2VxKDEwLDYwLGJ5PTEpCnl2YWxzPC1kSm9obnNvbih4dmFscyxqZGlzdCkKcHJpbnQoamRpc3QpCiNsaW5lcyh4dmFscyx5dmFscyxjb2w9InJlZCIp