")
cat("Your input")
cat(" | |
---|---|

") now <- system('date +%Y-%m-%d,%T', intern=TRUE) host <- system("echo $REMOTE_ADDR", intern=TRUE) cat("Request from: ") cat(" | ")
cat(host, " at", now," ") cat(" |

") cat("Score: ") cat(" | ") cat(score) cat(" |

") cat("n: ") cat(" | ") cat(n) cat(" |

") cat("Significant digits wanted: ") cat(" | ") cat(round) cat(" |

")
cat("Results")
cat(" | |

p = ",p)
cat(" | |

")
cat("Explanation")
cat(" | |

That's the probability of scoring",score,"or more from ",n,"possible matches is",p,". If this is <0.05 then you are entitled on the conventional logic of significance testing, to say that the test has shown statistically significant evidence that something better than chance is at work here. You can set a more stringent criterion if you like, of course.")
cat("It's as simple as that: if you achieved a score with a probability lower than your preset criterion probability (alpha) then you have statistical support that something non-random is at work in your matching the objects back to their origins. Whether that will impress people will depend on whet her you have eliminated trite ways in which you might have achieved this.") cat(" If you are using this method and publishing the results or finding them of use to you, then ") linkto("I'd love to hear about what you're doing.", "http://www.psyctc.org/cgi-bin/mailto.pl?chris") ; br() cat(" | |

")
cat("Technicalities")
cat(" | |

")
cat("Calculation done in R by a CGI script written by Chris Evans using David Firth's excellent GGIwithR package.")
cat(" Neither myself, David Firth or anyone in the R team accept responsibility for the results or consequences of their use, the maths of this is fairly easy and spelled out in my publication on the topic and R is very reliable but it's always possible for things to go wrong so think through what you see here carefully before doing anything serious with it.")
linkto("Contact me if something isn't working right ...", "http://www.psyctc.org/cgi-bin/mailto.pl?chris") ; br()
lf()
cat(" ") cat("Output produced at ", date()) cat(" ") lf() cat(" |