Practical: Phenotypic Factor Analysis

Size: px
Start display at page:

Download "Practical: Phenotypic Factor Analysis"

Transcription

1 Practical: Phnotypic Factor Analysis Big 5 dimnsions Nuroticism & Extravrsion in 361 fmal UvA studnts - Exploratory Factor Analysis (EFA) using R (factanal) with Varimax and Promax rotation - Confirmatory Factor Analysis (CFA) using OpnMx Dolan & Abdllaoui Bouldr Workshop 2016

2 Nuroticism: idntifis individuals who ar pron to psychological distrss n1 - Anxity: lvl of fr floating anxity n2 - Angry Hostility: tndncy to xprinc angr and rlatd stats such as frustration and bittrnss n3 - Dprssion: tndncy to xprinc flings of guilt, sadnss, dspondncy and lonlinss n4 - Slf-Consciousnss: shynss or social anxity n5 - Impulsivnss: tndncy to act on cravings and urgs rathr than dlaying gratification n6 - Vulnrability: gnral suscptibility to strss Extravrsion: quantity and intnsity of nrgy dirctd outwards into th social world 1 - Warmth: intrst in and frindlinss towards othrs 2 - Grgariousnss: prfrnc for th company of othrs 3 - Assrtivnss: social ascndancy and forcfulnss of xprssion 4 - Activity: pac of living 5 - Excitmnt Sking: nd for nvironmntal stimulation 6 - Positiv Emotions: tndncy to xprinc positiv motions 2

3 # Part 1: rad th data EFA # clar th mmory rm(list=ls(all=true)) # load OpnMx library(opnmx) # st workingdirctory stwd( YOUR_WORKING_DIRECTORY") # rad th data datb5=rad.tabl('rdataf') # rad th fmal data # assign variabl nams varlabs=c('sx', 'n1', 'n2', 'n3', 'n4', 'n5', 'n6', '1', '2', '3', '4', '5', '6', 'o1', 'o2', 'o3', 'o4', 'o5', 'o6', 'a1', 'a2', 'a3', 'a4', 'a5', 'a6', 'c1', 'c2', 'c3', 'c4', 'c5', 'c6') colnams(datb5)=varlabs # slct th variabls of intrst isl=c(2:13) # slction of variabls n1-n6, 1-6 datb2=datb5[,isl] # th data fram that w'll us blow. 3

4 # Part 2: summary statistics EFA Ss1=cov(datb2[,1:12]) # calculat th covarianc matrix in fmals print(round(ss1,1)) Rs1=cov2cor(Ss1) print(round(rs1,2)) # convrt to corrlation matrix Ms1=apply(datb2[,1:12],2,man) print(round(ms1,2)) # End of part 2 # fmals mans > print(round(ms1,2)) n1 n2 n3 n4 n5 n > print(round(rs1,2)) n1 n2 n3 n4 n5 n n n n n n n

5 EFA SCREEPLOT ignvalus How many factors? Ambiguous... - Eignvalus > 1 suggsts 3 factors - Elbow critrion suggsts 2 factors :12 5

6 EFA S y = Y t + Q 1 1 N E n1 n2 n3 n4 n5 n

7 S y = Y t + Q EFA Matrix (12x2) of factor loadings (not: loading <.1 not shown): Factor1 Factor2 n n n n n n n1 n2 n3 1 1 N n4 n5 n E Factor Corrlation matrix (2x2) Y: Factor1 Factor2 Factor Factor Diagonal covarianc matrix (12x12) of rsiduals (Q): n1 n2 n3 n4 n5 n

8 unrotatd rotatd 8

9 EFA Goodnss of fit of EFA 2 factor modl. Tst of th hypothsis that 2 factors ar sufficint. Th chi squar statistic is on 43 dgrs of frdom. Th p-valu is By this statistical critrion th modl is judgd to b accptabl if th p-valu is gratr than th chosn alpha (.g. alpha=.05). By th statistical critrion, w d rjct th modl! 9

10 # Part 4A: saturatd modl CFA ny=12 # numbr of indicators n=2 # xpctd numbr of common factors varnams=colnams(datb2) # var nams ### fit th saturatd modl ### # dfin th mans and covarianc matrix in OpnMx to obtain th saturatd modl logliklihood Rs1=mxMatrix(typ='Stand',nrow=ny,ncol=ny,fr=TRUE,valu=.05, lbound=-.9,ubound=.9,nam='cor') Sds1=mxMatrix(typ='Diag',nrow=ny,ncol=ny,fr=TRUE,valu=5,nam='sds') Man1=mxMatrix(typ='Full',nrow=1,ncol=ny,fr=TRUE,valu=25,nam='man1') MkS1=mxAlgbra(xprssion=sds%*%cor%*%sds,nam='Ssat1') # 12x12 corrlation matrix # 12x12 diagonal matrix (st dvs) # 1x12 vctor mans # xpctd covarianc matrix satmodls1=mxmodl('part1',rs1, Sds1, Man1,MkS1) # assmbl th modl # data + stimation function satdats1=mxmodl("part2", mxdata( obsrvd=datb2, typ="raw"), # th data mxexpctationnormal( covarianc="part1.ssat1", mans="part1.man1", dimnams=varnams), # th fit function mxfitfunctionml() ) # data & xpctd cov/mans # fit th saturatd modl... Modls1 <- mxmodl("modls1", satmodls1, satdats1, mxalgbra(part2.objctiv, nam="minus2logliklihood"), mxfitfunctionalgbra("minus2logliklihood")) Modls1_out <- mxrun(modls1) 10

11 CFA > summary(modls1_out) obsrvd statistics: 4332 stimatd paramtrs: 90 dgrs of frdom: log liklihood: numbr of obsrvations:

12 CFA S y = Y t + Q 1 1 N E n1 n2 n3 n4 n5 n

