;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;; ;;; Model and Batch Run Code for Psychomotor Vigilance Task (PVT) ;;; ;;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;; ;;; Model and results described in: ;;; ;;; Gunzelmann, G., Gross, J. B., Gluck, K. A., & Dinges, D. F. (in press). ;;; Sleep deprivation and sustained attention performance: Integrating ;;; mathematical and cognitive modeling. Cognitive Science. ;;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;; ;;; Author : Glenn Gunzelmann ;;; Copyright : Work of the US Government ;;; Availability: ;;; Address : AFRL/RHAC ;;; : 6030 South Kent St. ;;; : Mesa, AZ 85212-6061 ;;; : glenn.gunzelmann@us.af.mil ;;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;; ;;; Model Running Instructions: ;;; ;;; 1. Software Requirements ;;; A. Model was developed in ACL 6.2 for Windows ;;; i. Portability to other lisp environments/versions is unknown ;;; B. Model was developed using ACT-R 5 ;;; i. ACT-R version 5.1 ;;; ii. PM version 2.2 ;;; ;;; 2. Function Calls ;;; A. To run the model using current ACT-R parameter values ;;; 1. Call (pvt DURATION :SLEEP-THRESHOLD :VISIBLE :GDECREMENT) ;;; a. DURATION - Length of the session in seconds ;;; b. :SLEEP-THRESHOLD - Number of milliseconds allotted for the model to respond ;;; i. After this time, the model will be halted, a sleep attack recorded, ;;; and the next trial initiated ;;; ii. Default is 30000 (30 s) ;;; c. :VISIBLE - Determines whether real or virtual ACT-R windows are used ;;; i. T/NIL, default is T ;;; d. :GDECREMENT - Decrement assessed to G when a micro-lapse occurs ;;; i. Default is 0 ;;; ii. Value used in model fits is 0.035 ;;; B. To run the model through multiple sessions or across multiple parameter value combinations ;;; 1. Call (pvt-batch N :SUMMARY-FILE :SUMMARY-2 :IND-DATA :VISIBLE :FALSE-START-THRESHOLD ;;; :INTERVAL :LAPSE-THRESHOLD :SLEEP-THRESHOLD ;;; :SESSION-LENGTH :GDECREMENT :PARAMS ;;; :LAPSE-INTERVAL :HUMAN-DATA ;;; a. N - Number of iterations of the model to run ;;; b. :SUMMARY-FILE - Name of the file to generate to store summary-level data ;;; i. Default is "batch-data.txt" ;;; ii. Aggregation is 1 line of data per run ;;; c. :SUMMARY-2 - Name of the file to generate to store more reduced summary data ;;; i. Default is "stat-data.txt" ;;; ii. Aggregation is 1 line of data per parameter combination (average over runs) ;;; d. :IND-DATA - Flag to indicate whether or not to record individual model files ;;; i. File contains each individual trial, including the delay and reaction time ;;; ii. T/NIL; Default is NIL ;;; e. :VISIBLE - Whether to use a real or virtual ACT-R window ;;; i. Default is NIL ;;; f. :FALSE-START-THRESHOLD - Lower bound on alert reaction time ;;; i. Default is 150 ms ;;; g. :INTERVAL - For alert responses, how large to make the reaction time bins for computing proportions ;;; i. Default is 10 ms ;;; h. :LAPSE-THRESHOLD - Upper bound on alert reaction time ;;; i. Default is 500 ms ;;; i. :SLEEP-THRESHOLD - Number of milliseconds allotted for the model to respond ;;; i. After this time, the model will be halted, a sleep attack recorded, ;;; and the next trial initiated ;;; ii. Default is 30000 (30 s) ;;; j. :SESSION-LENGTH - Number of simulated seconds for each session ;;; i. Default is 600 ;;; k. :GDECREMENT - Decrement assessed to G when a micro-lapse occurs ;;; i. Default is 0 ;;; ii. Value used in model fits is 0.035 ;;; l. :PARAMS - List of parameters and ranges to run ;;; i. Format is '((NAME LOCATION MIN MAX STEP) (N2 L2 MN2 MX2 S2) ...) ;;; a. NAME - The parameter name (e.g., :g, :ut, or *gdecrement*) ;;; b. LOCATION - ACTR, RPM, or NONE ;;; -Determines how parameter value is set (e.g., sgp, rpm-set-params, or setf) ;;; c. MIN - Minimum value to run ;;; d. MAX - Maximum value to run ;;; e. STEP - Step size between values of the parameter ;;; m. :LAPSE-INTERVAL and :HUMAN-DATA are depricated ;;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;; Global parameters for running experiment ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;Main Experiment Window (defparameter *w* nil) ;;Variable to hold response information (defparameter *response* nil) ;;Parameters to support fatigue mechanisms (defparameter *actr-parameter-list* nil) (defparameter *rpm-parameter-list* nil) (defparameter *gdecrement* 0) ;;Hash table for collecting data across runs (defparameter *all* (make-hash-table :test #'equalp)) ;;Path information for saving data files (defparameter *path* (make-pathname :device (pathname-device *load-pathname*) :directory (pathname-directory *load-pathname*))) (defparameter *ind-path* (format nil "~aind-data/" *path*)) (defparameter *all-path* (format nil "~aall-data/" *path*)) ;;Create data directories if they do not exist (if (null (directory *ind-path*)) (make-directory *ind-path*)) (if (null (directory *all-path*)) (make-directory *all-path*)) ;;Structure for passing parameter space information into model for multiple runs (defstruct (parameter (:type list) (:constructor mk-parameter (keyword system min max interval))) keyword system min max interval) ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;; ;;; PVT Data from Penn - 88hr total sleep deprivation protocol ;;; -Averaged over each day ;;; -Proportion of responses ;;; -Data match what is presented in Figure 1 of paper ;;; ;;; *Data presented for information. Use requires permission from ;;; Dr. David Dinges, University of Pennsylvania (dinges@mail.med.upenn.edu) ;;; ;;; *Data from study in: ;;; -Doran, S. M., Van Dongen, H. P., & Dinges, D. F (2001). ;;; Sustained attention performance during sleep deprivation: ;;; Evidence of state instability. Archives of Italian Biology: ;;; Neuroscience 139(3), 253-267. ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; (defparameter *baseline* '(0.037126471 0.001006137 0.004527618 0.009960761 0.021531341 0.034812355 0.037830768 0.041452862 0.034007445 0.045074957 0.046886005 0.045175571 0.047992756 0.046886005 0.043968206 0.04165409 0.037528926 0.042157159 0.037528926 0.029681054 0.028171848 0.026964483 0.025757118 0.022939934 0.017003723 0.014689607 0.01559513 0.012979173 0.01398531 0.012677332 0.009759533 0.008753396 0.009357078 0.00563437 0.006338666 0.005533756 0.076969514 0.000100614)) (defparameter *day1* '(0.054139903 0.000282346 0.001764664 0.006564551 0.014823181 0.019270135 0.024705301 0.025622927 0.028305216 0.029999294 0.034587421 0.037340298 0.03854027 0.038893203 0.037693231 0.036140326 0.033599209 0.034799181 0.036422672 0.031905132 0.027175831 0.027740524 0.02287005 0.02308181 0.021105386 0.019481895 0.017081951 0.013058516 0.013270276 0.010023294 0.011011506 0.009246841 0.009317428 0.009388014 0.007058657 0.007482177 0.18380744 0.002399944)) (defparameter *day2* '(0.078253159 0.000591148 0.001330082 0.00406414 0.011084017 0.014852583 0.01736496 0.01995123 0.02187246 0.023645903 0.02660164 0.030074632 0.027857829 0.031995862 0.035025493 0.034064878 0.030296313 0.034212665 0.026971108 0.028670657 0.030222419 0.025493239 0.022906968 0.021798566 0.020616271 0.016921599 0.015517624 0.014039755 0.013153033 0.010640656 0.01093623 0.010345082 0.010418976 0.010418976 0.009901722 0.008276066 0.241631567 0.007980492)) (defparameter *day3* '(0.101803008 0.000157468 0.001181009 0.004802771 0.011652626 0.016770333 0.017793874 0.021179435 0.019211086 0.019998425 0.021179435 0.023305252 0.025509802 0.028737895 0.031336115 0.027399417 0.027793087 0.027793087 0.025431068 0.025352334 0.024958665 0.02267538 0.019604756 0.01889615 0.019053618 0.016061727 0.01511692 0.013069837 0.011652626 0.013384773 0.010314148 0.010865286 0.008503267 0.009133139 0.008109598 0.006219983 0.271553421 0.022439178)) ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;; ;;; PVT Data from Penn - 88hr total sleep deprivation protocol ;;; -Averaged by dependent measure for each session ;;; -Proportion of responses classified as false starts ;;; -Median of the alert responses ;;; -Proportion of responses classified as lapses ;;; -Proportion of responses classified as sleep attacks ;;; -Data match what is presented in Figure 6 of paper ;;; ;;; *Data presented for information. Use requires permission from ;;; Dr. David Dinges, University of Pennsylvania (dinges@mail.med.upenn.edu) ;;; ;;; *Data from study in: ;;; -Doran, S. M., Van Dongen, H. P., & Dinges, D. F (2001). ;;; Sustained attention performance during sleep deprivation: ;;; Evidence of state instability. Archives of Italian Biology: ;;; Neuroscience 139(3), 253-267. ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; (defparameter *false-starts* '(0.0234 0.0265 0.0300 0.0446 0.0316 0.0342 0.0544 0.0514 0.0321 0.0449 0.0570 0.0722 0.0658 0.0470 0.0567 0.0619 0.0579 0.0498 0.0515 0.0554 0.0568 0.0571 0.0849 0.1106 0.0643 0.0797 0.0967 0.0639 0.0723 0.0728 0.0864 0.0959 0.0780 0.0855 0.0957 0.1431 0.0947 0.0727 0.0589 0.1498 0.0967 0.1099 0.1552 0.0737)) (defparameter *median-rt* '(285.0 284.0 277.5 286.0 288.0 292.0 280.0 291.0 287.0 311.0 315.0 306.0 317.5 301.0 295.0 306.0 291.0 298.0 304.0 294.0 302.0 319.0 319.0 310.5 319.0 314.0 310.0 305.0 310.5 317.5 316.0 285.5 300.0 315.0 330.0 329.0 328.5 321.5 305.0 306.0 301.0 302.0 318.0 295.5)) (defparameter *lapses* '(0.0679 0.1014 0.0657 0.1098 0.0720 0.0766 0.0528 0.0711 0.1098 0.1629 0.2341 0.2263 0.2733 0.2597 0.1802 0.2122 0.1273 0.1762 0.1286 0.1369 0.1753 0.3161 0.2528 0.3271 0.2554 0.2647 0.2278 0.2233 0.2197 0.2401 0.2147 0.1993 0.2400 0.2875 0.3310 0.2930 0.3426 0.3293 0.2605 0.2326 0.2388 0.2085 0.2850 0.2394)) (defparameter *sleep-attacks* '(0.0000 0.0000 0.0000 0.0008 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0008 0.0078 0.0101 0.0046 0.0000 0.0017 0.0008 0.0000 0.0000 0.0042 0.0000 0.0045 0.0138 0.0139 0.0100 0.0082 0.0136 0.0035 0.0118 0.0072 0.0077 0.0034 0.0035 0.0118 0.0279 0.0319 0.0298 0.0540 0.0250 0.0226 0.0186 0.0138 0.0245 0.0147)) ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;; ;;; Alertness predictions from biomathematical models ;;; -Predictions are made at 2hr intervals over the sleep ;;; deprivation period, beginning at 0800 on the baseline day ;;; ;;; *CNPA: ;;; -Jewett, M. E., & Kronauer, R. E. (1999). Interactive ;;; mathematical models of subjective alertness and alertness ;;; in humans. Journal of Biological Rhythms, 14, 588-597. ;;; ;;; *SAFTE: ;;; -Hursh, S. R., Redmond, D. P., Johnson, M. L., Thorne, D. R., ;;; Belenky, G., Balkin, T. J., Storm, W. F., Miller, J. C., & ;;; Eddy, D. R. (2004). Fatigue models for applied research in ;;; warfighting. Aviation, Space, and Environmental Medicine, ;;; 75(3), A44-60. ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; (defparameter *cnpa* '(0.8890 0.9022 0.9204 0.9300 0.9463 0.9503 0.9348 0.8897 0.8118 0.7141 0.6189 0.5552 0.5212 0.5189 0.5458 0.5849 0.6152 0.6101 0.5748 0.5097 0.4182 0.3171 0.2274 0.1761 0.1591 0.1760 0.2233 0.2837 0.3363 0.3541 0.3411 0.2982 0.2277 0.1464 0.0745 0.0392 0.0375 0.0682 0.1285 0.2012 0.2659 0.2952 0.2921 0.2581)) (defparameter *safte* '(99.11 99.79 98.22 96.50 97.02 99.32 99.87 95.27 85.83 75.86 70.17 70.10 72.67 73.82 72.38 70.72 71.65 74.67 75.70 70.78 60.21 48.97 42.67 42.87 46.16 47.86 46.55 44.94 46.27 50.01 51.53 46.29 34.59 22.08 15.17 15.64 19.66 21.89 20.71 19.15 20.89 25.36 27.37 21.79)) ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;; ;;; ACT-R Model Run Results ;;; -Used to generate model data presented in Figures 4 and 6 ;;; -Parameter values for G and Tu are presented for the session fits ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; #| Model Fits - Unconstrained Parameter Search; Data Averaged over Whole Days Day EGS G UT GDEC Median_RT StdDev F_Starts LT160 LT170 LT180 LT190 LT200 LT210 LT220 LT230 LT240 LT250 LT260 LT270 LT280 LT290 LT300 LT310 LT320 LT330 LT340 LT350 LT360 LT370 LT380 LT390 LT400 LT410 LT420 LT430 LT440 LT450 LT460 LT470 LT480 LT490 LT500 Lapses S_Attacks Baseline 0.25 1.9799999 1.84 0.035 292 72.982994 0.039920654 2.48E-04 7.44E-04 0.002851475 0.007934541 0.015497149 0.022315893 0.03347384 0.035581455 0.048723035 0.04760724 0.053186215 0.0538061 0.054178033 0.054054055 0.046243493 0.04351599 0.038060997 0.03446566 0.034961566 0.034961566 0.025539301 0.023431689 0.020456236 0.021696009 0.017976692 0.01673692 0.013761468 0.012521696 0.013141582 0.013141582 0.011281924 0.008678403 0.007934541 0.008678403 0.008058517 0.07364245 9.92E-04 Day1TSD 0.25 1.7999998 1.78 0.035 313 77.97513 0.059145954 3.82E-04 0.002166985 0.002166985 0.005863607 0.007903123 0.014786488 0.021924794 0.024219248 0.031102613 0.029955385 0.034544297 0.0416826 0.03709369 0.041045252 0.034161888 0.033779476 0.03301466 0.03212237 0.028808158 0.02817081 0.030465264 0.023964308 0.021414913 0.024474187 0.017845761 0.018100701 0.019247929 0.015933717 0.017080944 0.014531549 0.01427661 0.012492033 0.012746973 0.008922881 0.010325048 0.18432122 0.009815169 Day2TSD 0.25 1.66 1.7 0.035 325 79.063416 0.08484118 9.38E-04 0.001072242 0.002546576 0.003618818 0.008845999 0.00911406 0.014877362 0.020104544 0.021042755 0.024929633 0.026135907 0.024527544 0.029218603 0.033641603 0.033373542 0.02774427 0.032971453 0.026537998 0.026403967 0.02707412 0.025867846 0.024929633 0.019970514 0.021846937 0.019032301 0.019702453 0.01863021 0.016083635 0.017557967 0.014609301 0.013134969 0.013939151 0.012464817 0.011794666 0.010856453 0.24299692 0.017021846 Day3TSD 0.25 1.5799999 1.64 0.035 327 80.37472 0.10919462 4.07E-04 6.79E-04 0.002444656 0.00421024 0.007469781 0.008284667 0.014396306 0.015890263 0.01874236 0.021187017 0.025261443 0.02621214 0.02838517 0.030422382 0.024854 0.028792612 0.024039114 0.023903301 0.02159446 0.024582371 0.02132283 0.022001902 0.018334918 0.021866087 0.01901399 0.018878175 0.015754448 0.016433518 0.017520033 0.014803749 0.01358142 0.013445606 0.011951651 0.011815836 0.010729322 0.27285075 0.01874236 Model Fits - Parameters Constrained Using CNPA Session EGS G UT GDEC Median_RT StdDev F_Starts LT160 LT170 LT180 LT190 LT200 LT210 LT220 LT230 LT240 LT250 LT260 LT270 LT280 LT290 LT300 LT310 LT320 LT330 LT340 LT350 LT360 LT370 LT380 LT390 LT400 LT410 LT420 LT430 LT440 LT450 LT460 LT470 LT480 LT490 LT500 Lapses S_Attacks 1 0.25 2.080 1.917 0.035 290 70.7113 0.0242 0.0002 0.0015 0.0044 0.0087 0.0147 0.0224 0.0330 0.0405 0.0469 0.0515 0.0584 0.0582 0.0575 0.0546 0.0494 0.0457 0.0434 0.0381 0.0332 0.0307 0.0296 0.0242 0.0229 0.0208 0.0173 0.0157 0.0157 0.0126 0.0124 0.0097 0.0092 0.0079 0.0072 0.0067 0.0069 0.0633 0.0007 2 0.25 2.088 1.922 0.035 290 70.6960 0.0233 0.0002 0.0012 0.0038 0.0095 0.0155 0.0228 0.0326 0.0418 0.0468 0.0528 0.0592 0.0581 0.0594 0.0559 0.0506 0.0485 0.0411 0.0357 0.0348 0.0310 0.0288 0.0247 0.0229 0.0181 0.0199 0.0162 0.0139 0.0123 0.0117 0.0109 0.0098 0.0085 0.0072 0.0077 0.0058 0.0564 0.0007 3 0.25 2.100 1.929 0.035 290 70.1722 0.0230 0.0001 0.0013 0.0040 0.0091 0.0139 0.0245 0.0332 0.0405 0.0483 0.0531 0.0575 0.0587 0.0564 0.0558 0.0534 0.0456 0.0414 0.0381 0.0365 0.0316 0.0281 0.0250 0.0230 0.0208 0.0176 0.0172 0.0136 0.0141 0.0118 0.0100 0.0095 0.0076 0.0072 0.0065 0.0057 0.0557 0.0006 4 0.25 2.107 1.932 0.035 289 69.6301 0.0232 0.0002 0.0013 0.0043 0.0086 0.0169 0.0249 0.0333 0.0409 0.0464 0.0578 0.0581 0.0591 0.0575 0.0576 0.0519 0.0476 0.0413 0.0382 0.0339 0.0319 0.0284 0.0253 0.0219 0.0192 0.0188 0.0155 0.0146 0.0125 0.0106 0.0094 0.0093 0.0080 0.0066 0.0057 0.0058 0.0532 0.0004 5 0.25 2.117 1.939 0.035 288 69.6683 0.0216 0.0003 0.0019 0.0043 0.0088 0.0157 0.0252 0.0340 0.0446 0.0505 0.0560 0.0568 0.0578 0.0600 0.0527 0.0524 0.0488 0.0413 0.0382 0.0349 0.0327 0.0285 0.0246 0.0206 0.0182 0.0175 0.0160 0.0152 0.0131 0.0107 0.0090 0.0097 0.0070 0.0066 0.0057 0.0058 0.0526 0.0004 6 0.25 2.120 1.940 0.035 287 69.6806 0.0205 0.0002 0.0018 0.0043 0.0092 0.0173 0.0238 0.0333 0.0417 0.0508 0.0569 0.0594 0.0621 0.0604 0.0572 0.0505 0.0457 0.0428 0.0356 0.0348 0.0309 0.0284 0.0246 0.0215 0.0201 0.0169 0.0161 0.0138 0.0115 0.0110 0.0101 0.0093 0.0074 0.0066 0.0063 0.0056 0.0508 0.0007 7 0.25 2.110 1.934 0.035 288 69.9684 0.0228 0.0000 0.0014 0.0041 0.0091 0.0144 0.0259 0.0339 0.0412 0.0491 0.0558 0.0578 0.0612 0.0578 0.0562 0.0503 0.0486 0.0418 0.0381 0.0341 0.0316 0.0276 0.0242 0.0219 0.0193 0.0159 0.0166 0.0139 0.0136 0.0113 0.0106 0.0088 0.0083 0.0071 0.0060 0.0054 0.0536 0.0004 8 0.25 2.080 1.917 0.035 291 70.6984 0.0255 0.0000 0.0014 0.0046 0.0093 0.0150 0.0241 0.0309 0.0384 0.0496 0.0536 0.0565 0.0592 0.0569 0.0533 0.0502 0.0442 0.0439 0.0391 0.0341 0.0310 0.0285 0.0262 0.0222 0.0206 0.0189 0.0169 0.0154 0.0137 0.0106 0.0105 0.0083 0.0082 0.0074 0.0070 0.0063 0.0576 0.0008 9 0.25 2.029 1.888 0.035 294 71.6000 0.0291 0.0002 0.0010 0.0041 0.0091 0.0131 0.0211 0.0295 0.0388 0.0451 0.0499 0.0519 0.0557 0.0548 0.0528 0.0506 0.0429 0.0423 0.0382 0.0352 0.0319 0.0300 0.0256 0.0231 0.0220 0.0196 0.0168 0.0155 0.0135 0.0116 0.0115 0.0108 0.0083 0.0081 0.0068 0.0072 0.0711 0.0012 10 0.25 1.965 1.852 0.035 297 73.9131 0.0375 0.0001 0.0010 0.0034 0.0074 0.0124 0.0193 0.0280 0.0354 0.0411 0.0474 0.0522 0.0517 0.0506 0.0492 0.0467 0.0422 0.0382 0.0337 0.0347 0.0324 0.0272 0.0270 0.0232 0.0223 0.0191 0.0183 0.0164 0.0135 0.0146 0.0124 0.0116 0.0095 0.0099 0.0095 0.0080 0.0907 0.0024 11 0.25 1.902 1.817 0.035 303 74.9541 0.0427 0.0003 0.0011 0.0033 0.0060 0.0124 0.0185 0.0234 0.0317 0.0380 0.0414 0.0454 0.0472 0.0471 0.0474 0.0436 0.0392 0.0370 0.0356 0.0302 0.0320 0.0307 0.0287 0.0244 0.0199 0.0220 0.0190 0.0176 0.0162 0.0139 0.0128 0.0116 0.0113 0.0098 0.0101 0.0084 0.1173 0.0028 12 0.25 1.860 1.793 0.035 306 76.2000 0.0460 0.0003 0.0010 0.0029 0.0056 0.0110 0.0186 0.0243 0.0276 0.0326 0.0378 0.0440 0.0436 0.0482 0.0425 0.0438 0.0391 0.0363 0.0343 0.0327 0.0295 0.0282 0.0260 0.0265 0.0218 0.0203 0.0182 0.0178 0.0166 0.0155 0.0142 0.0141 0.0129 0.0100 0.0105 0.0093 0.1319 0.0044 13 0.25 1.838 1.780 0.035 308 76.2455 0.0494 0.0002 0.0012 0.0031 0.0055 0.0105 0.0179 0.0205 0.0273 0.0320 0.0396 0.0386 0.0432 0.0448 0.0450 0.0413 0.0385 0.0372 0.0355 0.0345 0.0290 0.0295 0.0243 0.0239 0.0211 0.0224 0.0185 0.0187 0.0176 0.0140 0.0148 0.0133 0.0131 0.0118 0.0095 0.0099 0.1377 0.0051 14 0.25 1.836 1.779 0.035 308 75.9123 0.0512 0.0002 0.0012 0.0034 0.0067 0.0095 0.0161 0.0238 0.0273 0.0327 0.0367 0.0392 0.0462 0.0448 0.0424 0.0405 0.0396 0.0362 0.0348 0.0320 0.0331 0.0280 0.0254 0.0256 0.0211 0.0201 0.0203 0.0188 0.0171 0.0158 0.0147 0.0113 0.0117 0.0108 0.0104 0.0088 0.1370 0.0056 15 0.25 1.854 1.789 0.035 308 75.8956 0.0463 0.0002 0.0010 0.0034 0.0060 0.0108 0.0162 0.0216 0.0271 0.0333 0.0382 0.0416 0.0435 0.0450 0.0428 0.0430 0.0387 0.0380 0.0359 0.0325 0.0310 0.0290 0.0276 0.0228 0.0217 0.0212 0.0203 0.0175 0.0172 0.0150 0.0143 0.0131 0.0124 0.0111 0.0101 0.0091 0.1377 0.0040 16 0.25 1.880 1.804 0.035 304 74.9755 0.0433 0.0003 0.0011 0.0030 0.0063 0.0114 0.0163 0.0231 0.0303 0.0352 0.0405 0.0454 0.0445 0.0464 0.0480 0.0431 0.0396 0.0377 0.0343 0.0326 0.0324 0.0296 0.0258 0.0255 0.0229 0.0205 0.0175 0.0176 0.0159 0.0158 0.0135 0.0113 0.0112 0.0118 0.0092 0.0080 0.1258 0.0034 17 0.25 1.900 1.815 0.035 302 74.7768 0.0418 0.0001 0.0012 0.0032 0.0070 0.0122 0.0193 0.0250 0.0305 0.0362 0.0435 0.0448 0.0483 0.0483 0.0475 0.0441 0.0400 0.0384 0.0360 0.0327 0.0316 0.0288 0.0254 0.0238 0.0224 0.0207 0.0173 0.0161 0.0160 0.0142 0.0122 0.0123 0.0113 0.0102 0.0095 0.0081 0.1166 0.0031 18 0.25 1.896 1.813 0.035 303 75.2531 0.0438 0.0001 0.0012 0.0036 0.0066 0.0134 0.0164 0.0253 0.0307 0.0371 0.0410 0.0461 0.0454 0.0501 0.0446 0.0420 0.0421 0.0375 0.0354 0.0326 0.0302 0.0298 0.0255 0.0244 0.0220 0.0193 0.0198 0.0167 0.0152 0.0155 0.0128 0.0126 0.0126 0.0098 0.0096 0.0082 0.1183 0.0031 19 0.25 1.873 1.800 0.035 305 75.6424 0.0458 0.0003 0.0009 0.0033 0.0072 0.0123 0.0172 0.0244 0.0293 0.0349 0.0420 0.0436 0.0454 0.0446 0.0459 0.0400 0.0399 0.0389 0.0336 0.0332 0.0324 0.0282 0.0258 0.0244 0.0236 0.0202 0.0199 0.0180 0.0151 0.0145 0.0137 0.0116 0.0118 0.0106 0.0098 0.0091 0.1244 0.0044 20 0.25 1.830 1.776 0.035 308 76.5728 0.0481 0.0004 0.0012 0.0031 0.0058 0.0108 0.0166 0.0236 0.0282 0.0331 0.0391 0.0395 0.0428 0.0438 0.0410 0.0405 0.0390 0.0362 0.0326 0.0310 0.0320 0.0283 0.0268 0.0224 0.0225 0.0224 0.0194 0.0168 0.0162 0.0166 0.0145 0.0128 0.0122 0.0112 0.0097 0.0092 0.1457 0.0051 21 0.25 1.770 1.742 0.035 311 77.3759 0.0599 0.0004 0.0012 0.0029 0.0061 0.0092 0.0140 0.0209 0.0251 0.0294 0.0371 0.0372 0.0382 0.0419 0.0413 0.0379 0.0356 0.0325 0.0319 0.0308 0.0283 0.0285 0.0264 0.0235 0.0232 0.0197 0.0196 0.0168 0.0177 0.0157 0.0151 0.0133 0.0125 0.0110 0.0107 0.0104 0.1664 0.0077 22 0.25 1.704 1.704 0.035 317 78.2746 0.0695 0.0004 0.0010 0.0023 0.0047 0.0082 0.0137 0.0172 0.0232 0.0253 0.0292 0.0332 0.0346 0.0391 0.0363 0.0341 0.0323 0.0320 0.0291 0.0287 0.0274 0.0261 0.0254 0.0238 0.0227 0.0203 0.0198 0.0183 0.0175 0.0150 0.0150 0.0140 0.0129 0.0123 0.0112 0.0112 0.2026 0.0106 23 0.25 1.645 1.671 0.035 321 79.2536 0.0838 0.0005 0.0010 0.0026 0.0048 0.0074 0.0119 0.0155 0.0192 0.0214 0.0290 0.0304 0.0291 0.0327 0.0334 0.0332 0.0311 0.0284 0.0278 0.0288 0.0262 0.0235 0.0241 0.0224 0.0208 0.0202 0.0168 0.0184 0.0171 0.0151 0.0163 0.0127 0.0136 0.0138 0.0103 0.0124 0.2306 0.0137 24 0.25 1.611 1.652 0.035 325 79.4712 0.0892 0.0003 0.0014 0.0022 0.0043 0.0068 0.0097 0.0138 0.0195 0.0236 0.0257 0.0277 0.0276 0.0316 0.0304 0.0273 0.0288 0.0281 0.0281 0.0282 0.0250 0.0257 0.0234 0.0209 0.0216 0.0208 0.0194 0.0170 0.0163 0.0140 0.0147 0.0133 0.0139 0.0129 0.0128 0.0110 0.2487 0.0140 25 0.25 1.600 1.645 0.035 325 80.1092 0.0917 0.0006 0.0012 0.0028 0.0042 0.0080 0.0111 0.0116 0.0174 0.0218 0.0256 0.0259 0.0282 0.0308 0.0310 0.0289 0.0301 0.0260 0.0275 0.0259 0.0251 0.0236 0.0219 0.0219 0.0202 0.0190 0.0181 0.0192 0.0169 0.0161 0.0149 0.0132 0.0133 0.0120 0.0141 0.0116 0.2540 0.0144 26 0.25 1.611 1.652 0.035 323 79.8000 0.0906 0.0006 0.0010 0.0029 0.0040 0.0080 0.0104 0.0149 0.0165 0.0205 0.0261 0.0285 0.0281 0.0321 0.0314 0.0310 0.0286 0.0291 0.0249 0.0247 0.0261 0.0228 0.0215 0.0228 0.0208 0.0193 0.0185 0.0178 0.0166 0.0167 0.0153 0.0135 0.0137 0.0140 0.0105 0.0118 0.2491 0.0155 27 0.25 1.642 1.669 0.035 324 78.7605 0.0843 0.0006 0.0009 0.0022 0.0051 0.0079 0.0120 0.0147 0.0190 0.0225 0.0251 0.0286 0.0293 0.0313 0.0307 0.0308 0.0303 0.0295 0.0289 0.0304 0.0277 0.0255 0.0235 0.0232 0.0222 0.0202 0.0173 0.0178 0.0186 0.0151 0.0146 0.0149 0.0123 0.0122 0.0124 0.0101 0.2340 0.0143 28 0.25 1.682 1.692 0.035 320 78.7444 0.0743 0.0003 0.0015 0.0024 0.0052 0.0085 0.0129 0.0157 0.0226 0.0243 0.0269 0.0315 0.0343 0.0362 0.0330 0.0324 0.0313 0.0317 0.0311 0.0276 0.0280 0.0258 0.0253 0.0231 0.0227 0.0192 0.0203 0.0175 0.0164 0.0157 0.0156 0.0150 0.0144 0.0115 0.0115 0.0104 0.2117 0.0122 29 0.25 1.716 1.711 0.035 316 78.7254 0.0681 0.0002 0.0009 0.0029 0.0051 0.0087 0.0149 0.0175 0.0249 0.0277 0.0298 0.0321 0.0360 0.0374 0.0372 0.0342 0.0332 0.0310 0.0319 0.0288 0.0280 0.0266 0.0251 0.0211 0.0230 0.0224 0.0161 0.0179 0.0176 0.0161 0.0155 0.0146 0.0127 0.0133 0.0117 0.0106 0.1951 0.0098 30 0.25 1.728 1.718 0.035 317 78.1414 0.0668 0.0004 0.0009 0.0025 0.0058 0.0090 0.0138 0.0157 0.0220 0.0277 0.0304 0.0357 0.0357 0.0362 0.0385 0.0348 0.0327 0.0328 0.0273 0.0312 0.0280 0.0280 0.0268 0.0254 0.0222 0.0198 0.0207 0.0195 0.0170 0.0147 0.0146 0.0134 0.0134 0.0127 0.0121 0.0107 0.1922 0.0089 31 0.25 1.720 1.713 0.035 317 77.8407 0.0652 0.0004 0.0013 0.0027 0.0051 0.0089 0.0133 0.0178 0.0208 0.0280 0.0309 0.0342 0.0340 0.0360 0.0360 0.0347 0.0350 0.0345 0.0318 0.0279 0.0284 0.0266 0.0257 0.0245 0.0229 0.0221 0.0192 0.0195 0.0184 0.0154 0.0155 0.0120 0.0126 0.0127 0.0104 0.0108 0.1953 0.0095 32 0.25 1.691 1.697 0.035 319 79.2006 0.0763 0.0006 0.0010 0.0023 0.0053 0.0075 0.0136 0.0170 0.0211 0.0243 0.0279 0.0329 0.0349 0.0337 0.0335 0.0322 0.0330 0.0333 0.0297 0.0272 0.0252 0.0262 0.0252 0.0252 0.0200 0.0201 0.0185 0.0187 0.0177 0.0176 0.0157 0.0144 0.0132 0.0124 0.0119 0.0115 0.2088 0.0102 33 0.25 1.645 1.671 0.035 322 79.0401 0.0828 0.0003 0.0014 0.0031 0.0045 0.0067 0.0116 0.0159 0.0203 0.0230 0.0273 0.0293 0.0316 0.0303 0.0335 0.0305 0.0304 0.0310 0.0288 0.0263 0.0253 0.0265 0.0243 0.0217 0.0213 0.0204 0.0190 0.0184 0.0177 0.0166 0.0142 0.0134 0.0139 0.0122 0.0102 0.0118 0.2308 0.0137 34 0.25 1.592 1.641 0.035 324 80.1357 0.0952 0.0003 0.0011 0.0021 0.0049 0.0078 0.0105 0.0148 0.0167 0.0210 0.0231 0.0270 0.0277 0.0285 0.0307 0.0292 0.0298 0.0274 0.0257 0.0265 0.0246 0.0214 0.0218 0.0222 0.0192 0.0175 0.0177 0.0164 0.0170 0.0167 0.0147 0.0151 0.0128 0.0118 0.0124 0.0110 0.2614 0.0165 35 0.25 1.544 1.614 0.035 328 80.3938 0.1054 0.0004 0.0012 0.0024 0.0045 0.0064 0.0100 0.0121 0.0156 0.0191 0.0209 0.0240 0.0268 0.0267 0.0291 0.0261 0.0253 0.0247 0.0252 0.0222 0.0253 0.0217 0.0214 0.0197 0.0205 0.0193 0.0180 0.0164 0.0162 0.0150 0.0142 0.0127 0.0153 0.0129 0.0110 0.0112 0.2854 0.0157 36 0.25 1.521 1.601 0.035 332 80.7744 0.1086 0.0006 0.0012 0.0024 0.0035 0.0053 0.0100 0.0118 0.0143 0.0190 0.0199 0.0227 0.0217 0.0256 0.0261 0.0243 0.0242 0.0244 0.0227 0.0253 0.0222 0.0220 0.0213 0.0199 0.0190 0.0185 0.0179 0.0184 0.0163 0.0146 0.0143 0.0142 0.0142 0.0133 0.0109 0.0120 0.2988 0.0187 37 0.25 1.520 1.600 0.035 331 80.1629 0.1133 0.0003 0.0006 0.0019 0.0034 0.0052 0.0096 0.0122 0.0161 0.0179 0.0208 0.0228 0.0239 0.0257 0.0256 0.0251 0.0253 0.0231 0.0238 0.0209 0.0220 0.0210 0.0216 0.0228 0.0192 0.0190 0.0164 0.0167 0.0157 0.0165 0.0146 0.0140 0.0130 0.0137 0.0114 0.0092 0.2979 0.0179 38 0.25 1.540 1.611 0.035 329 80.2284 0.1054 0.0005 0.0013 0.0021 0.0041 0.0067 0.0092 0.0124 0.0171 0.0194 0.0205 0.0234 0.0273 0.0278 0.0250 0.0245 0.0260 0.0256 0.0246 0.0255 0.0216 0.0220 0.0209 0.0208 0.0194 0.0198 0.0176 0.0171 0.0185 0.0149 0.0141 0.0146 0.0114 0.0128 0.0109 0.0111 0.2867 0.0176 39 0.25 1.580 1.634 0.035 326 78.8715 0.0985 0.0008 0.0011 0.0027 0.0035 0.0065 0.0097 0.0139 0.0155 0.0206 0.0236 0.0272 0.0293 0.0290 0.0263 0.0283 0.0277 0.0281 0.0276 0.0259 0.0236 0.0263 0.0218 0.0208 0.0214 0.0210 0.0187 0.0183 0.0174 0.0146 0.0127 0.0145 0.0113 0.0117 0.0108 0.0115 0.2597 0.0183 40 0.25 1.628 1.661 0.035 323 79.1340 0.0842 0.0004 0.0009 0.0025 0.0045 0.0075 0.0101 0.0139 0.0180 0.0240 0.0265 0.0301 0.0308 0.0321 0.0306 0.0307 0.0289 0.0296 0.0291 0.0275 0.0251 0.0243 0.0236 0.0238 0.0205 0.0204 0.0192 0.0194 0.0153 0.0158 0.0156 0.0143 0.0134 0.0117 0.0130 0.0118 0.2378 0.0133 41 0.25 1.670 1.685 0.035 321 78.5729 0.0759 0.0004 0.0008 0.0022 0.0045 0.0087 0.0134 0.0154 0.0187 0.0239 0.0267 0.0297 0.0350 0.0342 0.0338 0.0318 0.0319 0.0309 0.0309 0.0285 0.0271 0.0281 0.0235 0.0241 0.0211 0.0213 0.0199 0.0176 0.0168 0.0141 0.0146 0.0141 0.0140 0.0124 0.0111 0.0126 0.2189 0.0116 42 0.25 1.689 1.696 0.035 319 78.9743 0.0728 0.0003 0.0011 0.0023 0.0048 0.0088 0.0126 0.0177 0.0218 0.0247 0.0283 0.0307 0.0330 0.0345 0.0364 0.0337 0.0320 0.0307 0.0307 0.0280 0.0270 0.0262 0.0248 0.0229 0.0215 0.0200 0.0193 0.0181 0.0178 0.0158 0.0147 0.0149 0.0139 0.0133 0.0118 0.0106 0.2120 0.0104 43 0.25 1.687 1.695 0.035 320 78.7174 0.0756 0.0005 0.0010 0.0029 0.0046 0.0076 0.0134 0.0172 0.0220 0.0234 0.0283 0.0311 0.0326 0.0356 0.0315 0.0335 0.0314 0.0317 0.0319 0.0295 0.0291 0.0254 0.0248 0.0221 0.0206 0.0205 0.0175 0.0183 0.0174 0.0170 0.0153 0.0144 0.0145 0.0124 0.0116 0.0102 0.2133 0.0103 44 0.25 1.665 1.682 0.035 320 78.3369 0.0800 0.0006 0.0011 0.0023 0.0051 0.0073 0.0126 0.0161 0.0196 0.0235 0.0282 0.0301 0.0344 0.0332 0.0333 0.0333 0.0306 0.0302 0.0308 0.0274 0.0283 0.0253 0.0245 0.0224 0.0211 0.0210 0.0193 0.0190 0.0169 0.0150 0.0156 0.0133 0.0133 0.0123 0.0117 0.0098 0.2195 0.0121 Model Fits - Parameters Constrained Using SAFTE Iteration EGS G UT GDEC Median_RT StdDev F_Starts LT160 LT170 LT180 LT190 LT200 LT210 LT220 LT230 LT240 LT250 LT260 LT270 LT280 LT290 LT300 LT310 LT320 LT330 LT340 LT350 LT360 LT370 LT380 LT390 LT400 LT410 LT420 LT430 LT440 LT450 LT460 LT470 LT480 LT490 LT500 Lapses S_Attacks 1 0.25 2.115 1.937 0.035 288 69.7346 0.0218 0.0001 0.0015 0.0048 0.0086 0.0157 0.0233 0.0341 0.0423 0.0501 0.0566 0.0590 0.0603 0.0597 0.0532 0.0528 0.0484 0.0422 0.0372 0.0344 0.0315 0.0277 0.0254 0.0202 0.0191 0.0168 0.0170 0.0142 0.0136 0.0103 0.0095 0.0084 0.0071 0.0066 0.0063 0.0063 0.0534 0.0005 2 0.25 2.119 1.940 0.035 288 69.7526 0.0209 0.0000 0.0014 0.0037 0.0092 0.0159 0.0262 0.0332 0.0421 0.0508 0.0524 0.0609 0.0601 0.0591 0.0556 0.0541 0.0456 0.0436 0.0378 0.0355 0.0313 0.0284 0.0238 0.0204 0.0201 0.0174 0.0149 0.0131 0.0124 0.0117 0.0117 0.0078 0.0084 0.0069 0.0059 0.0059 0.0512 0.0005 3 0.25 2.108 1.933 0.035 288 69.8096 0.0219 0.0000 0.0014 0.0039 0.0088 0.0154 0.0268 0.0334 0.0432 0.0489 0.0546 0.0599 0.0583 0.0589 0.0548 0.0518 0.0464 0.0439 0.0370 0.0347 0.0311 0.0268 0.0253 0.0224 0.0196 0.0173 0.0157 0.0135 0.0128 0.0108 0.0089 0.0089 0.0081 0.0059 0.0071 0.0061 0.0556 0.0005 4 0.25 2.096 1.926 0.035 289 70.2600 0.0245 0.0003 0.0013 0.0047 0.0081 0.0158 0.0234 0.0332 0.0409 0.0487 0.0555 0.0581 0.0580 0.0580 0.0565 0.0520 0.0471 0.0392 0.0362 0.0359 0.0312 0.0274 0.0238 0.0232 0.0200 0.0188 0.0151 0.0138 0.0128 0.0114 0.0112 0.0096 0.0083 0.0071 0.0060 0.0055 0.0567 0.0006 5 0.25 2.100 1.929 0.035 289 70.0280 0.0239 0.0002 0.0010 0.0048 0.0098 0.0147 0.0232 0.0346 0.0411 0.0495 0.0530 0.0569 0.0597 0.0590 0.0560 0.0495 0.0479 0.0416 0.0395 0.0338 0.0309 0.0280 0.0251 0.0219 0.0188 0.0176 0.0166 0.0149 0.0129 0.0123 0.0088 0.0080 0.0073 0.0072 0.0072 0.0057 0.0565 0.0005 6 0.25 2.116 1.938 0.035 287 69.8066 0.0205 0.0001 0.0016 0.0049 0.0085 0.0183 0.0246 0.0337 0.0411 0.0494 0.0568 0.0577 0.0628 0.0599 0.0568 0.0512 0.0468 0.0428 0.0376 0.0331 0.0315 0.0263 0.0258 0.0217 0.0179 0.0166 0.0152 0.0141 0.0117 0.0109 0.0109 0.0094 0.0073 0.0065 0.0064 0.0054 0.0531 0.0006 7 0.25 2.120 1.940 0.035 287 69.3153 0.0211 0.0000 0.0015 0.0045 0.0097 0.0164 0.0246 0.0339 0.0428 0.0485 0.0538 0.0596 0.0597 0.0610 0.0588 0.0528 0.0477 0.0426 0.0360 0.0359 0.0313 0.0276 0.0251 0.0211 0.0193 0.0159 0.0161 0.0132 0.0119 0.0113 0.0098 0.0091 0.0079 0.0070 0.0057 0.0050 0.0515 0.0005 8 0.25 2.087 1.922 0.035 290 70.5281 0.0244 0.0000 0.0013 0.0040 0.0087 0.0154 0.0246 0.0343 0.0431 0.0481 0.0552 0.0554 0.0562 0.0577 0.0547 0.0516 0.0470 0.0419 0.0374 0.0330 0.0312 0.0279 0.0244 0.0237 0.0209 0.0177 0.0168 0.0152 0.0132 0.0117 0.0106 0.0095 0.0085 0.0065 0.0066 0.0052 0.0555 0.0009 9 0.25 2.021 1.884 0.035 294 72.2638 0.0297 0.0000 0.0008 0.0040 0.0084 0.0136 0.0218 0.0291 0.0398 0.0444 0.0508 0.0521 0.0541 0.0520 0.0514 0.0510 0.0432 0.0406 0.0398 0.0375 0.0305 0.0286 0.0242 0.0241 0.0206 0.0187 0.0167 0.0156 0.0153 0.0118 0.0114 0.0105 0.0095 0.0078 0.0073 0.0080 0.0747 0.0007 10 0.25 1.950 1.844 0.035 299 74.0954 0.0372 0.0003 0.0010 0.0039 0.0079 0.0125 0.0204 0.0248 0.0336 0.0409 0.0447 0.0482 0.0497 0.0487 0.0496 0.0469 0.0426 0.0395 0.0368 0.0360 0.0320 0.0282 0.0230 0.0237 0.0218 0.0196 0.0178 0.0163 0.0142 0.0160 0.0121 0.0126 0.0095 0.0090 0.0086 0.0086 0.0995 0.0023 11 0.25 1.910 1.821 0.035 302 74.7477 0.0417 0.0001 0.0013 0.0042 0.0075 0.0111 0.0185 0.0240 0.0293 0.0382 0.0421 0.0463 0.0477 0.0477 0.0465 0.0438 0.0397 0.0382 0.0353 0.0322 0.0296 0.0297 0.0291 0.0253 0.0211 0.0209 0.0185 0.0171 0.0154 0.0151 0.0120 0.0117 0.0121 0.0096 0.0092 0.0075 0.1178 0.0027 12 0.25 1.909 1.820 0.035 303 74.2088 0.0406 0.0003 0.0007 0.0034 0.0063 0.0123 0.0187 0.0262 0.0313 0.0365 0.0434 0.0454 0.0470 0.0453 0.0441 0.0485 0.0442 0.0373 0.0390 0.0341 0.0319 0.0302 0.0259 0.0245 0.0235 0.0210 0.0161 0.0159 0.0159 0.0151 0.0139 0.0122 0.0101 0.0099 0.0087 0.0083 0.1089 0.0035 13 0.25 1.927 1.831 0.035 302 74.6560 0.0387 0.0001 0.0011 0.0031 0.0066 0.0133 0.0179 0.0264 0.0318 0.0375 0.0420 0.0466 0.0500 0.0502 0.0455 0.0457 0.0414 0.0417 0.0365 0.0323 0.0297 0.0299 0.0273 0.0228 0.0207 0.0206 0.0196 0.0175 0.0173 0.0140 0.0132 0.0121 0.0108 0.0110 0.0092 0.0082 0.1050 0.0028 21 0.25 1.839 1.781 0.035 306 76.0416 0.0479 0.0002 0.0008 0.0030 0.0071 0.0104 0.0176 0.0239 0.0283 0.0326 0.0374 0.0437 0.0433 0.0430 0.0444 0.0416 0.0388 0.0371 0.0326 0.0323 0.0296 0.0275 0.0262 0.0225 0.0217 0.0206 0.0187 0.0168 0.0169 0.0156 0.0131 0.0126 0.0115 0.0118 0.0100 0.0082 0.1448 0.0057 14 0.25 1.935 1.835 0.035 301 73.4974 0.0359 0.0001 0.0012 0.0031 0.0076 0.0119 0.0188 0.0248 0.0323 0.0396 0.0428 0.0457 0.0515 0.0474 0.0497 0.0440 0.0418 0.0402 0.0380 0.0356 0.0341 0.0280 0.0274 0.0255 0.0223 0.0205 0.0168 0.0173 0.0145 0.0139 0.0119 0.0119 0.0119 0.0094 0.0081 0.0068 0.1050 0.0022 15 0.25 1.925 1.830 0.035 302 74.7955 0.0374 0.0003 0.0012 0.0035 0.0066 0.0121 0.0198 0.0248 0.0310 0.0390 0.0426 0.0484 0.0497 0.0498 0.0464 0.0435 0.0426 0.0418 0.0374 0.0291 0.0330 0.0284 0.0267 0.0240 0.0212 0.0210 0.0186 0.0176 0.0155 0.0143 0.0127 0.0126 0.0098 0.0112 0.0092 0.0093 0.1055 0.0026 16 0.25 1.914 1.823 0.035 301 74.6559 0.0385 0.0002 0.0011 0.0030 0.0077 0.0117 0.0183 0.0232 0.0324 0.0380 0.0411 0.0480 0.0480 0.0509 0.0503 0.0444 0.0411 0.0382 0.0359 0.0341 0.0318 0.0259 0.0254 0.0243 0.0218 0.0206 0.0171 0.0168 0.0162 0.0147 0.0126 0.0118 0.0115 0.0113 0.0090 0.0081 0.1120 0.0029 17 0.25 1.920 1.827 0.035 301 74.4512 0.0404 0.0001 0.0012 0.0033 0.0068 0.0119 0.0173 0.0251 0.0313 0.0381 0.0442 0.0457 0.0490 0.0506 0.0465 0.0460 0.0406 0.0388 0.0355 0.0340 0.0312 0.0284 0.0256 0.0235 0.0209 0.0219 0.0185 0.0159 0.0147 0.0140 0.0128 0.0120 0.0112 0.0094 0.0094 0.0092 0.1124 0.0026 18 0.25 1.941 1.839 0.035 300 74.7565 0.0379 0.0001 0.0012 0.0038 0.0078 0.0129 0.0204 0.0280 0.0322 0.0393 0.0462 0.0465 0.0478 0.0500 0.0464 0.0441 0.0424 0.0408 0.0365 0.0321 0.0308 0.0284 0.0256 0.0244 0.0219 0.0189 0.0179 0.0179 0.0155 0.0160 0.0131 0.0120 0.0098 0.0093 0.0092 0.0081 0.1023 0.0025 19 0.25 1.949 1.843 0.035 300 74.4482 0.0357 0.0001 0.0009 0.0036 0.0078 0.0123 0.0200 0.0255 0.0323 0.0407 0.0463 0.0496 0.0487 0.0494 0.0501 0.0454 0.0418 0.0398 0.0336 0.0330 0.0323 0.0291 0.0265 0.0249 0.0229 0.0203 0.0180 0.0168 0.0154 0.0150 0.0124 0.0119 0.0113 0.0103 0.0085 0.0085 0.0978 0.0016 20 0.25 1.914 1.823 0.035 302.5 74.7469 0.0376 0.0002 0.0012 0.0031 0.0067 0.0126 0.0177 0.0253 0.0305 0.0374 0.0417 0.0450 0.0500 0.0491 0.0466 0.0466 0.0422 0.0388 0.0375 0.0333 0.0324 0.0284 0.0272 0.0228 0.0230 0.0211 0.0182 0.0182 0.0158 0.0140 0.0142 0.0133 0.0114 0.0099 0.0092 0.0085 0.1063 0.0032 22 0.25 1.759 1.736 0.035 315 77.8416 0.0602 0.0003 0.0007 0.0032 0.0053 0.0102 0.0129 0.0189 0.0236 0.0300 0.0324 0.0369 0.0371 0.0374 0.0361 0.0383 0.0356 0.0367 0.0328 0.0304 0.0280 0.0261 0.0240 0.0244 0.0245 0.0209 0.0210 0.0186 0.0177 0.0158 0.0152 0.0136 0.0129 0.0130 0.0116 0.0103 0.1766 0.0072 23 0.25 1.715 1.710 0.035 317 77.9913 0.0677 0.0004 0.0013 0.0020 0.0062 0.0083 0.0123 0.0184 0.0226 0.0264 0.0307 0.0339 0.0348 0.0364 0.0366 0.0349 0.0335 0.0326 0.0312 0.0297 0.0291 0.0258 0.0260 0.0230 0.0217 0.0210 0.0174 0.0198 0.0155 0.0151 0.0151 0.0146 0.0130 0.0121 0.0102 0.0110 0.2001 0.0098 24 0.25 1.716 1.711 0.035 317 78.4657 0.0692 0.0005 0.0014 0.0032 0.0054 0.0082 0.0123 0.0186 0.0228 0.0275 0.0310 0.0326 0.0364 0.0353 0.0352 0.0348 0.0340 0.0311 0.0293 0.0288 0.0284 0.0273 0.0257 0.0215 0.0224 0.0208 0.0206 0.0191 0.0185 0.0174 0.0159 0.0134 0.0119 0.0135 0.0103 0.0102 0.1948 0.0109 25 0.25 1.740 1.724 0.035 316 77.6092 0.0614 0.0004 0.0011 0.0027 0.0055 0.0088 0.0136 0.0192 0.0227 0.0291 0.0324 0.0342 0.0357 0.0364 0.0373 0.0370 0.0355 0.0332 0.0336 0.0312 0.0302 0.0263 0.0268 0.0228 0.0229 0.0207 0.0215 0.0179 0.0169 0.0163 0.0141 0.0138 0.0135 0.0121 0.0110 0.0101 0.1835 0.0088 26 0.25 1.752 1.731 0.035 313 77.8173 0.0623 0.0003 0.0013 0.0028 0.0049 0.0104 0.0148 0.0173 0.0236 0.0291 0.0343 0.0350 0.0372 0.0390 0.0397 0.0365 0.0357 0.0339 0.0316 0.0305 0.0287 0.0273 0.0254 0.0227 0.0196 0.0194 0.0199 0.0181 0.0173 0.0176 0.0144 0.0142 0.0127 0.0125 0.0109 0.0099 0.1811 0.0080 27 0.25 1.742 1.726 0.035 314 77.8212 0.0636 0.0002 0.0008 0.0024 0.0058 0.0084 0.0151 0.0190 0.0243 0.0287 0.0319 0.0345 0.0360 0.0407 0.0360 0.0375 0.0323 0.0358 0.0326 0.0303 0.0292 0.0261 0.0262 0.0225 0.0227 0.0199 0.0181 0.0189 0.0161 0.0161 0.0145 0.0132 0.0136 0.0120 0.0117 0.0103 0.1832 0.0094 28 0.25 1.731 1.720 0.035 316 77.6943 0.0665 0.0003 0.0009 0.0028 0.0060 0.0090 0.0121 0.0178 0.0215 0.0272 0.0305 0.0345 0.0356 0.0394 0.0383 0.0351 0.0363 0.0311 0.0318 0.0302 0.0300 0.0273 0.0240 0.0247 0.0228 0.0210 0.0196 0.0184 0.0156 0.0153 0.0155 0.0145 0.0140 0.0130 0.0114 0.0095 0.1882 0.0083 29 0.25 1.740 1.725 0.035 315 77.4377 0.0651 0.0003 0.0009 0.0033 0.0044 0.0090 0.0144 0.0188 0.0255 0.0260 0.0310 0.0359 0.0359 0.0380 0.0358 0.0354 0.0356 0.0349 0.0319 0.0313 0.0297 0.0272 0.0244 0.0235 0.0227 0.0215 0.0185 0.0176 0.0166 0.0155 0.0166 0.0146 0.0126 0.0117 0.0108 0.0091 0.1862 0.0080 30 0.25 1.767 1.740 0.035 313 77.1510 0.0590 0.0006 0.0008 0.0030 0.0049 0.0091 0.0145 0.0202 0.0260 0.0283 0.0342 0.0364 0.0390 0.0389 0.0406 0.0351 0.0359 0.0348 0.0294 0.0337 0.0284 0.0286 0.0277 0.0240 0.0213 0.0192 0.0185 0.0198 0.0158 0.0161 0.0157 0.0128 0.0129 0.0116 0.0096 0.0102 0.1749 0.0084 31 0.25 1.778 1.746 0.035 312 76.8173 0.0579 0.0002 0.0014 0.0024 0.0050 0.0092 0.0144 0.0201 0.0266 0.0306 0.0352 0.0381 0.0385 0.0434 0.0396 0.0385 0.0364 0.0365 0.0347 0.0309 0.0278 0.0260 0.0265 0.0244 0.0228 0.0213 0.0202 0.0178 0.0165 0.0158 0.0146 0.0137 0.0121 0.0114 0.0110 0.0103 0.1607 0.0074 32 0.25 1.740 1.725 0.035 315 78.2394 0.0680 0.0001 0.0013 0.0029 0.0058 0.0093 0.0138 0.0183 0.0225 0.0284 0.0285 0.0357 0.0398 0.0377 0.0383 0.0355 0.0337 0.0330 0.0314 0.0309 0.0286 0.0266 0.0247 0.0241 0.0222 0.0210 0.0210 0.0169 0.0159 0.0158 0.0149 0.0148 0.0127 0.0128 0.0111 0.0114 0.1825 0.0079 33 0.25 1.658 1.678 0.035 320 78.5700 0.0800 0.0004 0.0014 0.0024 0.0045 0.0074 0.0122 0.0150 0.0198 0.0238 0.0283 0.0299 0.0335 0.0325 0.0334 0.0324 0.0322 0.0300 0.0305 0.0272 0.0278 0.0241 0.0232 0.0235 0.0218 0.0185 0.0187 0.0159 0.0174 0.0179 0.0141 0.0130 0.0135 0.0115 0.0112 0.0115 0.2278 0.0117 34 0.25 1.569 1.628 0.035 326 79.8051 0.0998 0.0005 0.0009 0.0031 0.0040 0.0068 0.0099 0.0125 0.0172 0.0222 0.0219 0.0246 0.0275 0.0267 0.0283 0.0269 0.0302 0.0263 0.0272 0.0257 0.0258 0.0222 0.0213 0.0201 0.0200 0.0187 0.0174 0.0177 0.0167 0.0160 0.0141 0.0133 0.0131 0.0109 0.0109 0.0127 0.2688 0.0180 35 0.25 1.520 1.600 0.035 330 80.8923 0.1114 0.0004 0.0014 0.0021 0.0042 0.0063 0.0087 0.0116 0.0147 0.0168 0.0214 0.0239 0.0249 0.0272 0.0262 0.0226 0.0259 0.0249 0.0245 0.0222 0.0218 0.0223 0.0204 0.0205 0.0202 0.0171 0.0185 0.0159 0.0152 0.0146 0.0156 0.0138 0.0135 0.0132 0.0122 0.0112 0.2939 0.0188 36 0.25 1.523 1.602 0.035 333 80.7717 0.1097 0.0004 0.0011 0.0021 0.0035 0.0061 0.0097 0.0124 0.0135 0.0174 0.0211 0.0218 0.0213 0.0257 0.0251 0.0245 0.0239 0.0237 0.0234 0.0241 0.0228 0.0239 0.0217 0.0186 0.0191 0.0188 0.0176 0.0176 0.0163 0.0153 0.0138 0.0135 0.0146 0.0132 0.0120 0.0116 0.3008 0.0182 37 0.25 1.552 1.618 0.035 328 80.5716 0.1026 0.0006 0.0007 0.0029 0.0038 0.0064 0.0098 0.0128 0.0156 0.0197 0.0222 0.0262 0.0272 0.0279 0.0248 0.0265 0.0268 0.0249 0.0258 0.0247 0.0215 0.0208 0.0220 0.0201 0.0195 0.0201 0.0172 0.0189 0.0178 0.0150 0.0157 0.0139 0.0136 0.0126 0.0119 0.0109 0.2802 0.0166 38 0.25 1.568 1.627 0.035 328 80.2848 0.1002 0.0004 0.0012 0.0024 0.0035 0.0064 0.0108 0.0135 0.0178 0.0178 0.0219 0.0238 0.0270 0.0266 0.0272 0.0274 0.0264 0.0261 0.0267 0.0257 0.0224 0.0222 0.0225 0.0212 0.0197 0.0190 0.0186 0.0182 0.0172 0.0147 0.0152 0.0144 0.0124 0.0133 0.0119 0.0118 0.2756 0.0166 39 0.25 1.559 1.622 0.035 330 79.6560 0.1035 0.0006 0.0009 0.0019 0.0037 0.0057 0.0087 0.0127 0.0163 0.0194 0.0232 0.0244 0.0274 0.0270 0.0273 0.0271 0.0273 0.0255 0.0255 0.0262 0.0245 0.0230 0.0233 0.0212 0.0198 0.0187 0.0183 0.0189 0.0170 0.0162 0.0146 0.0144 0.0137 0.0121 0.0109 0.0119 0.2682 0.0188 40 0.25 1.548 1.616 0.035 328 79.4461 0.1042 0.0004 0.0010 0.0019 0.0031 0.0062 0.0087 0.0127 0.0164 0.0189 0.0213 0.0262 0.0267 0.0278 0.0246 0.0284 0.0276 0.0277 0.0239 0.0242 0.0248 0.0224 0.0221 0.0224 0.0199 0.0179 0.0166 0.0161 0.0150 0.0136 0.0145 0.0138 0.0141 0.0117 0.0131 0.0108 0.2814 0.0181 41 0.25 1.561 1.623 0.035 326 80.2042 0.1030 0.0007 0.0011 0.0021 0.0049 0.0068 0.0099 0.0140 0.0169 0.0202 0.0246 0.0246 0.0262 0.0277 0.0275 0.0268 0.0266 0.0266 0.0261 0.0260 0.0234 0.0213 0.0230 0.0219 0.0185 0.0180 0.0188 0.0166 0.0157 0.0169 0.0143 0.0139 0.0128 0.0114 0.0120 0.0108 0.2706 0.0178 42 0.25 1.592 1.641 0.035 326 79.8899 0.0897 0.0008 0.0012 0.0023 0.0046 0.0076 0.0109 0.0140 0.0168 0.0197 0.0229 0.0261 0.0301 0.0300 0.0301 0.0274 0.0278 0.0274 0.0265 0.0264 0.0236 0.0249 0.0238 0.0239 0.0212 0.0204 0.0180 0.0157 0.0169 0.0151 0.0139 0.0138 0.0137 0.0134 0.0108 0.0125 0.2593 0.0166 43 0.25 1.606 1.649 0.035 326 80.1060 0.0914 0.0003 0.0012 0.0019 0.0042 0.0078 0.0116 0.0134 0.0184 0.0202 0.0272 0.0277 0.0306 0.0274 0.0294 0.0273 0.0281 0.0291 0.0265 0.0275 0.0256 0.0229 0.0223 0.0222 0.0213 0.0191 0.0186 0.0170 0.0183 0.0171 0.0158 0.0151 0.0133 0.0127 0.0132 0.0112 0.2473 0.0157 44 0.25 1.567 1.627 0.035 329 80.5002 0.1030 0.0005 0.0006 0.0023 0.0044 0.0058 0.0105 0.0129 0.0167 0.0195 0.0219 0.0254 0.0275 0.0264 0.0291 0.0280 0.0267 0.0236 0.0247 0.0235 0.0243 0.0215 0.0223 0.0210 0.0203 0.0193 0.0171 0.0191 0.0178 0.0146 0.0170 0.0142 0.0131 0.0129 0.0127 0.0112 0.2730 0.0156 |# ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;; ;;; Experiment Delivery Code ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; (defun pvt (duration &key (sleep-threshold 30000) (visible t) (gdecrement 0)) (let ((data nil) (time-list nil) (start-block (if *actr-enabled-p* (pm-get-time) (get-internal-real-time))) (delay 0) (response nil) (start-time nil) (gval (no-output (car (sgp :g)))) (box (make-box 60 100 140 140))) ;;Set up experiment window (setf *w* (open-exp-window "Task Window" :x 200 :y 200 :width 200 :height 200 :visible visible)) (pm-install-device *w*) (if *actr-enabled-p* (pm-run 1 :full-time t) (sleep 1)) ;;Main trial loop - run for the duration of the session (until (> (current-interval start-block) duration) (clear-exp-window) (pm-proc-display) ;; Make sure there are no carryover effects from previous trial (no-output (eval (list 'mod-chunk (get-safe-wme (goal-focus)) 'state 'wait))) (no-output (eval (list 'sgp ':g gval))) (clear (vis-m *mp*)) ;; Prepare for trial (setf delay (+ 2 (random 8))) (setf *response* nil) ;; (format t "~%******************* DELAY STARTING *******************") (setf start-time (if *actr-enabled-p* (pm-get-time) (get-internal-real-time))) ;; (format t "~%Delay - ~a" (pm-get-time)) (if *actr-enabled-p* (loop (pm-run (- delay (- (mp-time *mp*) start-time))) (cond ((>= (- (* (mp-time *mp*) 1000) start-time) (* delay 1000)) (return t)) (*response* (return t))) ; Tell the system that cognition is busy (setf (ps-state (cog-m *mp*)) :pending-execution) ; explicitly schedule the next conflict-resolution for the later time that you want (queue-command :time (rand-time (first (no-output (sgp :dat)))) :where :COGNITION :command 'select-production :params (list 1) :priority -1)) (loop (cond ((or *response* (> (- (get-internal-real-time) start-time) (* delay 1000))) (return t))) (sleep .01))) (when *response* (setf *response* start-time)) (when (null *response*) (add-text-to-exp-window :text (format nil "START") :color 'red :x 100 :y 100 :width 80 :height 80) (setf *gdecrement* gdecrement) (setf start-time (if *actr-enabled-p* (pm-get-time) (get-internal-real-time))) (pm-proc-display) ;; (format t "~%******************* DISPLAY STIM *******************") (if *actr-enabled-p* (loop (pm-run (- sleep-threshold (- (mp-time *mp*) start-time))) (cond ((>= (- (* (mp-time *mp*) 1000) start-time) (* sleep-threshold)) (return t)) (*response* (return t))) ; Decrement G (no-output (eval (list 'sgp ':g (max .000001 (- (car (sgp :g)) *gdecrement*))))) ; Tell the system that cognition is busy (setf (ps-state (cog-m *mp*)) :pending-execution) ; explicitly schedule the next conflict-resolution for the later time that you want (queue-command :time (rand-time (first (no-output (sgp :dat)))) :where :COGNITION :command 'select-production :params (list 1) :priority -1)) (loop (cond (*response* (return t))) (process-single-event))) (clear-exp-window) (pm-proc-display)) (if (and (null *response*) *actr-enabled-p*) (setf *response* (pm-get-time))) (setf data (cons (list delay (round (- *response* start-time))) data)) (add-text-to-exp-window :text (format nil "~a" (round (- *response* start-time))) :color 'blue :x 80 :y 90 :width 40 :height 20) (if *actr-enabled-p* (pm-run .5 :full-time t) (sleep .5))) (clear-exp-window) (add-text-to-exp-window :text "Done" :color 'red :x 80 :y 90 :width 40 :height 20) (reverse data) )) (defun current-interval (start) (let ((cur-time (if *actr-enabled-p* (pm-get-time) (get-internal-real-time)))) (/ (- cur-time start) 1000) )) ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;; ;;; Batch runs ;;; -Multiple Runs ;;; -Multiple Parameter Combinations ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; (defvar *pm-parameter-list* nil) (defun pvt-batch (n &key (summary-file "batch-data.txt") (summary-2 "statistics.txt") (visible nil) (params nil) (interval 10) (ind-data nil) (session-length 600) (lapse-threshold 500) (lapse-interval 500) (sleep-threshold 30000) (false-start-threshold 150) (human-data *baseline*) (gdecrement 0)) (setf *all* (make-hash-table :test #'equalp)) (setf *gdecrement* gdecrement) (setf (gethash (list 'fstarts) *all*) 0) (setf (gethash (list 'lapses) *all*) 0) (setf (gethash (list 'sleeps) *all*) 0) (setf (gethash (list 'total) *all*) 0) (let ((start-batch (get-internal-real-time)) (params-list nil) (iterations nil) (cur-run-data nil)) (setf summary-file (format nil "~a~a" *all-path* summary-file)) (setf summary-2 (format nil "~a~a" *all-path* summary-2)) (loop (cond ((null params) (return t))) (setf params-list (cons (make-list-of-values (car params)) params-list)) (setf params (cdr params))) (setf iterations (make-all-instances params-list)) (if (null iterations) (setf iterations '(()))) (with-open-file (summary summary-file :direction :output :if-exists :append :if-does-not-exist :create) (with-open-file (stats summary-2 :direction :output :if-exists :append :if-does-not-exist :create) (format summary "Iteration~c" #\tab) (format stats "Iteration~c" #\tab) (mapcar #'(lambda (parameter) (format summary "~a~c" (first parameter) #\tab) (format stats "~a~c" (first parameter) #\tab) ) (first iterations)) (format summary "Count~cALL_Mean_RT~cNRML-Mean_RT~cMedian_RT~cStdDev~cF_Starts~c" #\tab #\tab #\tab #\tab #\tab #\tab) (format stats "Count~cALL_Mean_RT~cNRML-Mean_RT~cMedian_RT~cStdDev~cF_Starts~c" #\tab #\tab #\tab #\tab #\tab #\tab) (dotimes (x (- (round (/ lapse-threshold interval)) (round (/ 150 interval)))) (format summary "LT~a~c" (* interval (+ (1+ x) (/ 150 interval))) #\tab) (format stats "LT~a~c" (* interval (+ (1+ x) (/ 150 interval))) #\tab)) (format summary "Lapses~cS_Attacks" #\tab) (format stats "Lapses~cS_Attacks" #\tab) (format summary "~%") (format stats "~%") (dotimes (x (length iterations)) (setf *actr-parameter-list* nil) (clrhash *all*) (setf (gethash (list 'fstarts) *all*) 0) (setf (gethash (list 'lapses) *all*) 0) (setf (gethash (list 'sleeps) *all*) 0) (setf (gethash (list 'total) *all*) 0) (dotimes (y n) (setf *rpm-parameter-list* nil) (mapcar #'(lambda (parameter) (cond ((equal (third parameter) 'rpm) (setf *rpm-parameter-list* (append (list (first parameter) (second parameter)) *rpm-parameter-list*))) ((equal (third parameter) 'actr) (setf *actr-parameter-list* (append (list (first parameter) (second parameter)) *actr-parameter-list*))) (t (eval (list 'setf (first parameter) (second parameter))))) ) (nth x iterations)) (reset) (setf cur-run-data (pvt session-length :sleep-threshold sleep-threshold :visible visible :gdecrement *gdecrement*)) (if ind-data (make-single-subject-file cur-run-data (nth x iterations) (1+ y))) (add-data-to-summary-file cur-run-data (nth x iterations) (1+ y) summary interval lapse-threshold sleep-threshold false-start-threshold) ) (format t "~a~c~c" (1+ x) #\. #\tab) (format stats "~a~c" (1+ x) #\tab) (do-stats interval lapse-threshold sleep-threshold false-start-threshold human-data (nth x iterations) stats) ;; (print-data-from-hash interval lapse-threshold sleep-threshold false-start-threshold lapse-interval) ) ; (print-data-from-hash interval lapse-threshold sleep-threshold false-start-threshold lapse-interval) ; (do-stats interval lapse-threshold sleep-threshold false-start-threshold human-data) (Format t "Run Time:~c~a~%~%" #\tab (- (get-internal-real-time) start-batch)) )) )) (defun make-list-of-values (conditions) ; (pprint conditions) (let ((keyword (parameter-keyword conditions)) (which-system (parameter-system conditions)) (value-list nil) (value (parameter-min conditions))) (loop (cond ((> value (parameter-max conditions)) (return t))) (setf value-list (cons value value-list)) (incf value (parameter-interval conditions)) (setf value (/ (round value .00001) 100000.0))) (setf value-list (cons which-system (reverse value-list))) (cons keyword value-list) )) (defun make-all-instances (parameters) (let ((keyword nil) (which-system nil) (full-set nil)) (dotimes (x (length parameters)) ; (pprint (nth x parameters)) (setf keyword (first (nth x parameters))) (setf which-system (second (nth x parameters))) (let ((copy-set (copy-list full-set)) (temp nil)) (setf full-set nil) (cond ((null copy-set) (mapcar #'(lambda (value) (if (numberp value) (setf full-set (cons (list (list keyword value which-system)) full-set)))) (nth x parameters))) (t (mapcar #'(lambda (value) (mapcar #'(lambda (instance) (setf temp (cons (list keyword value which-system) instance)) (if (numberp value) (setf full-set (cons temp full-set)))) copy-set)) (nth x parameters))) ) )) full-set )) ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;; ;;; Data Management ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; (defun make-single-subject-file (data settings iteration) (let ((filename (format nil "~ap~a" *ind-path* iteration)) (cur-param nil) (trial-num 0)) (dotimes (x (length settings)) (setf cur-param (format nil "~a" (first (nth x settings)))) ;(setf (schar cur-param 0) #\-) (setf filename (format nil "~a-~a~a" filename cur-param (second (nth x settings)))) ) (setf filename (format nil "~a~ctxt" filename #\.)) (with-open-file (dataout filename :direction :output :if-exists :append :if-does-not-exist :create) (format dataout "Trial~cDelay~cRT~%" #\tab #\tab) (mapcar #'(lambda (trial) (format dataout "~a~c~a~c~a~%" (incf trial-num) #\tab (first trial) #\tab (second trial)) ) data) ) )) ;;Write our summary data for single model run (1 session) (defun add-data-to-summary-file (data settings iteration summary interval lapse-threshold sleep-threshold false-start-threshold) (let ((rt 0) (cur-count 0) (responses 0) (false-starts 0) (lapses 0) (sleeps 0) (last-count 0) (criterion 0) (nrml-rspns nil) (stdev 0)) (format summary "~a~c" iteration #\tab) (dotimes (x (length settings)) (setf cur-param (format nil "~a" (first (nth x settings)))) (format summary "~a~c" (second (nth x settings)) #\tab) ) (mapcar #'(lambda (trial) (if (gethash (list (second trial)) *all*) (setf (gethash (list (second trial)) *all*) (1+ (gethash (list (second trial)) *all*))) (setf (gethash (list (second trial)) *all*) 1)) (incf responses) (incf rt (second trial)) (cond ((>= (second trial) sleep-threshold) (incf sleeps)) ((>= (second trial) lapse-threshold) (incf lapses)) ((< (second trial) false-start-threshold) (incf false-starts)) (t (setf nrml-rspns (cons (second trial) nrml-rspns)))) ) data) (setf stdev (stdev nrml-rspns)) (format summary "~a~c~a~c~a~c~a~c~a~c~a~c" responses #\tab (/ rt responses 1.0) #\tab (get-average nrml-rspns) #\tab (median nrml-rspns) #\tab stdev #\tab (/ false-starts responses 1.0) #\tab) (setf (gethash (list 'fstarts) *all*) (+ (gethash (list 'fstarts) *all*) false-starts)) (setf (gethash (list 'lapses) *all*) (+ (gethash (list 'lapses) *all*) lapses)) (setf (gethash (list 'sleeps) *all*) (+ (gethash (list 'sleeps) *all*) sleeps)) (setf (gethash (list 'total) *all*) (+ (gethash (list 'total) *all*) responses)) (setf last-count false-starts) (setf criterion (+ false-start-threshold interval)) (loop (mapcar #'(lambda (trial) (if (< (second trial) criterion) (incf cur-count)) ) data) (format summary "~a~c" (/ (- cur-count last-count) responses 1.0) #\tab) (setf last-count cur-count) (setf cur-count 0) (cond ((>= criterion lapse-threshold) (return t))) (setf criterion (+ criterion interval))) (format summary "~a~c~a~c~%" (/ lapses responses 1.0) #\tab (/ sleeps responses 1.0) #\tab) )) ;;Format summary data to standard out (defun print-data-from-hash (interval lapse-threshold sleep-threshold false-start-threshold lapse-interval) (let ((value (+ false-start-threshold interval)) (current 0)) (format t "False Starts:~c~a~%" #\tab (if (numberp (gethash (list 'fstarts) *all*)) (/ (gethash (list 'fstarts) *all*) (gethash (list 'total) *all*) 1.0) 0)) (loop (format t "~c~a~c~c~a~%" #\< value #\: #\tab (if (numberp (gethash (list value) *all*)) (/ (gethash (list value) *all*) (gethash (list 'total) *all*) 1.0) 0)) (setf value (+ value interval)) (cond ((> value lapse-threshold) (return t)))) (format t "Lapses:~c~a~%" #\tab (if (numberp (gethash (list 'lapses) *all*)) (/ (gethash (list 'lapses) *all*) (gethash (list 'total) *all*) 1.0) 0)) (format t "Sleep Attacks:~c~a~%" #\tab (if (numberp (gethash (list 'sleeps) *all*)) (/ (gethash (list 'sleeps) *all*) (gethash (list 'total) *all*) 1.0) 0)) (format t "~%~%") )) ;;Write out summary data for multiple model runs with matching parameter values (N runs in pvt-batch) (defun do-stats (interval lapse-threshold sleep-threshold false-start-threshold human-data settings stat-file) (let ((list-of-nrml-data nil) (val 0) (value false-start-threshold) (stop-value) (cumulative-count 0) (total (gethash (list 'total) *all*)) (stdev 0)) (loop (cond ((> value lapse-threshold) (return t))) (setf val (gethash (list value) *all*)) (if (numberp val) (dotimes (x val) (setf list-of-nrml-data (cons value list-of-nrml-data)))) (incf value)) (dotimes (x (length settings)) (setf cur-param (format nil "~a" (first (nth x settings)))) (format stat-file "~a~c" (second (nth x settings)) #\tab) (format t "~a: ~a;~c" (first (nth x settings)) (second (nth x settings)) #\tab) ) (setf stdev (stdev list-of-nrml-data)) (format stat-file "~a~c~a~c~a~c~a~c~a~c~a~c" total #\tab 'NC #\tab (get-average list-of-nrml-data) #\tab (median list-of-nrml-data) #\tab stdev #\tab (if (numberp (gethash (list 'fstarts) *all*)) (/ (gethash (list 'fstarts) *all*) total 1.0) 0) #\tab) (setf value false-start-threshold) (setf cumulative-count (+ cumulative-count (if (numberp (gethash (list 'fstarts) *all*)) (gethash (list 'fstarts) *all*) 0))) (setf val 0) (loop (cond ((>= value lapse-threshold) (return t))) (setf cumulative-count (+ cumulative-count val)) (setf val 0) (setf stop-value (+ value interval)) (loop (setf val (+ val (if (numberp (gethash (list value) *all*)) (gethash (list value) *all*) 0))) (incf value) (cond ((>= value stop-value) (return t)))) (format stat-file "~a~c" (/ val total 1.0) #\tab) ) (format stat-file "~a~c~a~c" (if (numberp (gethash (list 'lapses) *all*)) (/ (gethash (list 'lapses) *all*) total 1.0) 0) #\tab (if (numberp (gethash (list 'sleeps) *all*)) (/ (gethash (list 'sleeps) *all*) total 1.0) 0) #\tab) (format t "Responses:~c~a~%~%" #\tab (gethash (list 'total) *all*)) (format stat-file "~%") )) ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;; Statistical Functions ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; (defun median (values) (if values (let* ((srtd (sort values '<)) (middle (/ (length srtd) 2))) (if (oddp (length srtd)) (nth (floor middle) srtd) (/ (+ (nth middle srtd) (nth (1- middle) srtd)) 2.0)) ) 'x)) (defun stdev (values) (if (and values (> (length values) 1)) (let ((avgrt (get-average values)) (summation 0)) (mapcar #'(lambda (value) (setf summation (+ summation (* (- value avgrt) (- value avgrt))))) values) (sqrt (/ summation (1- (length values))))) 'x)) (defun get-average (values) (if values (let ((sum 0)) (mapcar #'(lambda (value) (setf sum (+ value sum))) values) (/ sum (length values) 1.0)) 'x)) ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;; ;;; ACT-R Code ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;; Support Code ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;;Revise key event handler to capture response information (defmethod rpm-window-key-event-handler ((win rpm-window) key) (setf *response* (if *actr-enabled-p* (pm-get-time) (get-internal-real-time))) t) ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;; Initialization ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; (clear-all) (pm-reset) (setf *actr-enabled-p* t) ;;ACT-R Globals (sgp :er t :era t :esc t :v nil :pm t :egs .25); :pl t) ; ;;PM Globals (pm-set-params :real-time nil :show-focus nil) (pm-set-params :randomize-time 2) ;;Apply the randomize-time function to cognitive cycle times (setf *rand-ps* t) ;;Apply ACT-R fatigue parameters, if any ;;-These are constructed in the task management code (if *actr-parameter-list* (sgp-fct *actr-parameter-list*)) ;;Apply PM fatigue parameters, if any ;;-These are constructed in the task management code (eval (cons 'pm-set-params *rpm-parameter-list*)) ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;; Chunk-types and declarative memory ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; (chunk-type do-pvt state) (add-dm (trial isa do-pvt state wait) (wait isa chunk) (looking isa chunk) (done isa chunk)) (goal-focus trial) ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; ;; Productions ;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;; (P wait =goal> ISA do-pvt state wait =manual-state> ISA module-state modality free =visual-state> ISA module-state modality free !eval! (null *visual-location*) ==> =goal> ) (P attend =goal> ISA do-pvt state wait =visual-location> ISA visual-location =visual-state> ISA module-state modality free ==> +visual> ISA visual-object screen-pos =visual-location =goal> state looking ) (P respond =goal> ISA do-pvt state =state - state done =visual-location> isa visual-location =visual> ISA visual-object screen-pos =visual-location =manual-state> ISA module-state modality free =visual-state> ISA module-state modality free ==> +manual> ISA punch hand right finger index -visual> -visual-location> =goal> ISA do-pvt state done ) (P press-key =goal> ISA do-pvt state =state - state done =manual-state> ISA module-state modality free =visual-state> ISA module-state modality free ==> +manual> ISA punch hand right finger index -visual> -visual-location> =goal> state done ) (spp press-key :p 0) #| Sample calls for model runs... *These runs replicate the settings used to fit the human data aggregated across days -Figure 4. (pvt-batch 100 :params '((:egs actr 0.25 0.25 1) (:g actr 1.98 1.98 1) (:ut actr 1.84 1.84 1)) :interval 10 :lapse-threshold 500 :gdecrement 0.035) (pvt-batch 100 :params '((:egs actr 0.25 0.25 1) (:g actr 1.80 1.80 1) (:ut actr 1.78 1.78 1)) :interval 10 :lapse-threshold 500 :gdecrement 0.035) (pvt-batch 100 :params '((:egs actr 0.25 0.25 1) (:g actr 1.66 1.66 1) (:ut actr 1.70 1.70 1)) :gdecrement 0.035) (pvt-batch 100 :params '((:egs actr 0.25 0.25 1) (:g actr 1.58 1.58 1) (:ut actr 1.64 1.64 1)) :gdecrement 0.035) |#