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)
postscript()
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<-5000000
randomobservations<-rJohnson(observationnumber, jdist)
newobservations<-rep(observation,times=round(observationweight*observationnumber/100))
randomobservations<-append(randomobservations,newobservations)
jdist<-JohnsonFit(randomobservations)
#
 
xvals<-seq(10,60,by=1)
yvals<-dJohnson(xvals,jdist)
print(jdist)
#lines(xvals,yvals,col="red")

bGlicmFyeShTdXBwRGlzdHMpCmpkaXN0PC1saXN0KGdhbW1hPS0wLjcsIGRlbHRhPTAuOCwgeGk9MjYsIGxhbWJkYT0xLjQsIHR5cGU9IlNVIikKeHZhbHM8LXNlcSgxMCw2MCxieT0xKQp5dmFsczwtZEpvaG5zb24oeHZhbHMsamRpc3QpCnByaW50KGpkaXN0KQpwb3N0c2NyaXB0KCkKcGxvdCh4dmFscyx5dmFscyxjb2w9ImJsdWUiLHR5cGU9ImwiKQogCiN0aGlzIGlzIHNvbWUgdHJpYWwgb2JzZXJ2YXRpb24gLSBzYXkgd2Ugc2VlIGFuIG9ic2VydmF0aW9uIG9mICI0MCIgYW5kIGdpdmUgaXQgYSB3ZWlnaHRpbmcgb2Ygb25lLXRoaXJkCm9ic2VydmF0aW9uPC00MApvYnNlcnZhdGlvbndlaWdodDwtMS8zCiAKI3lvdXIgY29kZSBoZXJlIC0gbW9kaWZpZXMgamRpc3QgZm9yIGV4YW1wbGUKb2JzZXJ2YXRpb25udW1iZXI8LTUwMDAwMDAKcmFuZG9tb2JzZXJ2YXRpb25zPC1ySm9obnNvbihvYnNlcnZhdGlvbm51bWJlciwgamRpc3QpCm5ld29ic2VydmF0aW9uczwtcmVwKG9ic2VydmF0aW9uLHRpbWVzPXJvdW5kKG9ic2VydmF0aW9ud2VpZ2h0Km9ic2VydmF0aW9ubnVtYmVyLzEwMCkpCnJhbmRvbW9ic2VydmF0aW9uczwtYXBwZW5kKHJhbmRvbW9ic2VydmF0aW9ucyxuZXdvYnNlcnZhdGlvbnMpCmpkaXN0PC1Kb2huc29uRml0KHJhbmRvbW9ic2VydmF0aW9ucykKIwogCnh2YWxzPC1zZXEoMTAsNjAsYnk9MSkKeXZhbHM8LWRKb2huc29uKHh2YWxzLGpkaXN0KQpwcmludChqZGlzdCkKI2xpbmVzKHh2YWxzLHl2YWxzLGNvbD0icmVkIik=