13 Ly=mxMatrix(typ='Full',nrow=ny,ncol=n, fr=matrix(c( T,F, T,F, T,F, T,F, T,F, T,F, F,T, F,T, F,T, F,T, F,T, F,T),ny,n,byrow=T), valus=c(4,4,4,4,4,4,0,0,0,0,0,0, 0,0,0,0,0,0,4,4,4,4,4,4), # rad colunm-wis labls=matrix(c( 'f1_1','f1_2', 'f2_1','f2_2', 'f3_1','f3_2', 'f4_1','f4_2', 'f5_1','f5_2', 'f6_1','f6_2', 'f7_1','f7_2', 'f8_1','f8_2', 'f9_1','f9_2', 'f10_1','f10_2', 'f11_1','f11_2', 'f12_1','f12_2'),ny,n,byrow=t),nam='ly'), n1 n2 n3 Dfin factor loading matrix S y = Y t + Q 1 1 N n4 n5 n E CFA

14 CFA T=mxMatrix(typ='Diag',nrow=ny,ncol=ny, labls=c('rn1','rn2','rn3','rn4','rn5','rn6', 'r1','r2','r3','r4','r5','r6'), fr=true,valu=10,nam='t') Dfin covarianc matrix of rsiduals Q S y = Y t + Q 1 1 N E n1 n2 n3 n4 n5 n

15 CFA ## latnt corrlation matrix Ps=mxMatrix(typ='Symm',nrow=n,ncol=n, fr=c(false,true,false), labls=c('v1_0','r12_0', 'v2_0'), valus=c(1,.0,1),nam='ps') Dfin covarianc matrix of factors Y S y = Y t + Q NOTE: scaling of th common factors by fixing th variancs to qual 1. Y is a corrlation matrix! 1 1 N E n1 n2 n3 n4 n5 n

16 CFA # mans Tau=mxMatrix(typ='Full',nrow=1,ncol=ny,fr=TRUE,valu=25, labls=c('mn1','mn2','mn3','mn4','mn5','mn6', 'm1','m2','m3','m4','m5','m6'), nam='mans') 16

17 CFA S y = Y t + Q MKS=mxAlgbra(xprssion=Ly%*%(Ps)%*%t(Ly)+ T%*%t(T),nam='Sigma'), MKM=mxAlgbra(xprssion=Mans,nam='mans')... assmbl th modl and run CFM1_out = mxrun(cfamodl1) 1 1 N E n1 n2 n3 n4 n5 n

18 CFA > mxcompar(modls1_out,cfm1_out) bas comparison p minus2ll df AIC diffll diffdf p 1 modls1 <NA> NA NA NA 2 modls1 CFM This modl dos not fit (but w alrady know this from th EFA rsults). 18

19 CFA > print(round(fa1s1pr$loadings[,1:2],4)) Factor1 Factor2 n n n n n n > round(st_ly,3) [,1] [,2] [1,] [2,] [3,] [4,] [5,] [6,] [7,] [8,] [9,] [10,] [11,] [12,] To do: fr th cross loadings Ly[5,2] and Ly[9,1] 19

20 To do: fr th cross loadings Ly[5,2] and Ly[9,1] CFA r 1 1 N E n1 n2 n3 n4 n5 n Tst whthr th cross loadings ar unqual zro using th liklihood ratio tst 20

21 CFA bas comparison p minus2ll df AIC diffll diffdf p 1 CFM1 <NA> NA NA NA 2 CFM1 CFM Givn alpha=.05, w rjct th hypothsis Ly[5,2] = Ly[9,1] = 0 Thrfor ithr on or both ar not qual to zro. print(round(st_ly,2)) [,1] [,2] [1,] [2,] [3,] [4,] [5,] [6,] [7,] [8,] [9,] [10,] [11,] [12,]

22 Inspct th corrlation matrix Y CFA N E n1 n2 n3 n4 n5 n

23 Rliability of th indicators > round(diag(sfit1)/diag(sfit),2) [1] r 1 1 N E st_ly=mxeval(cfamodl1.ly,cfm2_out) st_ps=mxeval(cfamodl1.ps,cfm2_out) st_t=mxeval(cfamodl1.t,cfm2_out) st_t=st_t^2 Sfit1=st_Ly%*%st_Ps%*%t(st_Ly) Sfit=Sfit1+st_T rl=diag(sfit1)/diag(sfit) print(round(rl,3)) Varianc of n1 du to N dividd by th total varianc of n1: CFA / =.74. n1 n2 n3 n4 n5 n Th common factor N xplains 74% of th varianc in itm 1. 23

Chapter 13 GMM for Linear Factor Models in Discount Factor form. GMM on the pricing errors gives a crosssectional

Chapter 13 GMM for Linear Factor Models in Discount Factor form. GMM on the pricing errors gives a crosssectional Chaptr 13 GMM for Linar Factor Modls in Discount Factor form GMM on th pricing rrors givs a crosssctional rgrssion h cas of xcss rturns Hors rac sting for charactristic sting for pricd factors: lambdas

More information

Solution of Assignment #2

Solution of Assignment #2 olution of Assignmnt #2 Instructor: Alirza imchi Qustion #: For simplicity, assum that th distribution function of T is continuous. Th distribution function of R is: F R ( r = P( R r = P( log ( T r = P(log

More information

Applied Statistics II - Categorical Data Analysis Data analysis using Genstat - Exercise 2 Logistic regression

Applied Statistics II - Categorical Data Analysis Data analysis using Genstat - Exercise 2 Logistic regression Applid Statistics II - Catgorical Data Analysis Data analysis using Gnstat - Exrcis 2 Logistic rgrssion Analysis 2. Logistic rgrssion for a 2 x k tabl. Th tabl blow shows th numbr of aphids aliv and dad

More information

Observer Bias and Reliability By Xunchi Pu

Observer Bias and Reliability By Xunchi Pu Obsrvr Bias and Rliability By Xunchi Pu Introduction Clarly all masurmnts or obsrvations nd to b mad as accuratly as possibl and invstigators nd to pay carful attntion to chcking th rliability of thir

More information

Chapter 6 Student Lecture Notes 6-1

Chapter 6 Student Lecture Notes 6-1 Chaptr 6 Studnt Lctur Nots 6-1 Chaptr Goals QM353: Busnss Statstcs Chaptr 6 Goodnss-of-Ft Tsts and Contngncy Analyss Aftr compltng ths chaptr, you should b abl to: Us th ch-squar goodnss-of-ft tst to dtrmn

More information

AS 5850 Finite Element Analysis

AS 5850 Finite Element Analysis AS 5850 Finit Elmnt Analysis Two-Dimnsional Linar Elasticity Instructor Prof. IIT Madras Equations of Plan Elasticity - 1 displacmnt fild strain- displacmnt rlations (infinitsimal strain) in matrix form

More information

Review Statistics review 14: Logistic regression Viv Bewick 1, Liz Cheek 1 and Jonathan Ball 2

Review Statistics review 14: Logistic regression Viv Bewick 1, Liz Cheek 1 and Jonathan Ball 2 Critical Car Fbruary 2005 Vol 9 No 1 Bwick t al. Rviw Statistics rviw 14: Logistic rgrssion Viv Bwick 1, Liz Chk 1 and Jonathan Ball 2 1 Snior Lcturr, School of Computing, Mathmatical and Information Scincs,

More information

ECE602 Exam 1 April 5, You must show ALL of your work for full credit.

ECE602 Exam 1 April 5, You must show ALL of your work for full credit. ECE62 Exam April 5, 27 Nam: Solution Scor: / This xam is closd-book. You must show ALL of your work for full crdit. Plas rad th qustions carfully. Plas chck your answrs carfully. Calculators may NOT b

More information

TEMASEK JUNIOR COLLEGE, SINGAPORE. JC 2 Preliminary Examination 2017

TEMASEK JUNIOR COLLEGE, SINGAPORE. JC 2 Preliminary Examination 2017 TEMASEK JUNIOR COLLEGE, SINGAPORE JC Prliminary Eamination 7 MATHEMATICS 886/ Highr 9 August 7 Additional Matrials: Answr papr hours List of Formula (MF6) READ THESE INSTRUCTIONS FIRST Writ your Civics

More information

Errata. Items with asterisks will still be in the Second Printing

Errata. Items with asterisks will still be in the Second Printing Errata Itms with astrisks will still b in th Scond Printing Author wbsit URL: http://chs.unl.du/edpsych/rjsit/hom. P7. Th squar root of rfrrd to σ E (i.., σ E is rfrrd to not Th squar root of σ E (i..,

More information

Laboratory work # 8 (14) EXPERIMENTAL ESTIMATION OF CRITICAL STRESSES IN STRINGER UNDER COMPRESSION

Laboratory work # 8 (14) EXPERIMENTAL ESTIMATION OF CRITICAL STRESSES IN STRINGER UNDER COMPRESSION Laboratory wor # 8 (14) XPRIMNTAL STIMATION OF CRITICAL STRSSS IN STRINGR UNDR COMPRSSION At action of comprssing ffort on a bar (column, rod, and stringr) two inds of loss of stability ar possibl: 1)

More information

Nonparametric Methods: Goodness-of-Fit Tests

Nonparametric Methods: Goodness-of-Fit Tests Nonparamtric Mthods: Goodnss-o-Fit Tsts Chaptr 15 McGraw-Hill/Irwin Copyright 013 by Th McGraw-Hill Companis, Inc. All rights rsrvd. LEARNING OBJECTIVES LO 15-1 Conduct a tst o hypothsis comparing an obsrvd

More information

APPENDIX: STATISTICAL TOOLS

APPENDIX: STATISTICAL TOOLS I. Nots o radom samplig Why do you d to sampl radomly? APPENDI: STATISTICAL TOOLS I ordr to masur som valu o a populatio of orgaisms, you usually caot masur all orgaisms, so you sampl a subst of th populatio.

More information

Linear Non-Gaussian Structural Equation Models

Linear Non-Gaussian Structural Equation Models IMPS 8, Durham, NH Linar Non-Gaussian Structural Equation Modls Shohi Shimizu, Patrik Hoyr and Aapo Hyvarinn Osaka Univrsity, Japan Univrsity of Hlsinki, Finland Abstract Linar Structural Equation Modling

More information

Estimation of apparent fraction defective: A mathematical approach

Estimation of apparent fraction defective: A mathematical approach Availabl onlin at www.plagiarsarchlibrary.com Plagia Rsarch Library Advancs in Applid Scinc Rsarch, 011, (): 84-89 ISSN: 0976-8610 CODEN (USA): AASRFC Estimation of apparnt fraction dfctiv: A mathmatical

More information

Cramér-Rao Inequality: Let f(x; θ) be a probability density function with continuous parameter

Cramér-Rao Inequality: Let f(x; θ) be a probability density function with continuous parameter WHEN THE CRAMÉR-RAO INEQUALITY PROVIDES NO INFORMATION STEVEN J. MILLER Abstract. W invstigat a on-paramtr family of probability dnsitis (rlatd to th Parto distribution, which dscribs many natural phnomna)

More information

That is, we start with a general matrix: And end with a simpler matrix:

That is, we start with a general matrix: And end with a simpler matrix: DIAGON ALIZATION OF THE STR ESS TEN SOR INTRO DUCTIO N By th us of Cauchy s thorm w ar abl to rduc th numbr of strss componnts in th strss tnsor to only nin valus. An additional simplification of th strss

More information

Text: WMM, Chapter 5. Sections , ,

Text: WMM, Chapter 5. Sections , , Lcturs 6 - Continuous Probabilit Distributions Tt: WMM, Chaptr 5. Sctions 6.-6.4, 6.6-6.8, 7.-7. In th prvious sction, w introduc som of th common probabilit distribution functions (PDFs) for discrt sampl

More information

Abstract Interpretation: concrete and abstract semantics

Abstract Interpretation: concrete and abstract semantics Abstract Intrprtation: concrt and abstract smantics Concrt smantics W considr a vry tiny languag that manags arithmtic oprations on intgrs valus. Th (concrt) smantics of th languags cab b dfind by th funzcion

More information

Self-Adjointness and Its Relationship to Quantum Mechanics. Ronald I. Frank 2016

Self-Adjointness and Its Relationship to Quantum Mechanics. Ronald I. Frank 2016 Ronald I. Frank 06 Adjoint https://n.wikipdia.org/wiki/adjoint In gnral thr is an oprator and a procss that dfin its adjoint *. It is thn slf-adjoint if *. Innr product spac https://n.wikipdia.org/wiki/innr_product_spac

More information

Chapter 3 Lecture 14 Longitudinal stick free static stability and control 3 Topics

Chapter 3 Lecture 14 Longitudinal stick free static stability and control 3 Topics Chaptr 3 Lctur 14 Longitudinal stick fr static stability and control 3 Topics 3.4.4 Rquirmnt for propr stick forc variation 3.4.5 Fl of th stability lvl by th pilot Exampl 3.3 3.5 Dtrmination of stick-fr

More information

Introduction to Condensed Matter Physics

Introduction to Condensed Matter Physics Introduction to Condnsd Mattr Physics pcific hat M.P. Vaughan Ovrviw Ovrviw of spcific hat Hat capacity Dulong-Ptit Law Einstin modl Dby modl h Hat Capacity Hat capacity h hat capacity of a systm hld at

More information

Introduction to Arithmetic Geometry Fall 2013 Lecture #20 11/14/2013

Introduction to Arithmetic Geometry Fall 2013 Lecture #20 11/14/2013 18.782 Introduction to Arithmtic Gomtry Fall 2013 Lctur #20 11/14/2013 20.1 Dgr thorm for morphisms of curvs Lt us rstat th thorm givn at th nd of th last lctur, which w will now prov. Thorm 20.1. Lt φ:

More information

Constants and Conversions:

Constants and Conversions: EXAM INFORMATION Radial Distribution Function: P 2 ( r) RDF( r) Br R( r ) 2, B is th normalization constant. Ordr of Orbital Enrgis: Homonuclar Diatomic Molculs * * * * g1s u1s g 2s u 2s u 2 p g 2 p g

More information

Recursive Estimation of Dynamic Time-Varying Demand Models

Recursive Estimation of Dynamic Time-Varying Demand Models Intrnational Confrnc on Computr Systms and chnologis - CompSysch 06 Rcursiv Estimation of Dynamic im-varying Dmand Modls Alxandr Efrmov Abstract: h papr prsnts an implmntation of a st of rcursiv algorithms

More information

ph People Grade Level: basic Duration: minutes Setting: classroom or field site

ph People Grade Level: basic Duration: minutes Setting: classroom or field site ph Popl Adaptd from: Whr Ar th Frogs? in Projct WET: Curriculum & Activity Guid. Bozman: Th Watrcours and th Council for Environmntal Education, 1995. ph Grad Lvl: basic Duration: 10 15 minuts Stting:

More information

INFLUENCE OF GROUND SUBSIDENCE IN THE DAMAGE TO MEXICO CITY S PRIMARY WATER SYSTEM DUE TO THE 1985 EARTHQUAKE

INFLUENCE OF GROUND SUBSIDENCE IN THE DAMAGE TO MEXICO CITY S PRIMARY WATER SYSTEM DUE TO THE 1985 EARTHQUAKE 13 th World Confrnc on Earthquak Enginring Vancouvr, B.C., Canada August 1-6, 2004 Papr No. 2165 INFLUENCE OF GROUND SUBSIDENCE IN THE DAMAGE TO MEXICO CITY S PRIMARY WATER SYSTEM DUE TO THE 1985 EARTHQUAKE

More information

UNTYPED LAMBDA CALCULUS (II)

UNTYPED LAMBDA CALCULUS (II) 1 UNTYPED LAMBDA CALCULUS (II) RECALL: CALL-BY-VALUE O.S. Basic rul Sarch ruls: (\x.) v [v/x] 1 1 1 1 v v CALL-BY-VALUE EVALUATION EXAMPLE (\x. x x) (\y. y) x x [\y. y / x] = (\y. y) (\y. y) y [\y. y /

More information

Content Skills Assessments Lessons. Identify, classify, and apply properties of negative and positive angles.

Content Skills Assessments Lessons. Identify, classify, and apply properties of negative and positive angles. Tachr: CORE TRIGONOMETRY Yar: 2012-13 Cours: TRIGONOMETRY Month: All Months S p t m b r Angls Essntial Qustions Can I idntify draw ngativ positiv angls in stard position? Do I hav a working knowldg of

More information

Differentiation of Exponential Functions

Differentiation of Exponential Functions Calculus Modul C Diffrntiation of Eponntial Functions Copyright This publication Th Northrn Albrta Institut of Tchnology 007. All Rights Rsrvd. LAST REVISED March, 009 Introduction to Diffrntiation of

More information

2. Laser physics - basics

2. Laser physics - basics . Lasr physics - basics Spontanous and stimulatd procsss Einstin A and B cofficints Rat quation analysis Gain saturation What is a lasr? LASER: Light Amplification by Stimulatd Emission of Radiation "light"

More information

Exam 1. It is important that you clearly show your work and mark the final answer clearly, closed book, closed notes, no calculator.

Exam 1. It is important that you clearly show your work and mark the final answer clearly, closed book, closed notes, no calculator. Exam N a m : _ S O L U T I O N P U I D : I n s t r u c t i o n s : It is important that you clarly show your work and mark th final answr clarly, closd book, closd nots, no calculator. T i m : h o u r

More information

Analyzing Frequencies

Analyzing Frequencies Frquncy (# ndvduals) Frquncy (# ndvduals) /3/16 H o : No dffrnc n obsrvd sz frquncs and that prdctd by growth modl How would you analyz ths data? 15 Obsrvd Numbr 15 Expctd Numbr from growth modl 1 1 5

More information

Search for the Dark Photon at Belle for 0.27 < m A < 3 GeV/c 2

Search for the Dark Photon at Belle for 0.27 < m A < 3 GeV/c 2 Sarch for th Dark Photon at Bll for.7 < m A < 3 GV/c Igal Jagl Univrsity of Hawai i at Mānoa for th Bll Collaboration QCD, Montpllir, 6 Juillt BELLE Igal Jagl (UH) Sarch for th Dark photon at Bll QCD /

More information

VII. Quantum Entanglement

VII. Quantum Entanglement VII. Quantum Entanglmnt Quantum ntanglmnt is a uniqu stat of quantum suprposition. It has bn studid mainly from a scintific intrst as an vidnc of quantum mchanics. Rcntly, it is also bing studid as a basic

More information

Elements of Statistical Thermodynamics

Elements of Statistical Thermodynamics 24 Elmnts of Statistical Thrmodynamics Statistical thrmodynamics is a branch of knowldg that has its own postulats and tchniqus. W do not attmpt to giv hr vn an introduction to th fild. In this chaptr,

More information

Extracting Common Factors to Classify Companies Listed in the Stock Exchange of Thailand by Using an Accounting Based Model

Extracting Common Factors to Classify Companies Listed in the Stock Exchange of Thailand by Using an Accounting Based Model Extracting Common Factors to Classify Companis Listd in th Stock Exchang of Thailand by Using an Accounting Basd Modl Krisada Khruachal Faculty of Businss Administration, Siam Tchnology Collg, Thailand

More information

Appendix. Kalman Filter

Appendix. Kalman Filter Appndix A Kalman Filtr OPTIMAL stimation thory has a vry broad rang of applications which vary from stimation of rivr ows to satllit orbit stimation and nuclar ractor paramtr idntication. In this appndix

More information

CS 361 Meeting 12 10/3/18

CS 361 Meeting 12 10/3/18 CS 36 Mting 2 /3/8 Announcmnts. Homwork 4 is du Friday. If Friday is Mountain Day, homwork should b turnd in at my offic or th dpartmnt offic bfor 4. 2. Homwork 5 will b availabl ovr th wknd. 3. Our midtrm

More information

Probability and Stochastic Processes: A Friendly Introduction for Electrical and Computer Engineers Roy D. Yates and David J.

Probability and Stochastic Processes: A Friendly Introduction for Electrical and Computer Engineers Roy D. Yates and David J. Probability and Stochastic Procsss: A Frindly Introduction for Elctrical and Computr Enginrs Roy D. Yats and David J. Goodman Problm Solutions : Yats and Goodman,4.3. 4.3.4 4.3. 4.4. 4.4.4 4.4.6 4.. 4..7

More information

Schrodinger Equation in 3-d

Schrodinger Equation in 3-d Schrodingr Equation in 3-d ψ( xyz,, ) ψ( xyz,, ) ψ( xyz,, ) + + + Vxyz (,, ) ψ( xyz,, ) = Eψ( xyz,, ) m x y z p p p x y + + z m m m + V = E p m + V = E E + k V = E Infinit Wll in 3-d V = x > L, y > L,

More information

Sara Godoy del Olmo Calculation of contaminated soil volumes : Geostatistics applied to a hydrocarbons spill Lac Megantic Case

Sara Godoy del Olmo Calculation of contaminated soil volumes : Geostatistics applied to a hydrocarbons spill Lac Megantic Case wwwnvisol-canadaca Sara Godoy dl Olmo Calculation of contaminatd soil volums : Gostatistics applid to a hydrocarbons spill Lac Mgantic Cas Gostatistics: study of a PH contamination CONTEXT OF THE STUDY

More information

What are those βs anyway? Understanding Design Matrix & Odds ratios

What are those βs anyway? Understanding Design Matrix & Odds ratios Ral paramtr stimat WILD 750 - Wildlif Population Analysis of 6 What ar thos βs anyway? Undrsting Dsign Matrix & Odds ratios Rfrncs Hosmr D.W.. Lmshow. 000. Applid logistic rgrssion. John Wily & ons Inc.

More information

Self-interaction mass formula that relates all leptons and quarks to the electron

Self-interaction mass formula that relates all leptons and quarks to the electron Slf-intraction mass formula that rlats all lptons and quarks to th lctron GERALD ROSEN (a) Dpartmnt of Physics, Drxl Univrsity Philadlphia, PA 19104, USA PACS. 12.15. Ff Quark and lpton modls spcific thoris

More information

Eigenvalue Distributions of Quark Matrix at Finite Isospin Chemical Potential

Eigenvalue Distributions of Quark Matrix at Finite Isospin Chemical Potential Tim: Tusday, 5: Room: Chsapak A Eignvalu Distributions of Quark Matri at Finit Isospin Chmical Potntial Prsntr: Yuji Sasai Tsuyama National Collg of Tchnology Co-authors: Grnot Akmann, Atsushi Nakamura

More information

MULTIVARIATE BAYESIAN REGRESSION ANALYSIS APPLIED TO PSEUDO-ACCELERATION ATENUATTION RELATIONSHIPS

MULTIVARIATE BAYESIAN REGRESSION ANALYSIS APPLIED TO PSEUDO-ACCELERATION ATENUATTION RELATIONSHIPS h 4 th World Confrnc on Earthquak Enginring Octobr -7, 8, Bijing, China MULIVARIAE BAYESIAN REGRESSION ANALYSIS APPLIED O PSEUDO-ACCELERAION AENUAION RELAIONSHIPS D. Arroyo and M. Ordaz Associat Profssor,

More information

General Notes About 2007 AP Physics Scoring Guidelines

General Notes About 2007 AP Physics Scoring Guidelines AP PHYSICS C: ELECTRICITY AND MAGNETISM 2007 SCORING GUIDELINES Gnral Nots About 2007 AP Physics Scoring Guidlins 1. Th solutions contain th most common mthod of solving th fr-rspons qustions and th allocation

More information

Hydrogen Atom and One Electron Ions

Hydrogen Atom and One Electron Ions Hydrogn Atom and On Elctron Ions Th Schrödingr quation for this two-body problm starts out th sam as th gnral two-body Schrödingr quation. First w sparat out th motion of th cntr of mass. Th intrnal potntial

More information

15. Stress-Strain behavior of soils

15. Stress-Strain behavior of soils 15. Strss-Strain bhavior of soils Sand bhavior Usually shard undr draind conditions (rlativly high prmability mans xcss por prssurs ar not gnratd). Paramtrs govrning sand bhaviour is: Rlativ dnsity Effctiv

More information

4.2 Design of Sections for Flexure

4.2 Design of Sections for Flexure 4. Dsign of Sctions for Flxur This sction covrs th following topics Prliminary Dsign Final Dsign for Typ 1 Mmbrs Spcial Cas Calculation of Momnt Dmand For simply supportd prstrssd bams, th maximum momnt

More information

SME 3033 FINITE ELEMENT METHOD. Bending of Prismatic Beams (Initial notes designed by Dr. Nazri Kamsah)

SME 3033 FINITE ELEMENT METHOD. Bending of Prismatic Beams (Initial notes designed by Dr. Nazri Kamsah) Bnding of Prismatic Bams (Initia nots dsignd by Dr. Nazri Kamsah) St I-bams usd in a roof construction. 5- Gnra Loading Conditions For our anaysis, w wi considr thr typs of oading, as iustratd bow. Not:

More information

Panel Data Analysis Introduction

Panel Data Analysis Introduction Panl Data Analysis Introduction Modl Rprsntation N-first or T-first rprsntation Poold Modl Fixd Effcts Modl Random Effcts Modl Asymptotic Thory N, or T N, T Panl-Robust Infrnc Panl Data Analysis Introduction

More information

EXST Regression Techniques Page 1

EXST Regression Techniques Page 1 EXST704 - Rgrssion Tchniqus Pag 1 Masurmnt rrors in X W hav assumd that all variation is in Y. Masurmnt rror in this variabl will not ffct th rsults, as long as thy ar uncorrlatd and unbiasd, sinc thy

More information

Linear-Phase FIR Transfer Functions. Functions. Functions. Functions. Functions. Functions. Let

Linear-Phase FIR Transfer Functions. Functions. Functions. Functions. Functions. Functions. Let It is impossibl to dsign an IIR transfr function with an xact linar-phas It is always possibl to dsign an FIR transfr function with an xact linar-phas rspons W now dvlop th forms of th linarphas FIR transfr

More information

Mor Tutorial at www.dumblittldoctor.com Work th problms without a calculator, but us a calculator to chck rsults. And try diffrntiating your answrs in part III as a usful chck. I. Applications of Intgration

More information

Class(ic) Scorecards

Class(ic) Scorecards Class(ic) Scorcards Slcting Charactristics and Attributs in Logistic Rgrssion Edinburgh Crdit Scoring Confrnc - 25 August 2011 Grard Scallan grard.scallan@scorplus.com 1 Class(ic) Scorcards Using th Statistics!

More information

There is an arbitrary overall complex phase that could be added to A, but since this makes no difference we set it to zero and choose A real.

There is an arbitrary overall complex phase that could be added to A, but since this makes no difference we set it to zero and choose A real. Midtrm #, Physics 37A, Spring 07. Writ your rsponss blow or on xtra pags. Show your work, and tak car to xplain what you ar doing; partial crdit will b givn for incomplt answrs that dmonstrat som concptual

More information

DIFFERENTIAL EQUATION

DIFFERENTIAL EQUATION MD DIFFERENTIAL EQUATION Sllabus : Ordinar diffrntial quations, thir ordr and dgr. Formation of diffrntial quations. Solution of diffrntial quations b th mthod of sparation of variabls, solution of homognous

More information

Coupled Pendulums. Two normal modes.

Coupled Pendulums. Two normal modes. Tim Dpndnt Two Stat Problm Coupld Pndulums Wak spring Two normal mods. No friction. No air rsistanc. Prfct Spring Start Swinging Som tim latr - swings with full amplitud. stationary M +n L M +m Elctron

More information

Data Assimilation 1. Alan O Neill National Centre for Earth Observation UK

Data Assimilation 1. Alan O Neill National Centre for Earth Observation UK Data Assimilation 1 Alan O Nill National Cntr for Earth Obsrvation UK Plan Motivation & basic idas Univariat (scalar) data assimilation Multivariat (vctor) data assimilation 3d-Variational Mthod (& optimal

More information

Ch. 24 Molecular Reaction Dynamics 1. Collision Theory

Ch. 24 Molecular Reaction Dynamics 1. Collision Theory Ch. 4 Molcular Raction Dynamics 1. Collision Thory Lctur 16. Diffusion-Controlld Raction 3. Th Matrial Balanc Equation 4. Transition Stat Thory: Th Eyring Equation 5. Transition Stat Thory: Thrmodynamic

More information

Electromagnetics Research Group A THEORETICAL MODEL OF A LOSSY DIELECTRIC SLAB FOR THE CHARACTERIZATION OF RADAR SYSTEM PERFORMANCE SPECIFICATIONS

Electromagnetics Research Group A THEORETICAL MODEL OF A LOSSY DIELECTRIC SLAB FOR THE CHARACTERIZATION OF RADAR SYSTEM PERFORMANCE SPECIFICATIONS Elctromagntics Rsarch Group THEORETICL MODEL OF LOSSY DIELECTRIC SLB FOR THE CHRCTERIZTION OF RDR SYSTEM PERFORMNCE SPECIFICTIONS G.L. Charvat, Prof. Edward J. Rothwll Michigan Stat Univrsit 1 Ovrviw of

More information

Prediction of the Pressure Signature of a Ship in a Seastate

Prediction of the Pressure Signature of a Ship in a Seastate Prdiction of th Prssur Signatur of a Ship in a Sastat Mark Hyman Thai Nguyn, Knnard Watson Hydromchanics Group, Cod R11 NSWC Coastal Systms Station Panama City, FL 347 phon: (85) 34-416 fax: (85) 35-5374

More information

Logistic, Poisson, and Nonlinear Regression Problems

Logistic, Poisson, and Nonlinear Regression Problems Logistic, Poisson, and Nonlinar Rgrssion Problms Logistic Rgrssion QA.1. A study is conductd to masur th ffcts of lvls of a hrbicid on th probability of dath for 2 wd varitis. Th rsarchr slcts 5 dosag

More information

AP Biology Lab 7 GENETICS OF ORGANISMS

AP Biology Lab 7 GENETICS OF ORGANISMS AP Biology Laboratory Dat: Nam and Priod: OVERVIEW AP Biology Lab 7 GENETICS OF ORGANISMS In this lab you will us living organism to do gntic cross3s. You will larn how to collct and manipulat th organisms,

More information

Properties of Phase Space Wavefunctions and Eigenvalue Equation of Momentum Dispersion Operator

Properties of Phase Space Wavefunctions and Eigenvalue Equation of Momentum Dispersion Operator Proprtis of Phas Spac Wavfunctions and Eignvalu Equation of Momntum Disprsion Oprator Ravo Tokiniaina Ranaivoson 1, Raolina Andriambololona 2, Hanitriarivo Rakotoson 3 raolinasp@yahoo.fr 1 ;jacqulinraolina@hotmail.com

More information

ST 524 NCSU - Fall 2008 One way Analysis of variance Variances not homogeneous

ST 524 NCSU - Fall 2008 One way Analysis of variance Variances not homogeneous ST 54 NCSU - Fall 008 On way Analyss of varanc Varancs not homognous On way Analyss of varanc Exampl (Yandll, 997) A plant scntst masurd th concntraton of a partcular vrus n plant sap usng ELISA (nzym-lnkd

More information

Analyzing genotype-by-environment interaction using curvilinear regression

Analyzing genotype-by-environment interaction using curvilinear regression 357 Analyzing Scintia GE Agricola using curvilinar rgrssion Analyzing gnotyp-by-nvironmnt intraction using curvilinar rgrssion Dulc Gamito Santinhos Prira 1 * Paulo Canas Rodrigus 24 wona Mjza 3 Stanislaw

More information

Properties of Quarks ( ) Isospin. π = 1, 1

Properties of Quarks ( ) Isospin. π = 1, 1 Proprtis of Quarks Isospin So far, w hav discussd thr familis of lptons but principally concntratd on on doublt of quarks, th u and d. W will now introduc othr typs of quarks, along with th nw quantum

More information

Search sequence databases 3 10/25/2016

Search sequence databases 3 10/25/2016 Sarch squnc databass 3 10/25/2016 Etrm valu distribution Ø Suppos X is a random variabl with probability dnsity function p(, w sampl a larg numbr S of indpndnt valus of X from this distribution for an

More information

Dealing with quantitative data and problem solving life is a story problem! Attacking Quantitative Problems

Dealing with quantitative data and problem solving life is a story problem! Attacking Quantitative Problems Daling with quantitati data and problm soling lif is a story problm! A larg portion of scinc inols quantitati data that has both alu and units. Units can sa your butt! Nd handl on mtric prfixs Dimnsional

More information

ARIMA Methods of Detecting Outliers in Time Series Periodic Processes

ARIMA Methods of Detecting Outliers in Time Series Periodic Processes Articl Intrnational Journal of Modrn Mathmatical Scincs 014 11(1): 40-48 Intrnational Journal of Modrn Mathmatical Scincs Journal hompag:www.modrnscintificprss.com/journals/ijmms.aspx ISSN:166-86X Florida

More information

A Propagating Wave Packet Group Velocity Dispersion

A Propagating Wave Packet Group Velocity Dispersion Lctur 8 Phys 375 A Propagating Wav Packt Group Vlocity Disprsion Ovrviw and Motivation: In th last lctur w lookd at a localizd solution t) to th 1D fr-particl Schrödingr quation (SE) that corrsponds to

More information

Pipe flow friction, small vs. big pipes

Pipe flow friction, small vs. big pipes Friction actor (t/0 t o pip) Friction small vs larg pips J. Chaurtt May 016 It is an intrsting act that riction is highr in small pips than largr pips or th sam vlocity o low and th sam lngth. Friction

More information

The general linear model for fmri

The general linear model for fmri h gnral linar modl for fmri Mthods and Modls in fmri, 7.0.07 Jakob Hinzl hinzl@biomd..thz.ch ranslational Nuromodling Unit (NU) Institut for Biomdical Enginring (IB) Univrsit and EH Zürich Man thanks to

More information

LEP Higgs Search Results. Chris Tully Weak Interactions and Neutrinos Workshop January 21-26, 2002

LEP Higgs Search Results. Chris Tully Weak Interactions and Neutrinos Workshop January 21-26, 2002 LEP Higgs Sarch Rsults Chris Tully Wak Intractions and Nutrinos Workshop January 1-6, 00 Viw of LEP Exprimnts Spontanous Symmtry Braking Higgs Fild 1-complx doublt Hyprcharg(Y) U(1) Y Symmtry Lft-Handd

More information

Week 3: Connected Subgraphs

Week 3: Connected Subgraphs Wk 3: Connctd Subgraphs Sptmbr 19, 2016 1 Connctd Graphs Path, Distanc: A path from a vrtx x to a vrtx y in a graph G is rfrrd to an xy-path. Lt X, Y V (G). An (X, Y )-path is an xy-path with x X and y

More information

What does the data look like? Logistic Regression. How can we apply linear model to categorical data like this? Linear Probability Model

What does the data look like? Logistic Regression. How can we apply linear model to categorical data like this? Linear Probability Model Logistic Rgrssion Almost sam to rgrssion ct that w hav binary rsons variabl (ys or no, succss or failur) W assum that th qualitativ rsons variabl has a discrt distribution lik binomial. What dos th data

More information

Full Waveform Inversion Using an Energy-Based Objective Function with Efficient Calculation of the Gradient

Full Waveform Inversion Using an Energy-Based Objective Function with Efficient Calculation of the Gradient Full Wavform Invrsion Using an Enrgy-Basd Objctiv Function with Efficint Calculation of th Gradint Itm yp Confrnc Papr Authors Choi, Yun Sok; Alkhalifah, ariq Ali Citation Choi Y, Alkhalifah (217) Full

More information

Lie Groups HW7. Wang Shuai. November 2015

Lie Groups HW7. Wang Shuai. November 2015 Li roups HW7 Wang Shuai Novmbr 015 1 Lt (π, V b a complx rprsntation of a compact group, show that V has an invariant non-dgnratd Hrmitian form. For any givn Hrmitian form on V, (for xampl (u, v = i u

More information

Robust Regression. Appendix to An R and S-PLUS Companion to Applied Regression. John Fox. January 2002

Robust Regression. Appendix to An R and S-PLUS Companion to Applied Regression. John Fox. January 2002 Robust Rgrssion Appndix to An R and S-PLUS Companion to Applid Rgrssion John Fox January 2002 1 M-Estimation Linar last-squars stimats can bhav badly whn th rror distribution is not normal, particularly

More information

22/ Breakdown of the Born-Oppenheimer approximation. Selection rules for rotational-vibrational transitions. P, R branches.

22/ Breakdown of the Born-Oppenheimer approximation. Selection rules for rotational-vibrational transitions. P, R branches. Subjct Chmistry Papr No and Titl Modul No and Titl Modul Tag 8/ Physical Spctroscopy / Brakdown of th Born-Oppnhimr approximation. Slction ruls for rotational-vibrational transitions. P, R branchs. CHE_P8_M

More information

Chapter 10. The Chi-Squared Test

Chapter 10. The Chi-Squared Test Chaptr 10 Th Chi-Squard Tst Forty for you, sixty for m And qual partnrs w will b. Grald Barzan Th Avonford Star Choos your hospital ward carfully! Statistics rlasd today as part of th nw govrnmnt initiativ

More information

Chapter 13 Aggregate Supply

Chapter 13 Aggregate Supply Chaptr 13 Aggrgat Supply 0 1 Larning Objctivs thr modls of aggrgat supply in which output dpnds positivly on th pric lvl in th short run th short-run tradoff btwn inflation and unmploymnt known as th Phillips

More information

Diploma Macro Paper 2

Diploma Macro Paper 2 Diploma Macro Papr 2 Montary Macroconomics Lctur 6 Aggrgat supply and putting AD and AS togthr Mark Hays 1 Exognous: M, G, T, i*, π Goods markt KX and IS (Y, C, I) Mony markt (LM) (i, Y) Labour markt (P,

More information

MCB137: Physical Biology of the Cell Spring 2017 Homework 6: Ligand binding and the MWC model of allostery (Due 3/23/17)

MCB137: Physical Biology of the Cell Spring 2017 Homework 6: Ligand binding and the MWC model of allostery (Due 3/23/17) MCB37: Physical Biology of th Cll Spring 207 Homwork 6: Ligand binding and th MWC modl of allostry (Du 3/23/7) Hrnan G. Garcia March 2, 207 Simpl rprssion In class, w drivd a mathmatical modl of how simpl

More information

Additional Math (4047) Paper 2 (100 marks) y x. 2 d. d d

Additional Math (4047) Paper 2 (100 marks) y x. 2 d. d d Aitional Math (07) Prpar b Mr Ang, Nov 07 Fin th valu of th constant k for which is a solution of th quation k. [7] Givn that, Givn that k, Thrfor, k Topic : Papr (00 marks) Tim : hours 0 mins Nam : Aitional

More information

A Bayesian criterion for simplicity in inverse problem parametrization

A Bayesian criterion for simplicity in inverse problem parametrization Gophys. J. Int. (000) 140, 67^85 A Baysian critrion for simplicity in invrs problm paramtrization Albrto Malinvrno Schlumbrgr-Doll Rsarch, Old Quarry Road, Ridg ld, CT 06877, USA. E-mail: albrto@ridg ld.sdr.slb.com

More information

Why is a E&M nature of light not sufficient to explain experiments?

Why is a E&M nature of light not sufficient to explain experiments? 1 Th wird world of photons Why is a E&M natur of light not sufficint to xplain xprimnts? Do photons xist? Som quantum proprtis of photons 2 Black body radiation Stfan s law: Enrgy/ ara/ tim = Win s displacmnt

More information

ME311 Machine Design

ME311 Machine Design ME311 Machin Dsign Lctur 4: Strss Concntrations; Static Failur W Dornfld 8Sp017 Fairfild Univrsit School of Enginring Strss Concntration W saw that in a curvd bam, th strss was distortd from th uniform

More information

Math 34A. Final Review

Math 34A. Final Review Math A Final Rviw 1) Us th graph of y10 to find approimat valus: a) 50 0. b) y (0.65) solution for part a) first writ an quation: 50 0. now tak th logarithm of both sids: log() log(50 0. ) pand th right

More information

INC 693, 481 Dynamics System and Modelling: Linear Graph Modeling II Dr.-Ing. Sudchai Boonto Assistant Professor

INC 693, 481 Dynamics System and Modelling: Linear Graph Modeling II Dr.-Ing. Sudchai Boonto Assistant Professor INC 69, 48 Dynamics Systm and Modlling: Linar Graph Modling II Dr.-Ing. Sudchai Boonto Assistant Profssor Dpartmnt of Control Systm and Instrumntation Enginring King Mongkut s Unnivrsity of Tchnology Thonuri

More information

Unit 30: Inference for Regression

Unit 30: Inference for Regression Unit 30: Infrnc for Rgrssion Summary of Vido In Unit 11, Fitting Lins to Data, w xamind th rlationship btwn wintr snowpack and spring runoff. Colorado rsourc managrs mad prdictions about th sasonal watr

More information

Lecture Outline. Skin Depth Power Flow 8/7/2018. EE 4347 Applied Electromagnetics. Topic 3e

Lecture Outline. Skin Depth Power Flow 8/7/2018. EE 4347 Applied Electromagnetics. Topic 3e 8/7/018 Cours Instructor Dr. Raymond C. Rumpf Offic: A 337 Phon: (915) 747 6958 E Mail: rcrumpf@utp.du EE 4347 Applid Elctromagntics Topic 3 Skin Dpth & Powr Flow Skin Dpth Ths & Powr nots Flow may contain

More information

Chapter 14 Aggregate Supply and the Short-run Tradeoff Between Inflation and Unemployment

Chapter 14 Aggregate Supply and the Short-run Tradeoff Between Inflation and Unemployment Chaptr 14 Aggrgat Supply and th Short-run Tradoff Btwn Inflation and Unmploymnt Modifid by Yun Wang Eco 3203 Intrmdiat Macroconomics Florida Intrnational Univrsity Summr 2017 2016 Worth Publishrs, all

More information

Finite element discretization of Laplace and Poisson equations

Finite element discretization of Laplace and Poisson equations Finit lmnt discrtization of Laplac and Poisson quations Yashwanth Tummala Tutor: Prof S.Mittal 1 Outlin Finit Elmnt Mthod for 1D Introduction to Poisson s and Laplac s Equations Finit Elmnt Mthod for 2D-Discrtization

More information

SECTION where P (cos θ, sin θ) and Q(cos θ, sin θ) are polynomials in cos θ and sin θ, provided Q is never equal to zero.

SECTION where P (cos θ, sin θ) and Q(cos θ, sin θ) are polynomials in cos θ and sin θ, provided Q is never equal to zero. SETION 6. 57 6. Evaluation of Dfinit Intgrals Exampl 6.6 W hav usd dfinit intgrals to valuat contour intgrals. It may com as a surpris to larn that contour intgrals and rsidus can b usd to valuat crtain

More information

Quasi-Classical States of the Simple Harmonic Oscillator

Quasi-Classical States of the Simple Harmonic Oscillator Quasi-Classical Stats of th Simpl Harmonic Oscillator (Draft Vrsion) Introduction: Why Look for Eignstats of th Annihilation Oprator? Excpt for th ground stat, th corrspondnc btwn th quantum nrgy ignstats

More information

Inflation and Unemployment

Inflation and Unemployment C H A P T E R 13 Aggrgat Supply and th Short-run Tradoff Btwn Inflation and Unmploymnt MACROECONOMICS SIXTH EDITION N. GREGORY MANKIW PowrPoint Slids by Ron Cronovich 2008 Worth Publishrs, all rights rsrvd

More information

MEMORIAL UNIVERSITY OF NEWFOUNDLAND

MEMORIAL UNIVERSITY OF NEWFOUNDLAND MEMORIAL UNIVERSITY OF NEWFOUNDLAND DEPARTMENT OF MATHEMATICS AND STATISTICS Midtrm Examination Statistics 2500 001 Wintr 2003 Nam: Studnt No: St by Dr. H. Wang OFFICE USE ONLY Mark: Instructions: 1. Plas

More information