dArrayNptsd = deque([0.]) plotWidth = 500 agents = [] totInput = [0.] minSep = .0025 rSigma = .3 dataFileName = '' # Initialize random seed seed(3) if doingLogging: writelog.init('run.log') # Codes for four-point Bezier spline codes = [Path.MOVETO, Path.CURVE4, Path.CURVE4, Path.CURVE4, ] #%%########[ probRaceEthncy ]####### def probDemographics(boroIndx): # # For each demographic category [Income, Age, Sex, Race-Ethnicity] # establish thresholds such that a random % variable rand()*100. is # likely to fall into the appropriate category matching demographics
lastX = 0. lastT = 0. dt = 1 dArray = deque([0.]) tArray = deque([0.]) plotWidth = 500 isFemale = 0 isYoung = 1 race = 'white' ethn = 'hispanic' income = 'lowrange' hadStressors = 0 nSandyTraumas = 0 PTSD_sx = 0 writelog.init('PTSD_Combo1.log') writelog.write('sex \t age \t race \t ethn \t\t income \t hadStr nSandTr exp(P)\n') writelog.write('=== \t === \t ==== \t ==== \t\t ====== \t ====== ======= ======\n\n') plt.close('all') fig = plt.figure(figsize=(12,6)) fig.canvas.set_window_title('PTSD_Distribution.py') # Create axes boxes in a linspace xorg line xorg = np.linspace(.05,.9,7) axSex = plt.axes([xorg[0],.6,.08,.2]); axSex.set_title('Sex') axAge = plt.axes([xorg[1],.6,.08,.2]); axAge.set_title('Age') axRace = plt.axes([xorg[2],.6,.08,.2]); axRace.set_title('Race') axEthn = plt.axes([xorg[3],.6,.08,.2]); axEthn.set_title('Ethn') axIncome = plt.axes([xorg[4],.6,.08,.2]); axIncome.set_title('Income') axStress = plt.axes([xorg[5],.6,.08,.2]); axStress.set_title('Stress')
lastT = 0. dt = 10 dArray = deque([0.]) tArray = deque([0.]) plotWidth = 500 isMale = 0 isOld = 1 isBlack = 0. race = 'white' ethn = 'hispanic' # income = 'lowrange' hadPrior = 0 logitCBT = 0 probCBT = 0 writelog.init('PTSD_Logistic.log') writelog.write('sex \t age \t race \t ethn \t\t nPrior \t logitCBT \t\t probCBT \n') writelog.write('=== \t === \t ==== \t ==== \t\t ====== \t ====== \t\t ======= \n\n') plt.close('all') fig = plt.figure(figsize=(12,6)) fig.canvas.set_window_title('PTSD_LogisticCombo.py') # Create axes boxes in a linspace xorg line xorg = np.linspace(.05,.9,5) axSex = plt.axes([xorg[0],.6,.08,.2]); axSex.set_title('Sex') axAge = plt.axes([xorg[1],.6,.08,.2]); axAge.set_title('Age') axRace = plt.axes([xorg[2],.6,.08,.2]); axRace.set_title('Race') axEthn = plt.axes([xorg[3],.6,.08,.2]); axEthn.set_title('Ethn') # axIncome = plt.axes([xorg[4],.6,.08,.2]); axIncome.set_title('Income') # axStress = plt.axes([xorg[5],.6,.08,.2]); axStress.set_title('Stress')