BBS in a Box 15

home *** CD-ROM | disk | FTP | other *** search
/ BBS in a Box 15 / BBS in a box XV-2.iso / Files II / Prog / U-Z / VideoToolbox 4.95.sit / VideoToolbox / Utilities / Quick3 / TestPsychometricFit.c < prev next >
MacBinary | 1994-09-06 | 4.9 KB | [TEXT/MMCC]
open in:
MacOS 8.1 extracted |
Win98 extracted |
DOS extracted
browse contents |
view JSON data |
view as text

This file was processed as: MacBinary (archive/macBinary).
Confidence	Program	Detection	Match Type	Support
`66%`	dexvert	Compact Compressed (Unix) `(archive/compact)`	ext	Supported
`10%`	dexvert	MacBinary `(archive/macBinary)`	fallback	Supported
`1%`	dexvert	Text File `(text/txt)`	fallback	Supported
`100%`	file	MacBinary II, Tue Sep 6 00:33:56 1994, modified Tue Sep 6 00:33:56 1994, creator 'MMCC', type ASCII, 4409 bytes "TestPsychometricFit.c" , at 0x11b9 382 bytes resource	default (weak)
`99%`	file	data	default
`74%`	TrID	Macintosh plain text (MacBinary)	default
`25%`	TrID	MacBinary 2	default (weak)
`100%`	siegfried	fmt/1762 MacBinary (II)	default
`100%`	lsar	MacBinary	default
id metadata
key	value
macFileType	`[TEXT]`
macFileCreator	`[MMCC]`
hex view
+--------+-------------------------+-------------------------+--------+--------+
|00000000| 00 15 54 65 73 74 50 73 | 79 63 68 6f 6d 65 74 72 |..TestPs|ychometr|
|00000010| 69 63 46 69 74 2e 63 00 | 00 00 00 00 00 00 00 00 |icFit.c.|........|
|00000020| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000030| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000040| 00 54 45 58 54 4d 4d 43 | 43 00 00 00 00 00 00 00 |.TEXTMMC|C.......|
|00000050| 00 00 00 00 00 11 39 00 | 00 01 7e aa 91 a1 b4 aa |......9.|..~.....|
|00000060| 91 a1 b4 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00000070| 00 00 00 00 00 00 00 00 | 00 00 81 81 9d 28 00 00 |........|.....(..|
|00000080| 2f 2a 0d 54 65 73 74 50 | 73 79 63 68 6f 6d 65 74 |/*.TestP|sychomet|
|00000090| 72 69 63 46 69 74 2e 63 | 0d 43 6f 70 79 72 69 67 |ricFit.c|.Copyrig|
|000000a0| 68 74 20 28 63 29 20 31 | 39 39 30 2d 31 39 39 32 |ht (c) 1|990-1992|
|000000b0| 20 44 65 6e 69 73 20 47 | 2e 20 50 65 6c 6c 69 0d | Denis G|. Pelli.|
|000000c0| 54 68 69 73 20 69 73 20 | 61 20 73 69 6d 70 6c 65 |This is |a simple|
|000000d0| 20 64 72 69 76 65 72 20 | 74 6f 20 73 68 6f 77 20 | driver |to show |
|000000e0| 74 68 61 74 20 50 73 79 | 63 68 6f 6d 65 74 72 69 |that Psy|chometri|
|000000f0| 63 46 69 74 2e 63 20 77 | 6f 72 6b 73 2e 0d 0d 49 |cFit.c w|orks...I|
|00000100| 20 75 73 65 20 74 68 65 | 20 57 65 69 62 75 6c 6c | use the| Weibull|
|00000110| 20 70 73 79 63 68 6f 6d | 65 74 72 69 63 20 66 75 | psychom|etric fu|
|00000120| 6e 63 74 69 6f 6e 20 74 | 6f 20 70 72 6f 76 69 64 |nction t|o provid|
|00000130| 65 20 74 68 65 20 70 72 | 6f 62 61 62 69 6c 69 74 |e the pr|obabilit|
|00000140| 79 20 61 74 20 65 61 63 | 68 20 63 6f 6e 74 72 61 |y at eac|h contra|
|00000150| 73 74 0d 61 6e 64 20 74 | 68 65 20 4e 75 6d 65 72 |st.and t|he Numer|
|00000160| 69 63 61 6c 20 52 65 63 | 69 70 65 73 20 62 6e 6c |ical Rec|ipes bnl|
|00000170| 64 65 76 28 29 20 62 69 | 6e 6f 6d 69 61 6c 20 64 |dev() bi|nomial d|
|00000180| 65 76 69 61 74 65 73 20 | 66 75 6e 63 74 69 6f 6e |eviates |function|
|00000190| 20 74 6f 20 73 69 6d 75 | 6c 61 74 65 20 74 68 65 | to simu|late the|
|000001a0| 20 61 70 70 72 6f 70 72 | 69 61 74 65 0d 6e 75 6d | appropr|iate.num|
|000001b0| 62 65 72 20 6f 66 20 74 | 72 69 61 6c 73 20 61 74 |ber of t|rials at|
|000001c0| 20 74 68 61 74 20 63 6f | 6e 74 72 61 73 74 2e 0d | that co|ntrast..|
|000001d0| 0d 54 68 65 6e 20 49 20 | 63 61 6c 6c 20 50 73 79 |.Then I |call Psy|
|000001e0| 63 68 6f 6d 65 74 72 69 | 63 46 69 74 28 29 20 61 |chometri|cFit() a|
|000001f0| 6e 64 20 61 73 6b 20 69 | 74 20 74 6f 20 66 69 74 |nd ask i|t to fit|
|00000200| 20 74 68 65 20 57 65 69 | 62 75 6c 6c 28 29 20 66 | the Wei|bull() f|
|00000210| 75 6e 63 74 69 6f 6e 20 | 74 6f 20 74 68 65 20 64 |unction |to the d|
|00000220| 61 74 61 2c 20 61 6e 64 | 0d 74 68 65 6e 20 49 20 |ata, and|.then I |
|00000230| 70 72 69 6e 74 20 6f 75 | 74 20 74 68 65 20 72 65 |print ou|t the re|
|00000240| 73 75 6c 74 73 2e 20 41 | 73 20 6f 6e 65 20 77 6f |sults. A|s one wo|
|00000250| 75 6c 64 20 65 78 70 65 | 63 74 2c 20 74 68 65 20 |uld expe|ct, the |
|00000260| 66 69 74 73 20 61 72 65 | 20 76 65 72 79 20 67 6f |fits are| very go|
|00000270| 6f 64 2e 20 49 74 20 73 | 65 65 6d 73 0d 72 65 61 |od. It s|eems.rea|
|00000280| 73 6f 6e 61 62 6c 65 20 | 74 6f 20 72 65 6a 65 63 |sonable |to rejec|
|00000290| 74 20 66 69 74 73 20 61 | 74 20 74 68 65 20 35 25 |t fits a|t the 5%|
|000002a0| 20 73 69 67 6e 69 66 69 | 63 61 6e 63 65 20 6c 65 | signifi|cance le|
|000002b0| 76 65 6c 2e 20 53 6f 20 | 65 78 70 65 63 74 20 74 |vel. So |expect t|
|000002c0| 6f 20 72 65 6a 65 63 74 | 20 61 62 6f 75 74 0d 35 |o reject| about.5|
|000002d0| 25 20 6f 66 20 79 6f 75 | 72 20 66 69 74 73 20 65 |% of you|r fits e|
|000002e0| 76 65 6e 20 69 66 20 61 | 6c 6c 20 74 68 65 20 61 |ven if a|ll the a|
|000002f0| 73 73 75 6d 70 74 69 6f | 6e 73 20 6f 66 20 74 68 |ssumptio|ns of th|
|00000300| 69 73 20 6d 6f 64 65 6c | 20 61 72 65 20 63 6f 72 |is model| are cor|
|00000310| 72 65 63 74 2e 0d 0d 4f | 70 74 69 6f 6e 61 6c 6c |rect...O|ptionall|
|00000320| 79 2c 20 61 6c 73 6f 20 | 63 6f 6d 70 61 72 65 73 |y, also |compares|
|00000330| 20 51 55 49 43 4b 33 20 | 77 69 74 68 20 51 55 45 | QUICK3 |with QUE|
|00000340| 53 54 2e 20 28 4e 6f 74 | 65 3a 20 51 55 45 53 54 |ST. (Not|e: QUEST|
|00000350| 20 69 73 20 6e 6f 74 20 | 70 61 72 74 20 6f 66 20 | is not |part of |
|00000360| 74 68 65 20 0d 56 69 64 | 65 6f 54 6f 6f 6c 62 6f |the .Vid|eoToolbo|
|00000370| 78 2e 20 53 6f 72 72 79 | 2e 29 0d 0d 48 49 53 54 |x. Sorry|.)..HIST|
|00000380| 4f 52 59 3a 0d 34 2f 36 | 2f 39 30 09 64 67 70 09 |ORY:.4/6|/90.dgp.|
|00000390| 77 72 6f 74 65 20 69 74 | 2e 20 53 65 65 6d 73 20 |wrote it|. Seems |
|000003a0| 74 6f 20 77 6f 72 6b 20 | 66 69 6e 65 20 66 6f 72 |to work |fine for|
|000003b0| 20 61 6c 6c 20 63 61 73 | 65 73 2c 20 31 20 74 6f | all cas|es, 1 to|
|000003c0| 20 31 30 30 2c 30 30 30 | 2c 30 30 30 20 74 72 69 | 100,000|,000 tri|
|000003d0| 61 6c 73 20 70 65 72 0d | 09 09 09 63 6f 6e 74 72 |als per.|...contr|
|000003e0| 61 73 74 2c 20 61 6e 64 | 20 66 65 77 20 61 6e 64 |ast, and| few and|
|000003f0| 20 6d 61 6e 79 20 63 6f | 6e 74 72 61 73 74 73 2e | many co|ntrasts.|
|00000400| 0d 31 31 2f 31 38 2f 39 | 32 20 64 67 70 20 61 64 |.11/18/9|2 dgp ad|
|00000410| 64 65 64 20 63 6f 6d 70 | 61 72 69 73 6f 6e 20 77 |ded comp|arison w|
|00000420| 69 74 68 20 51 55 45 53 | 54 2e 0d 32 2f 32 30 2f |ith QUES|T..2/20/|
|00000430| 39 33 09 64 67 70 09 61 | 64 64 65 64 20 63 61 6c |93.dgp.a|dded cal|
|00000440| 6c 20 74 6f 20 52 65 71 | 75 69 72 65 28 29 2e 0d |l to Req|uire()..|
|00000450| 39 2f 35 2f 39 34 20 64 | 67 70 20 72 65 6d 6f 76 |9/5/94 d|gp remov|
|00000460| 65 64 20 61 73 73 75 6d | 70 74 69 6f 6e 20 69 6e |ed assum|ption in|
|00000470| 20 70 72 69 6e 74 66 27 | 73 20 74 68 61 74 20 69 | printf'|s that i|
|00000480| 6e 74 3d 3d 73 68 6f 72 | 74 2e 0d 2a 2f 0d 23 69 |nt==shor|t..*/.#i|
|00000490| 6e 63 6c 75 64 65 20 22 | 56 69 64 65 6f 54 6f 6f |nclude "|VideoToo|
|000004a0| 6c 62 6f 78 2e 68 22 0d | 23 64 65 66 69 6e 65 20 |lbox.h".|#define |
|000004b0| 51 55 45 53 54 20 30 0d | 23 69 66 20 51 55 45 53 |QUEST 0.|#if QUES|
|000004c0| 54 0d 09 23 69 6e 63 6c | 75 64 65 20 22 51 75 65 |T..#incl|ude "Que|
|000004d0| 73 74 2e 68 22 0d 23 65 | 6e 64 69 66 0d 23 69 6e |st.h".#e|ndif.#in|
|000004e0| 63 6c 75 64 65 20 22 51 | 75 69 63 6b 33 2e 68 22 |clude "Q|uick3.h"|
|000004f0| 0d 23 69 6e 63 6c 75 64 | 65 20 3c 61 73 73 65 72 |.#includ|e <asser|
|00000500| 74 2e 68 3e 0d 23 69 6e | 63 6c 75 64 65 20 3c 74 |t.h>.#in|clude <t|
|00000510| 69 6d 65 2e 68 3e 0d 23 | 69 6e 63 6c 75 64 65 20 |ime.h>.#|include |
|00000520| 22 6e 72 2e 68 22 09 09 | 2f 2a 20 4e 75 6d 65 72 |"nr.h"..|/* Numer|
|00000530| 69 63 61 6c 20 52 65 63 | 69 70 65 73 20 69 6e 20 |ical Rec|ipes in |
|00000540| 43 20 2a 2f 0d 0d 76 6f | 69 64 20 6d 61 69 6e 28 |C */..vo|id main(|
|00000550| 76 6f 69 64 29 0d 7b 0d | 09 64 61 74 61 52 65 63 |void).{.|.dataRec|
|00000560| 6f 72 64 20 64 61 74 61 | 2c 6d 6f 6e 6f 74 6f 6e |ord data|,monoton|
|00000570| 69 63 44 61 74 61 3b 0d | 09 63 6f 6e 74 72 61 73 |icData;.|.contras|
|00000580| 74 52 65 63 6f 72 64 20 | 2a 63 50 74 72 3b 0d 09 |tRecord |*cPtr;..|
|00000590| 70 61 72 61 6d 52 65 63 | 6f 72 64 20 70 61 72 61 |paramRec|ord para|
|000005a0| 6d 73 3b 0d 09 69 6e 74 | 20 69 3b 0d 09 64 6f 75 |ms;..int| i;..dou|
|000005b0| 62 6c 65 20 63 68 69 53 | 71 75 61 72 65 2c 77 65 |ble chiS|quare,we|
|000005c0| 69 62 75 6c 6c 4c 4c 2c | 6d 6f 6e 6f 74 6f 6e 69 |ibullLL,|monotoni|
|000005d0| 63 4c 4c 3b 09 2f 2a 20 | 6c 6f 67 20 6c 69 6b 65 |cLL;./* |log like|
|000005e0| 6c 69 68 6f 6f 64 20 2a | 2f 0d 09 69 6e 74 20 63 |lihood *|/..int c|
|000005f0| 68 69 53 71 75 61 72 65 | 44 46 2c 77 65 69 62 75 |hiSquare|DF,weibu|
|00000600| 6c 6c 44 46 2c 6d 6f 6e | 6f 74 6f 6e 69 63 44 46 |llDF,mon|otonicDF|
|00000610| 3b 09 2f 2a 20 64 65 67 | 72 65 65 73 20 6f 66 20 |;./* deg|rees of |
|00000620| 66 72 65 65 64 6f 6d 20 | 2a 2f 0d 09 64 6f 75 62 |freedom |*/..doub|
|00000630| 6c 65 20 70 2c 6c 6f 67 | 43 2c 72 61 6e 67 65 3b |le p,log|C,range;|
|00000640| 0d 09 6c 6f 6e 67 20 69 | 64 75 6d 3b 0d 09 69 6e |..long i|dum;..in|
|00000650| 74 20 63 6f 6e 64 3d 30 | 3b 0d 09 23 69 66 20 51 |t cond=0|;..#if Q|
|00000660| 55 45 53 54 0d 09 09 51 | 75 65 73 74 20 2a 71 3b |UEST...Q|uest *q;|
|00000670| 0d 09 23 65 6e 64 69 66 | 0d 09 0d 09 52 65 71 75 |..#endif|....Requ|
|00000680| 69 72 65 28 30 29 3b 0d | 09 70 61 72 61 6d 73 2e |ire(0);.|.params.|
|00000690| 6c 6f 67 41 6c 70 68 61 | 3d 30 3b 0d 09 70 61 72 |logAlpha|=0;..par|
|000006a0| 61 6d 73 2e 62 65 74 61 | 3d 33 3b 0d 09 70 61 72 |ams.beta|=3;..par|
|000006b0| 61 6d 73 2e 67 61 6d 6d | 61 3d 30 2e 31 31 31 31 |ams.gamm|a=0.1111|
|000006c0| 37 33 3b 20 0d 09 70 61 | 72 61 6d 73 2e 64 65 6c |73; ..pa|rams.del|
|000006d0| 74 61 3d 30 2e 30 31 3b | 0d 09 0d 09 69 64 75 6d |ta=0.01;|....idum|
|000006e0| 3d 74 69 6d 65 28 4e 55 | 4c 4c 29 3b 0d 09 70 72 |=time(NU|LL);..pr|
|000006f0| 69 6e 74 66 28 22 52 61 | 6e 64 6f 6d 20 73 65 65 |intf("Ra|ndom see|
|00000700| 64 20 25 6c 64 5c 6e 22 | 2c 69 64 75 6d 29 3b 09 |d %ld\n"|,idum);.|
|00000710| 2f 2a 20 53 6f 20 77 65 | 20 63 61 6e 20 72 65 70 |/* So we| can rep|
|00000720| 72 6f 64 75 63 65 20 69 | 6e 74 65 72 65 73 74 69 |roduce i|nteresti|
|00000730| 6e 67 20 63 61 73 65 73 | 20 2a 2f 0d 09 64 61 74 |ng cases| */..dat|
|00000740| 61 2e 63 6f 6e 74 72 61 | 73 74 73 3d 35 3b 0d 09 |a.contra|sts=5;..|
|00000750| 72 61 6e 67 65 3d 30 2e | 35 3b 09 09 09 09 09 09 |range=0.|5;......|
|00000760| 09 2f 2a 20 6c 6f 67 20 | 63 6f 6e 74 72 61 73 74 |./* log |contrast|
|00000770| 20 72 61 6e 67 65 2c 20 | 63 65 6e 74 65 72 65 64 | range, |centered|
|00000780| 20 6f 6e 20 6c 6f 67 41 | 6c 70 68 61 20 2a 2f 0d | on logA|lpha */.|
|00000790| 09 72 61 6e 67 65 3d 32 | 3b 0d 09 66 6f 72 28 69 |.range=2|;..for(i|
|000007a0| 3d 30 3b 69 3c 64 61 74 | 61 2e 63 6f 6e 74 72 61 |=0;i<dat|a.contra|
|000007b0| 73 74 73 3b 69 2b 2b 29 | 7b 0d 09 09 6c 6f 67 43 |sts;i++)|{...logC|
|000007c0| 3d 70 61 72 61 6d 73 2e | 6c 6f 67 41 6c 70 68 61 |=params.|logAlpha|
|000007d0| 20 2b 20 72 61 6e 67 65 | 2a 28 69 2f 28 64 61 74 | + range|*(i/(dat|
|000007e0| 61 2e 63 6f 6e 74 72 61 | 73 74 73 2d 31 2e 30 29 |a.contra|sts-1.0)|
|000007f0| 20 2d 20 30 2e 35 29 3b | 0d 09 09 64 61 74 61 2e | - 0.5);|...data.|
|00000800| 63 5b 69 5d 2e 63 6f 6e | 74 72 61 73 74 3d 70 6f |c[i].con|trast=po|
|00000810| 77 28 31 30 2e 30 2c 6c | 6f 67 43 29 3b 0d 09 09 |w(10.0,l|ogC);...|
|00000820| 64 61 74 61 2e 63 5b 69 | 5d 2e 74 72 69 61 6c 73 |data.c[i|].trials|
|00000830| 3d 31 30 3b 09 09 09 09 | 2f 2a 20 74 72 69 61 6c |=10;....|/* trial|
|00000840| 73 20 61 74 20 74 68 69 | 73 20 63 6f 6e 74 72 61 |s at thi|s contra|
|00000850| 73 74 20 2a 2f 0d 09 09 | 64 61 74 61 2e 63 5b 69 |st */...|data.c[i|
|00000860| 5d 2e 63 6f 72 72 65 63 | 74 3d 62 6e 6c 64 65 76 |].correc|t=bnldev|
|00000870| 28 57 65 69 62 75 6c 6c | 28 64 61 74 61 2e 63 5b |(Weibull|(data.c[|
|00000880| 69 5d 2e 63 6f 6e 74 72 | 61 73 74 2c 26 70 61 72 |i].contr|ast,&par|
|00000890| 61 6d 73 29 2c 64 61 74 | 61 2e 63 5b 69 5d 2e 74 |ams),dat|a.c[i].t|
|000008a0| 72 69 61 6c 73 2c 26 69 | 64 75 6d 29 3b 0d 09 7d |rials,&i|dum);..}|
|000008b0| 0d 09 23 69 66 20 51 55 | 45 53 54 0d 09 09 2f 2a |..#if QU|EST.../*|
|000008c0| 20 51 75 65 73 74 20 50 | 61 72 61 6d 65 74 65 72 | Quest P|arameter|
|000008d0| 73 20 2a 2f 0d 09 09 71 | 3d 28 51 75 65 73 74 20 |s */...q|=(Quest |
|000008e0| 2a 29 6d 61 6c 6c 6f 63 | 28 73 69 7a 65 6f 66 28 |*)malloc|(sizeof(|
|000008f0| 51 75 65 73 74 29 29 3b | 0d 09 09 61 73 73 65 72 |Quest));|...asser|
|00000900| 74 28 71 21 3d 4e 55 4c | 4c 29 3b 0d 09 09 71 2d |t(q!=NUL|L);...q-|
|00000910| 3e 6e 43 6f 6e 64 73 3d | 31 3b 09 0d 09 09 71 2d |>nConds=|1;....q-|
|00000920| 3e 6e 4c 65 76 65 6c 73 | 3d 36 30 30 3b 0d 09 09 |>nLevels|=600;...|
|00000930| 71 2d 3e 6e 54 72 69 61 | 6c 73 3d 30 3b 0d 09 09 |q->nTria|ls=0;...|
|00000940| 71 2d 3e 67 72 61 69 6e | 3d 30 2e 30 31 3b 09 2f |q->grain|=0.01;./|
|00000950| 2a 20 73 74 65 70 20 73 | 69 7a 65 20 6f 66 20 67 |* step s|ize of g|
|00000960| 72 69 64 2c 20 69 6e 20 | 6c 6f 67 20 63 6f 6e 74 |rid, in |log cont|
|00000970| 72 61 73 74 20 2a 2f 0d | 09 09 71 2d 3e 69 6e 69 |rast */.|..q->ini|
|00000980| 74 69 61 6c 53 44 3d 31 | 3b 0d 09 20 09 71 2d 3e |tialSD=1|;.. .q->|
|00000990| 6e 52 65 73 70 6f 6e 73 | 65 73 3d 32 3b 0d 09 20 |nRespons|es=2;.. |
|000009a0| 09 71 2d 3e 71 75 61 6e | 74 69 6c 65 4f 72 64 65 |.q->quan|tileOrde|
|000009b0| 72 3d 4e 41 4e 3b 0d 09 | 20 09 71 2d 3e 66 61 6b |r=NAN;..| .q->fak|
|000009c0| 65 49 74 3d 30 3b 0d 09 | 20 09 71 2d 3e 66 75 6e |eIt=0;..| .q->fun|
|000009d0| 63 74 69 6f 6e 3d 57 65 | 69 62 75 6c 6c 50 52 65 |ction=We|ibullPRe|
|000009e0| 73 70 6f 6e 73 65 3b 0d | 09 09 71 2d 3e 62 65 74 |sponse;.|..q->bet|
|000009f0| 61 3d 70 61 72 61 6d 73 | 2e 62 65 74 61 3b 0d 09 |a=params|.beta;..|
|00000a00| 09 71 2d 3e 67 61 6d 6d | 61 3d 70 61 72 61 6d 73 |.q->gamm|a=params|
|00000a10| 2e 67 61 6d 6d 61 3b 0d | 09 09 71 2d 3e 64 65 6c |.gamma;.|..q->del|
|00000a20| 74 61 3d 70 61 72 61 6d | 73 2e 64 65 6c 74 61 3b |ta=param|s.delta;|
|00000a30| 0d 09 09 71 2d 3e 65 70 | 73 69 6c 6f 6e 3d 30 2e |...q->ep|silon=0.|
|00000a40| 30 3b 0d 09 09 63 6f 6e | 64 3d 30 3b 0d 09 09 71 |0;...con|d=0;...q|
|00000a50| 2d 3e 67 75 65 73 73 5b | 63 6f 6e 64 5d 3d 30 3b |->guess[|cond]=0;|
|00000a60| 0d 09 09 51 75 65 73 74 | 4f 70 65 6e 28 71 29 3b |...Quest|Open(q);|
|00000a70| 0d 09 09 66 6f 72 28 69 | 3d 30 3b 69 3c 64 61 74 |...for(i|=0;i<dat|
|00000a80| 61 2e 63 6f 6e 74 72 61 | 73 74 73 3b 69 2b 2b 29 |a.contra|sts;i++)|
|00000a90| 7b 0d 09 09 09 6c 6f 67 | 43 3d 6c 6f 67 31 30 28 |{....log|C=log10(|
|00000aa0| 64 61 74 61 2e 63 5b 69 | 5d 2e 63 6f 6e 74 72 61 |data.c[i|].contra|
|00000ab0| 73 74 29 3b 0d 09 09 09 | 66 6f 72 28 6a 3d 30 3b |st);....|for(j=0;|
|00000ac0| 6a 3c 64 61 74 61 2e 63 | 5b 69 5d 2e 63 6f 72 72 |j<data.c|[i].corr|
|00000ad0| 65 63 74 3b 6a 2b 2b 29 | 51 75 65 73 74 55 70 64 |ect;j++)|QuestUpd|
|00000ae0| 61 74 65 28 71 2c 63 6f | 6e 64 2c 6c 6f 67 43 2c |ate(q,co|nd,logC,|
|00000af0| 31 29 3b 0d 09 09 09 66 | 6f 72 28 3b 6a 3c 64 61 |1);....f|or(;j<da|
|00000b00| 74 61 2e 63 5b 69 5d 2e | 74 72 69 61 6c 73 3b 6a |ta.c[i].|trials;j|
|00000b10| 2b 2b 29 51 75 65 73 74 | 55 70 64 61 74 65 28 71 |++)Quest|Update(q|
|00000b20| 2c 63 6f 6e 64 2c 6c 6f | 67 43 2c 30 29 3b 0d 09 |,cond,lo|gC,0);..|
|00000b30| 09 7d 0d 09 09 71 5f 72 | 65 6d 6f 76 65 50 72 69 |.}...q_r|emovePri|
|00000b40| 6f 72 28 71 2d 3e 71 43 | 6f 6e 64 73 5b 63 6f 6e |or(q->qC|onds[con|
|00000b50| 64 5d 29 3b 0d 09 09 6d | 6f 64 65 20 3d 20 71 5f |d]);...m|ode = q_|
|00000b60| 6d 6f 64 65 28 71 2d 3e | 71 43 6f 6e 64 73 5b 63 |mode(q->|qConds[c|
|00000b70| 6f 6e 64 5d 29 3b 0d 09 | 09 51 75 65 73 74 43 6c |ond]);..|.QuestCl|
|00000b80| 6f 73 65 28 71 29 3b 0d | 09 23 65 6e 64 69 66 0d |ose(q);.|.#endif.|
|00000b90| 09 70 72 69 6e 74 66 28 | 22 54 65 73 74 69 6e 67 |.printf(|"Testing|
|00000ba0| 20 74 68 65 20 66 75 6e | 63 74 69 6f 6e 20 50 73 | the fun|ction Ps|
|00000bb0| 79 63 68 6f 6d 65 74 72 | 69 63 46 69 74 2e 5c 6e |ychometr|icFit.\n|
|00000bc0| 22 29 3b 0d 09 70 72 69 | 6e 74 66 28 22 53 69 6d |");..pri|ntf("Sim|
|00000bd0| 75 6c 61 74 69 6e 67 20 | 61 6e 20 6f 62 73 65 72 |ulating |an obser|
|00000be0| 76 65 72 20 77 69 74 68 | 20 61 20 57 65 69 62 75 |ver with| a Weibu|
|00000bf0| 6c 6c 20 70 73 79 63 68 | 6f 6d 65 74 72 69 63 20 |ll psych|ometric |
|00000c00| 66 75 6e 63 74 69 6f 6e | 2e 5c 6e 22 29 3b 0d 09 |function|.\n");..|
|00000c10| 77 65 69 62 75 6c 6c 44 | 46 3d 31 3b 09 2f 2a 20 |weibullD|F=1;./* |
|00000c20| 6e 75 6d 62 65 72 20 6f | 66 20 70 61 72 61 6d 65 |number o|f parame|
|00000c30| 74 65 72 73 20 74 6f 20 | 62 65 20 61 64 6a 75 73 |ters to |be adjus|
|00000c40| 74 65 64 20 69 6e 20 66 | 69 74 74 69 6e 67 20 2a |ted in f|itting *|
|00000c50| 2f 0d 09 70 72 69 6e 74 | 66 28 22 54 68 65 20 73 |/..print|f("The s|
|00000c60| 69 6d 75 6c 61 74 65 64 | 20 64 61 74 61 20 77 69 |imulated| data wi|
|00000c70| 6c 6c 20 62 65 20 66 69 | 74 20 75 73 69 6e 67 20 |ll be fi|t using |
|00000c80| 25 64 20 64 65 67 72 65 | 65 73 20 6f 66 20 66 72 |%d degre|es of fr|
|00000c90| 65 65 64 6f 6d 2e 5c 6e | 22 2c 28 69 6e 74 29 77 |eedom.\n|",(int)w|
|00000ca0| 65 69 62 75 6c 6c 44 46 | 29 3b 0d 09 70 72 69 6e |eibullDF|);..prin|
|00000cb0| 74 66 28 22 4f 62 73 65 | 72 76 65 72 3a 20 22 29 |tf("Obse|rver: ")|
|00000cc0| 3b 0d 09 70 72 69 6e 74 | 66 28 22 6c 6f 67 41 6c |;..print|f("logAl|
|00000cd0| 70 68 61 25 36 2e 32 66 | 2c 20 62 65 74 61 25 34 |pha%6.2f|, beta%4|
|00000ce0| 2e 31 66 2c 20 67 61 6d | 6d 61 25 35 2e 32 66 2c |.1f, gam|ma%5.2f,|
|00000cf0| 20 64 65 6c 74 61 25 35 | 2e 32 66 5c 6e 22 2c 70 | delta%5|.2f\n",p|
|00000d00| 61 72 61 6d 73 2e 6c 6f | 67 41 6c 70 68 61 2c 70 |arams.lo|gAlpha,p|
|00000d10| 61 72 61 6d 73 2e 62 65 | 74 61 2c 70 61 72 61 6d |arams.be|ta,param|
|00000d20| 73 2e 67 61 6d 6d 61 2c | 70 61 72 61 6d 73 2e 64 |s.gamma,|params.d|
|00000d30| 65 6c 74 61 29 3b 0d 09 | 70 3d 50 73 79 63 68 6f |elta);..|p=Psycho|
|00000d40| 6d 65 74 72 69 63 46 69 | 74 28 26 70 61 72 61 6d |metricFi|t(&param|
|00000d50| 73 2c 26 57 65 69 62 75 | 6c 6c 2c 26 64 61 74 61 |s,&Weibu|ll,&data|
|00000d60| 2c 26 77 65 69 62 75 6c | 6c 4c 4c 2c 77 65 69 62 |,&weibul|lLL,weib|
|00000d70| 75 6c 6c 44 46 2c 26 63 | 68 69 53 71 75 61 72 65 |ullDF,&c|hiSquare|
|00000d80| 2c 26 63 68 69 53 71 75 | 61 72 65 44 46 29 3b 0d |,&chiSqu|areDF);.|
|00000d90| 09 70 72 69 6e 74 66 28 | 22 46 69 74 3a 20 20 20 |.printf(|"Fit:   |
|00000da0| 20 20 20 22 29 3b 0d 09 | 70 72 69 6e 74 66 28 22 |   ");..|printf("|
|00000db0| 6c 6f 67 41 6c 70 68 61 | 25 36 2e 32 66 2c 20 62 |logAlpha|%6.2f, b|
|00000dc0| 65 74 61 25 34 2e 31 66 | 2c 20 67 61 6d 6d 61 25 |eta%4.1f|, gamma%|
|00000dd0| 35 2e 32 66 2c 20 64 65 | 6c 74 61 25 35 2e 32 66 |5.2f, de|lta%5.2f|
|00000de0| 2c 20 73 69 67 6e 69 66 | 69 63 61 6e 63 65 25 35 |, signif|icance%5|
|00000df0| 2e 32 66 5c 6e 22 2c 70 | 61 72 61 6d 73 2e 6c 6f |.2f\n",p|arams.lo|
|00000e00| 67 41 6c 70 68 61 2c 70 | 61 72 61 6d 73 2e 62 65 |gAlpha,p|arams.be|
|00000e10| 74 61 2c 70 61 72 61 6d | 73 2e 67 61 6d 6d 61 2c |ta,param|s.gamma,|
|00000e20| 70 61 72 61 6d 73 2e 64 | 65 6c 74 61 2c 70 29 3b |params.d|elta,p);|
|00000e30| 0d 09 23 69 66 20 51 55 | 45 53 54 0d 09 09 70 72 |..#if QU|EST...pr|
|00000e40| 69 6e 74 66 28 22 51 55 | 45 53 54 20 6d 6f 64 65 |intf("QU|EST mode|
|00000e50| 20 25 2e 32 66 5c 6e 22 | 2c 6d 6f 64 65 29 3b 0d | %.2f\n"|,mode);.|
|00000e60| 09 23 65 6e 64 69 66 0d | 0d 09 2f 2a 0d 09 57 65 |.#endif.|../*..We|
|00000e70| 27 72 65 20 64 6f 6e 65 | 2c 20 62 75 74 20 6a 75 |'re done|, but ju|
|00000e80| 73 74 20 74 6f 20 73 68 | 6f 77 20 6f 66 66 2c 20 |st to sh|ow off, |
|00000e90| 6c 65 74 27 73 20 70 72 | 69 6e 74 20 6f 75 74 20 |let's pr|int out |
|00000ea0| 65 76 65 72 79 74 68 69 | 6e 67 20 74 68 61 74 20 |everythi|ng that |
|00000eb0| 61 6e 79 6f 6e 65 0d 09 | 63 6f 75 6c 64 20 70 6f |anyone..|could po|
|00000ec0| 73 73 69 62 6c 79 20 77 | 61 6e 74 2e 20 49 6e 20 |ssibly w|ant. In |
|00000ed0| 72 65 61 6c 20 6c 69 66 | 65 20 49 20 77 6f 75 6c |real lif|e I woul|
|00000ee0| 64 20 73 6b 69 70 20 74 | 68 69 73 20 6a 75 6e 6b |d skip t|his junk|
|00000ef0| 2e 0d 09 2a 2f 0d 09 6d | 6f 6e 6f 74 6f 6e 69 63 |...*/..m|onotonic|
|00000f00| 44 61 74 61 3d 64 61 74 | 61 3b 0d 09 4d 6f 6e 6f |Data=dat|a;..Mono|
|00000f10| 74 6f 6e 69 63 46 69 74 | 28 26 6d 6f 6e 6f 74 6f |tonicFit|(&monoto|
|00000f20| 6e 69 63 44 61 74 61 2c | 26 6d 6f 6e 6f 74 6f 6e |nicData,|&monoton|
|00000f30| 69 63 4c 4c 2c 26 6d 6f | 6e 6f 74 6f 6e 69 63 44 |icLL,&mo|notonicD|
|00000f40| 46 29 3b 09 2f 2a 20 6f | 76 65 72 77 72 69 74 65 |F);./* o|verwrite|
|00000f50| 73 20 64 61 74 61 20 77 | 69 74 68 20 66 69 74 20 |s data w|ith fit |
|00000f60| 2a 2f 0d 09 70 72 69 6e | 74 66 28 22 5c 6e 63 6f |*/..prin|tf("\nco|
|00000f70| 6e 74 72 61 73 74 20 54 | 72 69 61 6c 73 20 52 69 |ntrast T|rials Ri|
|00000f80| 67 68 74 20 20 52 61 74 | 69 6f 20 20 57 65 69 62 |ght  Rat|io  Weib|
|00000f90| 75 6c 6c 20 4d 6f 6e 6f | 74 6f 6e 65 5c 6e 22 29 |ull Mono|tone\n")|
|00000fa0| 3b 0d 09 66 6f 72 28 69 | 3d 30 3b 69 3c 64 61 74 |;..for(i|=0;i<dat|
|00000fb0| 61 2e 63 6f 6e 74 72 61 | 73 74 73 3b 69 2b 2b 29 |a.contra|sts;i++)|
|00000fc0| 7b 0d 09 09 63 50 74 72 | 3d 26 64 61 74 61 2e 63 |{...cPtr|=&data.c|
|00000fd0| 5b 69 5d 3b 0d 09 09 70 | 72 69 6e 74 66 28 22 25 |[i];...p|rintf("%|
|00000fe0| 36 2e 33 66 20 20 20 25 | 35 6c 64 20 25 35 6c 64 |6.3f   %|5ld %5ld|
|00000ff0| 20 25 37 2e 33 66 20 25 | 37 2e 33 66 20 25 37 2e | %7.3f %|7.3f %7.|
|00001000| 33 66 5c 6e 22 2c 0d 09 | 09 09 63 50 74 72 2d 3e |3f\n",..|..cPtr->|
|00001010| 63 6f 6e 74 72 61 73 74 | 2c 63 50 74 72 2d 3e 74 |contrast|,cPtr->t|
|00001020| 72 69 61 6c 73 2c 63 50 | 74 72 2d 3e 63 6f 72 72 |rials,cP|tr->corr|
|00001030| 65 63 74 2c 0d 09 09 09 | 63 50 74 72 2d 3e 63 6f |ect,....|cPtr->co|
|00001040| 72 72 65 63 74 2f 28 64 | 6f 75 62 6c 65 29 63 50 |rrect/(d|ouble)cP|
|00001050| 74 72 2d 3e 74 72 69 61 | 6c 73 2c 0d 09 09 09 57 |tr->tria|ls,....W|
|00001060| 65 69 62 75 6c 6c 28 63 | 50 74 72 2d 3e 63 6f 6e |eibull(c|Ptr->con|
|00001070| 74 72 61 73 74 2c 26 70 | 61 72 61 6d 73 29 2c 0d |trast,&p|arams),.|
|00001080| 09 09 09 6d 6f 6e 6f 74 | 6f 6e 69 63 44 61 74 61 |...monot|onicData|
|00001090| 2e 63 5b 69 5d 2e 63 6f | 72 72 65 63 74 2f 28 64 |.c[i].co|rrect/(d|
|000010a0| 6f 75 62 6c 65 29 6d 6f | 6e 6f 74 6f 6e 69 63 44 |ouble)mo|notonicD|
|000010b0| 61 74 61 2e 63 5b 69 5d | 2e 74 72 69 61 6c 73 0d |ata.c[i]|.trials.|
|000010c0| 09 09 09 29 3b 0d 09 7d | 0d 0d 09 63 68 69 53 71 |...);..}|...chiSq|
|000010d0| 75 61 72 65 3d 32 2e 30 | 2a 28 6d 6f 6e 6f 74 6f |uare=2.0|*(monoto|
|000010e0| 6e 69 63 4c 4c 2d 77 65 | 69 62 75 6c 6c 4c 4c 29 |nicLL-we|ibullLL)|
|000010f0| 3b 0d 09 63 68 69 53 71 | 75 61 72 65 44 46 3d 6d |;..chiSq|uareDF=m|
|00001100| 6f 6e 6f 74 6f 6e 69 63 | 44 46 2d 77 65 69 62 75 |onotonic|DF-weibu|
|00001110| 6c 6c 44 46 3b 0d 09 70 | 3d 50 43 68 69 53 71 75 |llDF;..p|=PChiSqu|
|00001120| 61 72 65 28 63 68 69 53 | 71 75 61 72 65 2c 63 68 |are(chiS|quare,ch|
|00001130| 69 53 71 75 61 72 65 44 | 46 29 3b 0d 09 70 72 69 |iSquareD|F);..pri|
|00001140| 6e 74 66 28 22 5c 6e 43 | 68 69 20 73 71 75 61 72 |ntf("\nC|hi squar|
|00001150| 65 20 25 2e 31 66 20 77 | 69 74 68 20 25 64 20 64 |e %.1f w|ith %d d|
|00001160| 65 67 72 65 65 73 20 6f | 66 20 66 72 65 65 64 6f |egrees o|f freedo|
|00001170| 6d 2c 20 79 69 65 6c 64 | 69 6e 67 20 61 20 73 69 |m, yield|ing a si|
|00001180| 67 6e 69 66 69 63 61 6e | 63 65 20 6f 66 20 25 2e |gnifican|ce of %.|
|00001190| 32 66 5c 6e 22 0d 09 09 | 2c 63 68 69 53 71 75 61 |2f\n"...|,chiSqua|
|000011a0| 72 65 2c 28 69 6e 74 29 | 63 68 69 53 71 75 61 72 |re,(int)|chiSquar|
|000011b0| 65 44 46 2c 70 29 3b 0d | 7d 00 00 00 00 00 00 00 |eDF,p);.|}.......|
|000011c0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000011d0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000011e0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|000011f0| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00001200| 00 00 01 00 00 00 01 4c | 00 00 00 4c 00 00 00 32 |.......L|...L...2|
|00001210| 0d 0d 26 8c 10 14 4c ee | 18 00 ff f8 4e 5e 4e 75 |..&...L.|....N^Nu|
|00001220| 87 73 6b 69 70 5f 62 6c | 00 00 4e 56 ff fe 48 e7 |.skip_bl|..NV..H.|
|00001230| 15 54 65 73 74 50 73 79 | 63 68 6f 6d 65 74 72 69 |.TestPsy|chometri|
|00001240| 63 46 69 74 2e 63 02 00 | 00 00 50 61 72 74 53 49 |cFit.c..|..PartSI|
|00001250| 00 00 50 61 72 74 53 49 | 54 21 00 00 00 00 00 00 |..PartSI|T!......|
|00001260| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00001270| 00 00 ab bf ee 4f 00 00 | 00 00 00 00 01 7e 0c 47 |.....O..|.....~.G|
|00001280| 00 10 6d 02 5d 47 30 06 | c1 c5 d0 47 3c 00 52 94 |..m.]G0.|...G<.R.|
|00001290| 78 01 60 dc 52 47 4a 33 | 70 00 66 d8 4a 04 66 0c |x.`.RGJ3|p.f.J.f.|
|000012a0| 2f 14 48 6d f9 6a 4e ba | f8 10 50 4f 4a 6e ff fe |/.Hm.jN.|..POJn..|
|000012b0| 6c 08 48 c6 20 06 44 80 | 3c 00 48 c6 20 06 4c ee |l.H. .D.|<.H. .L.|
|000012c0| 18 f0 ff e6 4e 5e 4e 75 | 86 73 5f 61 74 6f 69 00 |....N^Nu|.s_atoi.|
|000012d0| 00 00 4e 56 ff f8 48 e7 | 01 18 26 6e 00 0c 28 6e |..NV..H.|..&n..(n|
|000012e0| 00 08 60 2e 10 1c 48 80 | 48 c0 2f 00 4e ad 01 d2 |..`...H.|H./.N...|
|000012f0| 2d 40 ff fc 10 1b 48 80 | 48 c0 2f 00 4e ad 01 d2 |-@....H.|H./.N...|
|00001300| 00 00 00 48 00 09 4d 6f | 6e 61 63 6f 00 00 00 00 |...H..Mo|naco....|
|00001310| 00 00 00 00 00 00 00 00 | 00 00 00 00 00 00 00 00 |........|........|
|00001320| 00 00 00 00 00 00 00 03 | 00 04 00 2d 00 04 01 8c |........|...-....|
|00001330| 02 7c 00 2d 00 04 01 8c | 02 7c aa 91 5b 69 00 00 |.|.-....|.|..[i..|
|00001340| 04 0b 00 00 04 0b 00 00 | 01 fd 01 00 00 00 01 00 |........|........|
|00001350| 00 00 01 4c 00 00 00 4c | 00 00 00 32 00 be 2e 7c |...L...L|...2...||
|00001360| 12 bc 00 00 00 1c 00 32 | 00 00 4d 50 53 52 00 00 |.......2|..MPSR..|
|00001370| 00 0a 03 ed ff ff 00 00 | 00 00 00 00 00 00 00 00 |........|........|
+--------+-------------------------+-------------------------+--------+--------+