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)
#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
jdist<-list(gamma=-0.6, delta=0.7, xi=27, lambda=1.5, type="SU")
#
 
xvals<-seq(10,60,by=1)
yvals<-dJohnson(xvals,jdist)
#lines(xvals,yvals,col="red")

bGlicmFyeShTdXBwRGlzdHMpCmpkaXN0PC1saXN0KGdhbW1hPS0wLjcsIGRlbHRhPTAuOCwgeGk9MjYsIGxhbWJkYT0xLjQsIHR5cGU9IlNVIikKeHZhbHM8LXNlcSgxMCw2MCxieT0xKQp5dmFsczwtZEpvaG5zb24oeHZhbHMsamRpc3QpCiNwbG90KHh2YWxzLHl2YWxzLGNvbD0iYmx1ZSIsdHlwZT0ibCIpCgojdGhpcyBpcyBzb21lIHRyaWFsIG9ic2VydmF0aW9uIC0gc2F5IHdlIHNlZSBhbiBvYnNlcnZhdGlvbiBvZiAiNDAiIGFuZCBnaXZlIGl0IGEgd2VpZ2h0aW5nIG9mIG9uZS10aGlyZApvYnNlcnZhdGlvbjwtNDAKb2JzZXJ2YXRpb253ZWlnaHQ8LTEvMwoKI3lvdXIgY29kZSBoZXJlIC0gbW9kaWZpZXMgamRpc3QgZm9yIGV4YW1wbGUKamRpc3Q8LWxpc3QoZ2FtbWE9LTAuNiwgZGVsdGE9MC43LCB4aT0yNywgbGFtYmRhPTEuNSwgdHlwZT0iU1UiKQojCgp4dmFsczwtc2VxKDEwLDYwLGJ5PTEpCnl2YWxzPC1kSm9obnNvbih4dmFscyxqZGlzdCkKI2xpbmVzKHh2YWxzLHl2YWxzLGNvbD0icmVkIik=