36-720: Log-Linear Models: Three-Way Tables
|
|
- Madeline Hensley
- 6 years ago
- Views:
Transcription
1 36-720: Log-Linear Models: Three-Way Tables Brian Junker September 10, 2007 Hierarhial Speifiation of Log-Linear Models Higher-Dimensional Tables Digression: Why Poisson If We Believe Multinomial? Produt Multinomial Log-linear Models September 10, 2007
2 Hierarhial Speifiation of Log-Linear Models A full set of log-linear models for 3-way tables that we have onsidered so far inlude: Log-Linear Model Generator M ( 1) log m i jk = u [0] M (0) log m i jk = u+u 1(i) + u 2( j) + u 3(k) [1][2][3] M (1) log m i jk = u+u 1(i) + u 2( j) + u 3(k) + u 23( jk) [1][23] M (2) log m i jk = u+u 1(i) + u 2( j) + u 3(k) + u 13(ik) [2][13] M (3) log m i jk = u+u 1(i) + u 2( j) + u 3(k) + u 12(i j) [3][12] M (4) log m i jk = u+u 1(i) + u 2( j) + u 3(k) + u 13(ik) + u 23( jk) [13][23] M (5) log m i jk = u+u 1(i) + u 2( j) + u 3(k) + u 12(i j) + u 23( jk) [12][23] M (6) log m i jk = u+u 1(i) + u 2( j) + u 3(k) + u 12(i j) + u 13(ik) [12][13] M (7) log m i jk = u+u 1(i) + u 2( j) + u 3(k) + u 12(i j) + u 13(ik) + u 23( jk) [12][13][23] M (8) log m i jk = u+u 1(i) + u 2( j) + u 3(k) + u 12(i j) + u 13(ik) + u 23( jk) + u 123(i jk) [123] September 10, 2007
3 These model speifiations have several noteworthy features: They obey the hierarhy priniple: If the k-way interation is in the model then every lower order interation and main effet is also in the model. For example, sine M (7) ontains u 12(ik) terms then we also know it ontains u 1(i) and u 2( j) terms). Under the hierarhy priniple, nested models are obtained by adding or dropping higher-order interations in the model. Following the hierarhy priniple, a model an be ompletely speified by speifying the indies of the highest-order interations (the generators). E.g.: [1][23] u 1(i) and u 23( jk) are in the model u+u 1(i) and u+u 2( j) + u 3(k) + u 23( jk) are in the model log m i jk = u+u 1(i) + u 2( j) + u 3(k) + u 23( jk) is the model This is how R s linear modeling notation works also:shl + Risk*Beh expands to1 + Shl + Risk + Beh + Risk:Beh. The generators are useful mnemonis for What the (onditional) independene model is, e.g. [12][13] 2 3 1; What the suffiient statistis are, e.g. for [12][13] the suffiient statistis are n 12(i j) and n 13(ik) September 10, 2007
4 These ten hierarhial models an be organized by nesting as follows: [123] [12][13][23] [13][23] [12][23] [13][12] [1][23] [2][13] [3][12] [1][2][3] [0] September 10, 2007
5 Example: Residual Devianes (LR vs [123]) in boys deviane data Take 1=Beh, 2=Risk, 3=Shl: [123] (0/0) [12][13][23] (0.9/2) [13][23] (1.9/3) [12][23] (4.1/4) [13][12] (11.3/4) [1][23] (5.6/5) [2][13] (12.8/5) [3][12] (14.9/6) [1][2][3] (16.4/7) [0] (100.7/11) September 10, 2007
6 Higher-Dimensional Tables Christensen provides data on the relationship between two drugs (k=1, 2) and musle tensions (l=1, 2) for two weights (i=1, 2) and types ( j=1, 2) of musles in mie. Drug (k) Tension (l) Weight (i) Musle ( j) Drug 1 Drug 2 High Type Type High Low Type Type Low High Type Type Low Type Type Whih log-linear models best desribe this data? September 10, 2007
7 It is already painful to write down the saturated model [1234] in hierarhial log-linear form: log m i jk = u+u 1(i) + u 2( j) + u 3(k) + u 4(l) + u 12(i j) + u 13(ik) + u 14(il) + u 23( jk) + u 24( jl) + u 34( jl) + u 123(i jk) + u 124(i jl) + u 134(ikl) + u 234( jkl) + u 1234(i jkl) However, the generator notation makes model speifiation easier: The independene model [1][2][3][4] onsists of the first line above. The model of no three-way interation [12][13][14][23][24][34] onsists of the first two lines above. The model of no four-way interation [123][124][134][234] onsists of the first three lines above. The model [123][4] (or [WMD][T]) spefiies that musle weight, type and drug are independent of tension. The model [TMD][WMD] speifies that (tension) (weight) (type,drug) September 10, 2007
8 As an illustration we fit three models n <- san()... Wt <- rep(rep(("h","l"),(4,4)),2) # [W]eight Mt <- rep(rep(("1","2"),(2,2)),4) # [M]usle type Dt <- rep(("1","2"),8) # [D]rug type Tn <- rep(("h","l"),(8,8)) # [T]ension mus <- data.frame(n,wt,mt,dt,tn) Independene: [T][W][M][D] glm(n Tn + Wt + Mt + Dt, data = mus, family = poisson) No three-way interation: [TW][TM][TD][WM][WD][MD] glm(n Tn*Wt + Tn*Mt + Tn*Dt + Wt*Mt + Wt*Dt + Mt*Dt,...) No four-way interation: [TWM][TWD][TMD][WMD] glm(n Tn*Wt*Mt*Dt - Tn:Wt:Mt:Dt,...) Model Resid. Df Resid. Dev P[χ 2 d f> Dev] Independene No three-way e-09 No four-way September 10, 2007
9 Digression: Why Poisson If We Believe Multinomial? Maximum likelihood for the Poisson model We re-index a table suh as n i jk to be just n, and we suppose that n, =1,..., C are independent Poisson ounts with means m. Then f (n m ) = e m m n n! = exp[n log m m log n!] = exp[n θ expθ log n!] whih is an exponential family model with natural parameterθ = log m. If we modelθ = log m linearly as θ C 1 = [log m ] C 1 = X C D β D 1 then the above is a generalized linear model for Poisson ounts with the log link funtionθ = logµ September 10, 2007
10 The log-likelihood may be written L(n β) = f (n β) = exp (n log m m ) log n! = exp nt Xβ exp([xβ] )+g(data) and, either expliitly setting derivatives equal to zero, or by using the general theory of glm s, we see that the likelihood equations redue to (n m )x d = 0, d= 1,..., D This is the usual result that the MLE in an exponential family model equates observed and expeted suffiient statistis n T X= m T X In log-linear models for tables, these are invariably appropriate marginal totals for the table (see example, next slide) September 10, 2007
11 Sub-digression: Where does X ome from? Consider the table n i jk with i=1, 2, j=1, 2, and k=1, 2. If we lay out the log expeted ell ounts in a olumn, the model of independene log m i jk = u+u 1(i) + u 2( j) + u 3(k) looks like this: θ= log m 111 log m 112 log m 121 log m 122 log m 211 log m 212 log m 221 log m 222 = u u 1(2) u 2(2) u 3(2) = Xβ Clearly the observed suffiient statistis for the u parameters are n T X= (n +++, n 2++, n +2+, n ++2 ) whih are equated by ML to the expeted suffiient statistis m T X= (m +++, m 2++, m +2+, m ++2 ) September 10, 2007
12 Bak to the Poisson Likelihood Further differentiation (or appliation of general results from glm s) shows that the information matrix forβis [ 2 ] = X T ŴX β d1 β d2 where W= diag( ˆm) and hene D D Var (β)=[x T ŴX] 1 Sine the MLE s under the saturated model (no relations among the m i s) are ˆm i = n i, the log-lr statisti for testing against the saturated model is [ ] 2[log L(n β) log L(n saturated)] = 2 n log(n / ˆm ) n + ˆm = 2 n log(n / ˆm ) sine n + = ˆm + as long as the log-linear model has the interept u in it September 10, 2007
13 Sine the individual ells ounts n i Poiss(m i ), with E[n i ]=Var (n i )=m i, it follows that the Pearson residuals r i = (n i ˆm i )/ ˆm i are approximately mean 0, variane 1 (and should be approximately normal for large n i ). This is why they are sensible residuals to use. Reall that the model deviane (LR statisti for testing against the saturated model) is G 2 = 2 n i log(n i / ˆm i )= The deviane residuals are i r D i = sgn (n i ˆm i ) d i and these are approximately normally distributed (sine they are signed square roots of approximately 1-df ontributions to G 2 ). i d i September 10, 2007
14 Maximum likelihood for the Multinomial model We again take a table like n i jk and re-index the ells as n 1,...,n C. The key observation here is that if we ondition on n + = n, C independent Poisson s beome a C-ell multinomial. Begin with a Poisson log-linear model with an intereptα, log m =α+ x T β (=u+[other u-terms] ) where (1, x T ) is the th row of X. Up to a funtion of the data g(data)= log n!, the Poisson log-likelihood is log L(n α, β) = n log m = m n (α+ xβ) T exp(α+ xβ) T = n + α+ n xβ τ T whereτ= m = exp(α+ x T β)=e α exp(x T β) September 10, 2007
15 Now sine logτ=α+log exp(x T β), it follows that log L(n α, β) = log L(n τ, β) = n xβ n T log exp(xβ) T + [n + logτ τ] = n xβ T n log exp(x T β) + [n + logτ τ] = n log p + [n + logτ τ] where p = m /m + = exp(α+ x T β)/ exp(α+ xt β)= exp(x T β)/ exp(xt β). The term in brakets is the multinomial log-likelihood, up to a funtion of the data h(data)=log n +! log n!. Thus, as log-likelihoods: Poisson α,β (n)={multinomial β (n) n + }+[Poisson τ (n + ) ] September 10, 2007
16 That is worth repeating: For the log-linear model log m =α+ x T β, Poisson α,β (n)={multinomial β (n) n + }+[Poisson τ (n + ) ]. as log-likelihoods. Therefore: There is a 1-1 orrespondene between Multinomial log-linear models and Poisson log-linear models with the same log-linear form. The Poisson model requires one more parameter (τ), orresponding to the grand total n +. The (non-interept) parametersβin the two models have the same MLE s and the same variane-ovariane matrix (sine theβ s only enter in the braketed expressions on the previous slide). The interept parameterαwill be different in the two models. The Pearson residuals are still variane-stabilized, and as long as the grand total n + is large relative to the ell ount n i, eah r i will be roughly N(0, 1). You an also show that the omputation and interpretation of G 2 and X 2 also do not hange between the models September 10, 2007
17 Produt Multinomial Log-linear Models A similar argument an be used to derive produt-multinomial MLE s from Poisson MLE s. The key to the equivalene between Poisson MLE s ˆβ and single-multinomial MLE s ˆβ was to inlude the intereptαorresponding to the fixed grand total n + in the log-linear model: log m =α+ x T β. Then as log-likelihoods Poisson α,β (n)={multinomial β (n) n + }+[Poisson τ (n + ) ]. The key to equivalene between Poisson MLE s and Produt Multinomial MLE s will be to inlude log-linear termsα (1),...,α (H) orresponding to the fixed margins n (1) +,..., n (H) +. Then as log-likelihoods Poisson α,β (n)= H h=1 { Multinomial β (n (h) ) n (h) + }+[Poisson τ (n (h) + ) ]. These results have been known for some time (e.g. Birh, 1963, JRSSB); a reent update/generalization is Lang (1996, JRSSB) September 10, 2007
18 Example Consider again the Aspirin and Heart Attak data: Myoardial Infartion Fatal Nonfatal Attak Attak Total Plaebo n 11 = 18 n 12 = 171 n 1+ = 189 Aspirin n 21 = 5 n 22 = 99 n 2+ = 104 Total n +1 = 23 n +2 = 270 n ++ = 293 If we onsider this to be observational data with only n ++ fixed and a multinomial model for n i j n ++, we an useglm(..., family=poisson) to fit log-linear models to this data as long as we inlude the interept u in eah log-linear model. If we onsider this to be a designed experiment so that the ells are produt multinomial with fixed row totals n 1+ and n 2+, then we an useglm(..., family=poisson) to fit log-linear models to this data as long as we inlude the u 1(i) terms in eah log-linear model September 10, 2007
19 Testing independene with the Single Multinomial Model: n <- san() Tx <- rep(("plaebo","aspirin"),(3,3)) Obs <- rep(("fatal","nonfatal","noattak"),2) aha.data <- data.frame(n,tx,obs) print(fit <- glm(n Tx + Obs,data=aha.data,family=poisson)) G 2 = on 2 d.f.; p , so again we rejet independene. Testing independene with the Produt Multinomial Model: Sine the model of independene already has log-linear terms (Tx, i.e. u 1(i) ) orresponding to the fixed row totals n i+, the same Poisson fit above also gives the results for the Produt Multinomial. In general, when H 0 already ontains terms orresponding to the fixed margins in the table, testing H 0 is idential under the Poisson, Multinomial, and Produt Multinomial sampling models September 10, 2007
20 Example Bak to the musle study in mie: Drug (k) Tension (l) Weight (i) Musle ( j) Drug 1 Drug 2 High Type Type High Low Type Type Low High Type Type Low Type Type We treated this data before as Poisson or single multinomial. In fat, it was a designed study with the total number of musles of eah type fixed in advane. So every log-linear model should have themterms (u 2( j) ) in it September 10, 2007
21 An advantage of hierarhial priniple is that the interesting models automatially ontain the main effets, so if the totals in one dimension are fixed, the sampling sheme doesn t matter for model fit/omparison: Model Resid. Df Resid. Dev P[χ 2 d f> Dev] [T][W][M][D] [TW][TM][TD][WM][WD][MD] e-09 [TWM][TWD][TMD][WMD] In the study, [T]ension and [W]eight of musle are measured on eah ombination of [M]usle type and [D]rug. Two more models of interest might be 1. [TMD][WMD]: (tension) (weight) (musle type,drug) 2. [TWM][DM]: (tension,weight) (drug) (musle type) In the real study, only the totals for musle [T]ype were fixed, so either of these models ould be fitted as well September 10, 2007
22 In another study of this type, perhaps musle [W]eight and [D]rug type would be fixed in advane. In that ase all models would need to inlude the [WD] interation: we ould fit and interpret [TMD][WMD] but not [TWM][DM] (without the [WD] interation, the totals for (weight) (drug) ombinations would not be fixed by the log-linear model; the produt multinomial model with [WD] margins fixed ould not be represented) September 10, 2007
Methods of evaluating tests
Methods of evaluating tests Let X,, 1 Xn be i.i.d. Bernoulli( p ). Then 5 j= 1 j ( 5, ) T = X Binomial p. We test 1 H : p vs. 1 1 H : p>. We saw that a LRT is 1 if t k* φ ( x ) =. otherwise (t is the observed
More informationAdvanced Computational Fluid Dynamics AA215A Lecture 4
Advaned Computational Fluid Dynamis AA5A Leture 4 Antony Jameson Winter Quarter,, Stanford, CA Abstrat Leture 4 overs analysis of the equations of gas dynamis Contents Analysis of the equations of gas
More information11. Categorical choice and survival models
11. Categorial hoie and survival models 11.1 Homogeneous models 11.2 Multinomial log models wh random effets 11.3 Transion (Markov) models 11.4 Survival models 11.1 Homogeneous models 11.1.1 Statistial
More informationON A PROCESS DERIVED FROM A FILTERED POISSON PROCESS
ON A PROCESS DERIVED FROM A FILTERED POISSON PROCESS MARIO LEFEBVRE and JEAN-LUC GUILBAULT A ontinuous-time and ontinuous-state stohasti proess, denoted by {Xt), t }, is defined from a proess known as
More informationChapter 8 Hypothesis Testing
Leture 5 for BST 63: Statistial Theory II Kui Zhang, Spring Chapter 8 Hypothesis Testing Setion 8 Introdution Definition 8 A hypothesis is a statement about a population parameter Definition 8 The two
More informationMaximum Entropy and Exponential Families
Maximum Entropy and Exponential Families April 9, 209 Abstrat The goal of this note is to derive the exponential form of probability distribution from more basi onsiderations, in partiular Entropy. It
More informationTENSOR FORM OF SPECIAL RELATIVITY
TENSOR FORM OF SPECIAL RELATIVITY We begin by realling that the fundamental priniple of Speial Relativity is that all physial laws must look the same to all inertial observers. This is easiest done by
More informationComputer Science 786S - Statistical Methods in Natural Language Processing and Data Analysis Page 1
Computer Siene 786S - Statistial Methods in Natural Language Proessing and Data Analysis Page 1 Hypothesis Testing A statistial hypothesis is a statement about the nature of the distribution of a random
More informationProbabilistic Graphical Models
Probabilisti Graphial Models David Sontag New York University Leture 12, April 19, 2012 Aknowledgement: Partially based on slides by Eri Xing at CMU and Andrew MCallum at UMass Amherst David Sontag (NYU)
More informationQuantum Mechanics: Wheeler: Physics 6210
Quantum Mehanis: Wheeler: Physis 60 Problems some modified from Sakurai, hapter. W. S..: The Pauli matries, σ i, are a triple of matries, σ, σ i = σ, σ, σ 3 given by σ = σ = σ 3 = i i Let stand for the
More informationTests of fit for symmetric variance gamma distributions
Tests of fit for symmetri variane gamma distributions Fragiadakis Kostas UADPhilEon, National and Kapodistrian University of Athens, 4 Euripidou Street, 05 59 Athens, Greee. Keywords: Variane Gamma Distribution,
More informationSupplementary Materials
Supplementary Materials Neural population partitioning and a onurrent brain-mahine interfae for sequential motor funtion Maryam M. Shanehi, Rollin C. Hu, Marissa Powers, Gregory W. Wornell, Emery N. Brown
More informationHankel Optimal Model Order Reduction 1
Massahusetts Institute of Tehnology Department of Eletrial Engineering and Computer Siene 6.245: MULTIVARIABLE CONTROL SYSTEMS by A. Megretski Hankel Optimal Model Order Redution 1 This leture overs both
More information36-720: Log-Linear Models: Three-Way Tables
36-720: Log-Linear Models: Three-Way Tables Brian Junker September 5, 2007 Aggregation and Association Modelling the Three-Way Table The Model of Complete Independence Models with One Factor Independent
More informationEvaluation of effect of blade internal modes on sensitivity of Advanced LIGO
Evaluation of effet of blade internal modes on sensitivity of Advaned LIGO T0074-00-R Norna A Robertson 5 th Otober 00. Introdution The urrent model used to estimate the isolation ahieved by the quadruple
More informationSOA/CAS MAY 2003 COURSE 1 EXAM SOLUTIONS
SOA/CAS MAY 2003 COURSE 1 EXAM SOLUTIONS Prepared by S. Broverman e-mail 2brove@rogers.om website http://members.rogers.om/2brove 1. We identify the following events:. - wathed gymnastis, ) - wathed baseball,
More informationDanielle Maddix AA238 Final Project December 9, 2016
Struture and Parameter Learning in Bayesian Networks with Appliations to Prediting Breast Caner Tumor Malignany in a Lower Dimension Feature Spae Danielle Maddix AA238 Final Projet Deember 9, 2016 Abstrat
More informationarxiv:math/ v1 [math.ca] 27 Nov 2003
arxiv:math/011510v1 [math.ca] 27 Nov 200 Counting Integral Lamé Equations by Means of Dessins d Enfants Sander Dahmen November 27, 200 Abstrat We obtain an expliit formula for the number of Lamé equations
More informationWhere as discussed previously we interpret solutions to this partial differential equation in the weak sense: b
Consider the pure initial value problem for a homogeneous system of onservation laws with no soure terms in one spae dimension: Where as disussed previously we interpret solutions to this partial differential
More informationSupporting Information
Supporting Information Olsman and Goentoro 10.1073/pnas.1601791113 SI Materials Analysis of the Sensitivity and Error Funtions. We now define the sensitivity funtion Sð, «0 Þ, whih summarizes the steepness
More information7 Max-Flow Problems. Business Computing and Operations Research 608
7 Max-Flow Problems Business Computing and Operations Researh 68 7. Max-Flow Problems In what follows, we onsider a somewhat modified problem onstellation Instead of osts of transmission, vetor now indiates
More informationGeneral Closed-form Analytical Expressions of Air-gap Inductances for Surfacemounted Permanent Magnet and Induction Machines
General Closed-form Analytial Expressions of Air-gap Indutanes for Surfaemounted Permanent Magnet and Indution Mahines Ronghai Qu, Member, IEEE Eletroni & Photoni Systems Tehnologies General Eletri Company
More informationPerforming Two-Way Analysis of Variance Under Variance Heterogeneity
Journal of Modern Applied Statistial Methods Volume Issue Artile 3 5--003 Performing Two-Way Analysis of Variane Under Variane Heterogeneity Sott J. Rihter University of North Carolina at Greensboro, sjriht@ung.edu
More informationA Characterization of Wavelet Convergence in Sobolev Spaces
A Charaterization of Wavelet Convergene in Sobolev Spaes Mark A. Kon 1 oston University Louise Arakelian Raphael Howard University Dediated to Prof. Robert Carroll on the oasion of his 70th birthday. Abstrat
More informationEstimation of Density Function Parameters with Censored Data from Product Life Tests
VŠB TU Ostrava Faulty of Mehanial Engineering Department of Control Systems and Instrumentation Estimation of Density Funtion Parameters with Censored Data from Produt Life Tests JiříTůma 7. listopadu
More informationEstimation for Unknown Parameters of the Extended Burr Type-XII Distribution Based on Type-I Hybrid Progressive Censoring Scheme
Estimation for Unknown Parameters of the Extended Burr Type-XII Distribution Based on Type-I Hybrid Progressive Censoring Sheme M.M. Amein Department of Mathematis Statistis, Faulty of Siene, Taif University,
More informationDeveloping Excel Macros for Solving Heat Diffusion Problems
Session 50 Developing Exel Maros for Solving Heat Diffusion Problems N. N. Sarker and M. A. Ketkar Department of Engineering Tehnology Prairie View A&M University Prairie View, TX 77446 Abstrat This paper
More informationLine Radiative Transfer
http://www.v.nrao.edu/ourse/astr534/ineradxfer.html ine Radiative Transfer Einstein Coeffiients We used armor's equation to estimate the spontaneous emission oeffiients A U for À reombination lines. A
More informationComplexity of Regularization RBF Networks
Complexity of Regularization RBF Networks Mark A Kon Department of Mathematis and Statistis Boston University Boston, MA 02215 mkon@buedu Leszek Plaskota Institute of Applied Mathematis University of Warsaw
More informationTaste for variety and optimum product diversity in an open economy
Taste for variety and optimum produt diversity in an open eonomy Javier Coto-Martínez City University Paul Levine University of Surrey Otober 0, 005 María D.C. Garía-Alonso University of Kent Abstrat We
More informationA variant of Coppersmith s Algorithm with Improved Complexity and Efficient Exhaustive Search
A variant of Coppersmith s Algorithm with Improved Complexity and Effiient Exhaustive Searh Jean-Sébastien Coron 1, Jean-Charles Faugère 2, Guénaël Renault 2, and Rina Zeitoun 2,3 1 University of Luxembourg
More informationAdvances in Radio Science
Advanes in adio Siene 2003) 1: 99 104 Copernius GmbH 2003 Advanes in adio Siene A hybrid method ombining the FDTD and a time domain boundary-integral equation marhing-on-in-time algorithm A Beker and V
More informationRemark 4.1 Unlike Lyapunov theorems, LaSalle s theorem does not require the function V ( x ) to be positive definite.
Leture Remark 4.1 Unlike Lyapunov theorems, LaSalle s theorem does not require the funtion V ( x ) to be positive definite. ost often, our interest will be to show that x( t) as t. For that we will need
More informationRobust Flight Control Design for a Turn Coordination System with Parameter Uncertainties
Amerian Journal of Applied Sienes 4 (7): 496-501, 007 ISSN 1546-939 007 Siene Publiations Robust Flight ontrol Design for a urn oordination System with Parameter Unertainties 1 Ari Legowo and Hiroshi Okubo
More informationRobust Recovery of Signals From a Structured Union of Subspaces
Robust Reovery of Signals From a Strutured Union of Subspaes 1 Yonina C. Eldar, Senior Member, IEEE and Moshe Mishali, Student Member, IEEE arxiv:87.4581v2 [nlin.cg] 3 Mar 29 Abstrat Traditional sampling
More informationn n=1 (air) n 1 sin 2 r =
Physis 55 Fall 7 Homework Assignment #11 Solutions Textbook problems: Ch. 7: 7.3, 7.4, 7.6, 7.8 7.3 Two plane semi-infinite slabs of the same uniform, isotropi, nonpermeable, lossless dieletri with index
More informationChapter Review of of Random Processes
Chapter.. Review of of Random Proesses Random Variables and Error Funtions Conepts of Random Proesses 3 Wide-sense Stationary Proesses and Transmission over LTI 4 White Gaussian Noise Proesses @G.Gong
More informationCounting Idempotent Relations
Counting Idempotent Relations Beriht-Nr. 2008-15 Florian Kammüller ISSN 1436-9915 2 Abstrat This artile introdues and motivates idempotent relations. It summarizes haraterizations of idempotents and their
More informationCONDITIONAL CONFIDENCE INTERVAL FOR THE SCALE PARAMETER OF A WEIBULL DISTRIBUTION. Smail Mahdi
Serdia Math. J. 30 (2004), 55 70 CONDITIONAL CONFIDENCE INTERVAL FOR THE SCALE PARAMETER OF A WEIBULL DISTRIBUTION Smail Mahdi Communiated by N. M. Yanev Abstrat. A two-sided onditional onfidene interval
More informationV. Interacting Particles
V. Interating Partiles V.A The Cumulant Expansion The examples studied in the previous setion involve non-interating partiles. It is preisely the lak of interations that renders these problems exatly solvable.
More information8.333: Statistical Mechanics I Problem Set # 4 Due: 11/13/13 Non-interacting particles
8.333: Statistial Mehanis I Problem Set # 4 Due: 11/13/13 Non-interating partiles 1. Rotating gas: Consider a gas of N idential atoms onfined to a spherial harmoni trap in three dimensions, i.e. the partiles
More informationLikelihood-confidence intervals for quantiles in Extreme Value Distributions
Likelihood-onfidene intervals for quantiles in Extreme Value Distributions A. Bolívar, E. Díaz-Franés, J. Ortega, and E. Vilhis. Centro de Investigaión en Matemátias; A.P. 42, Guanajuato, Gto. 36; Méxio
More information1 sin 2 r = 1 n 2 sin 2 i
Physis 505 Fall 005 Homework Assignment #11 Solutions Textbook problems: Ch. 7: 7.3, 7.5, 7.8, 7.16 7.3 Two plane semi-infinite slabs of the same uniform, isotropi, nonpermeable, lossless dieletri with
More informationSubject: Introduction to Component Matching and Off-Design Operation % % ( (1) R T % (
16.50 Leture 0 Subjet: Introdution to Component Mathing and Off-Design Operation At this point it is well to reflet on whih of the many parameters we have introdued (like M, τ, τ t, ϑ t, f, et.) are free
More informationMathacle. PSet Stats, Concepts In Statistics Level Number Name: Date: χ = npq
8.6. Chi-Square ( χ ) Test [MATH] From the DeMoivre -- Laplae Theorem, when ( 1 p) >> 1 1 m np b( n, p, m) ϕ npq npq np, where m is the observed number of suesses in n trials, the probability of suess
More informationST3241 Categorical Data Analysis I Generalized Linear Models. Introduction and Some Examples
ST3241 Categorical Data Analysis I Generalized Linear Models Introduction and Some Examples 1 Introduction We have discussed methods for analyzing associations in two-way and three-way tables. Now we will
More informationLecture 7: Sampling/Projections for Least-squares Approximation, Cont. 7 Sampling/Projections for Least-squares Approximation, Cont.
Stat60/CS94: Randomized Algorithms for Matries and Data Leture 7-09/5/013 Leture 7: Sampling/Projetions for Least-squares Approximation, Cont. Leturer: Mihael Mahoney Sribe: Mihael Mahoney Warning: these
More informationDiscrete Bessel functions and partial difference equations
Disrete Bessel funtions and partial differene equations Antonín Slavík Charles University, Faulty of Mathematis and Physis, Sokolovská 83, 186 75 Praha 8, Czeh Republi E-mail: slavik@karlin.mff.uni.z Abstrat
More information5.1 Composite Functions
SECTION. Composite Funtions 7. Composite Funtions PREPARING FOR THIS SECTION Before getting started, review the following: Find the Value of a Funtion (Setion., pp. 9 ) Domain of a Funtion (Setion., pp.
More information(a) We desribe physics as a sequence of events labelled by their space time coordinates: x µ = (x 0, x 1, x 2 x 3 ) = (c t, x) (12.
2 Relativity Postulates (a) All inertial observers have the same equations of motion and the same physial laws. Relativity explains how to translate the measurements and events aording to one inertial
More informationNonreversibility of Multiple Unicast Networks
Nonreversibility of Multiple Uniast Networks Randall Dougherty and Kenneth Zeger September 27, 2005 Abstrat We prove that for any finite direted ayli network, there exists a orresponding multiple uniast
More informationCS 687 Jana Kosecka. Uncertainty, Bayesian Networks Chapter 13, Russell and Norvig Chapter 14,
CS 687 Jana Koseka Unertainty Bayesian Networks Chapter 13 Russell and Norvig Chapter 14 14.1-14.3 Outline Unertainty robability Syntax and Semantis Inferene Independene and Bayes' Rule Syntax Basi element:
More informationThe Unified Geometrical Theory of Fields and Particles
Applied Mathematis, 014, 5, 347-351 Published Online February 014 (http://www.sirp.org/journal/am) http://dx.doi.org/10.436/am.014.53036 The Unified Geometrial Theory of Fields and Partiles Amagh Nduka
More informationGeneration of EM waves
Generation of EM waves Susan Lea Spring 015 1 The Green s funtion In Lorentz gauge, we obtained the wave equation: A 4π J 1 The orresponding Green s funtion for the problem satisfies the simpler differential
More informationThe Power of the Test in Three-Level Designs
DISCUSSION PAPER SERIES IZA DP No. 4 The Power of the Test in Three-Level Designs Spyros Konstantopoulos Otober 006 Forshungsinstitut zur Zukunft der Arbeit Institute for the Study of Labor The Power of
More informationThe Effectiveness of the Linear Hull Effect
The Effetiveness of the Linear Hull Effet S. Murphy Tehnial Report RHUL MA 009 9 6 Otober 009 Department of Mathematis Royal Holloway, University of London Egham, Surrey TW0 0EX, England http://www.rhul.a.uk/mathematis/tehreports
More informationHandling Uncertainty
Handling Unertainty Unertain knowledge Typial example: Diagnosis. Name Toothahe Cavity Can we ertainly derive the diagnosti rule: if Toothahe=true then Cavity=true? The problem is that this rule isn t
More information2. The Energy Principle in Open Channel Flows
. The Energy Priniple in Open Channel Flows. Basi Energy Equation In the one-dimensional analysis of steady open-hannel flow, the energy equation in the form of Bernoulli equation is used. Aording to this
More informationIMPEDANCE EFFECTS OF LEFT TURNERS FROM THE MAJOR STREET AT A TWSC INTERSECTION
09-1289 Citation: Brilon, W. (2009): Impedane Effets of Left Turners from the Major Street at A TWSC Intersetion. Transportation Researh Reord Nr. 2130, pp. 2-8 IMPEDANCE EFFECTS OF LEFT TURNERS FROM THE
More informationWave Propagation through Random Media
Chapter 3. Wave Propagation through Random Media 3. Charateristis of Wave Behavior Sound propagation through random media is the entral part of this investigation. This hapter presents a frame of referene
More information2 The Bayesian Perspective of Distributions Viewed as Information
A PRIMER ON BAYESIAN INFERENCE For the next few assignments, we are going to fous on the Bayesian way of thinking and learn how a Bayesian approahes the problem of statistial modeling and inferene. The
More informationEECS 120 Signals & Systems University of California, Berkeley: Fall 2005 Gastpar November 16, Solutions to Exam 2
EECS 0 Signals & Systems University of California, Berkeley: Fall 005 Gastpar November 6, 005 Solutions to Exam Last name First name SID You have hour and 45 minutes to omplete this exam. he exam is losed-book
More informationDesign and Development of Three Stages Mixed Sampling Plans for Variable Attribute Variable Quality Characteristics
International Journal of Statistis and Systems ISSN 0973-2675 Volume 12, Number 4 (2017), pp. 763-772 Researh India Publiations http://www.ripubliation.om Design and Development of Three Stages Mixed Sampling
More informationmax min z i i=1 x j k s.t. j=1 x j j:i T j
AM 221: Advaned Optimization Spring 2016 Prof. Yaron Singer Leture 22 April 18th 1 Overview In this leture, we will study the pipage rounding tehnique whih is a deterministi rounding proedure that an be
More informationGeometry of Transformations of Random Variables
Geometry of Transformations of Random Variables Univariate distributions We are interested in the problem of finding the distribution of Y = h(x) when the transformation h is one-to-one so that there is
More informationStructural Reconfiguration of Systems under Behavioral Adaptation
Strutural Reonfiguration of Systems under Behavioral Adaptation Carlos Canal a, Javier Cámara b, Gwen Salaün a Department of Computer Siene, University of Málaga, Spain b Department of Informatis Engineering,
More informationMAC Calculus II Summer All you need to know on partial fractions and more
MC -75-Calulus II Summer 00 ll you need to know on partial frations and more What are partial frations? following forms:.... where, α are onstants. Partial frations are frations of one of the + α, ( +
More informationGeneral Equilibrium. What happens to cause a reaction to come to equilibrium?
General Equilibrium Chemial Equilibrium Most hemial reations that are enountered are reversible. In other words, they go fairly easily in either the forward or reverse diretions. The thing to remember
More informationLecturer: Bengt E W Nilsson
9 4 3 Leturer: Bengt E W Nilsson 8::: Leturer absent, three students present. 8::9: Leturer present, three students. 8::: Six students. 8::8: Five students. Waiting a ouple of minutes to see if more ome.
More informationSome recent developments in probability distributions
Proeedings 59th ISI World Statistis Congress, 25-30 August 2013, Hong Kong (Session STS084) p.2873 Some reent developments in probability distributions Felix Famoye *1, Carl Lee 1, and Ayman Alzaatreh
More informationA model for measurement of the states in a coupled-dot qubit
A model for measurement of the states in a oupled-dot qubit H B Sun and H M Wiseman Centre for Quantum Computer Tehnology Centre for Quantum Dynamis Griffith University Brisbane 4 QLD Australia E-mail:
More informationMost results in this section are stated without proof.
Leture 8 Level 4 v2 he Expliit formula. Most results in this setion are stated without proof. Reall that we have shown that ζ (s has only one pole, a simple one at s =. It has trivial zeros at the negative
More informationExplicit formulae for integrals of s and p type GTFs
ELECTRONIC JOURNAL OF THEORETICAL CHEMISTRY, VOL., 66 70 (997) Expliit formulae for integrals of s and p type GTFs CAROL A. BAXTER AND DAVID B. COOK The Department of Chemistry, The University of Sheffield,
More informationMillennium Relativity Acceleration Composition. The Relativistic Relationship between Acceleration and Uniform Motion
Millennium Relativity Aeleration Composition he Relativisti Relationship between Aeleration and niform Motion Copyright 003 Joseph A. Rybzyk Abstrat he relativisti priniples developed throughout the six
More informationarxiv:gr-qc/ v2 6 Feb 2004
Hubble Red Shift and the Anomalous Aeleration of Pioneer 0 and arxiv:gr-q/0402024v2 6 Feb 2004 Kostadin Trenčevski Faulty of Natural Sienes and Mathematis, P.O.Box 62, 000 Skopje, Maedonia Abstrat It this
More informationQ2. [40 points] Bishop-Hill Model: Calculation of Taylor Factors for Multiple Slip
27-750, A.D. Rollett Due: 20 th Ot., 2011. Homework 5, Volume Frations, Single and Multiple Slip Crystal Plastiity Note the 2 extra redit questions (at the end). Q1. [40 points] Single Slip: Calulating
More informationeappendix for: SAS macro for causal mediation analysis with survival data
eappendix for: SAS maro for ausal mediation analysis with survival data Linda Valeri and Tyler J. VanderWeele 1 Causal effets under the ounterfatual framework and their estimators We let T a and M a denote
More informationMOLECULAR ORBITAL THEORY- PART I
5.6 Physial Chemistry Leture #24-25 MOLECULAR ORBITAL THEORY- PART I At this point, we have nearly ompleted our rash-ourse introdution to quantum mehanis and we re finally ready to deal with moleules.
More informationIntroduction to General and Generalized Linear Models
Introduction to General and Generalized Linear Models Generalized Linear Models - part II Henrik Madsen Poul Thyregod Informatics and Mathematical Modelling Technical University of Denmark DK-2800 Kgs.
More informationA simple expression for radial distribution functions of pure fluids and mixtures
A simple expression for radial distribution funtions of pure fluids and mixtures Enrio Matteoli a) Istituto di Chimia Quantistia ed Energetia Moleolare, CNR, Via Risorgimento, 35, 56126 Pisa, Italy G.
More informationAharonov-Bohm effect. Dan Solomon.
Aharonov-Bohm effet. Dan Solomon. In the figure the magneti field is onfined to a solenoid of radius r 0 and is direted in the z- diretion, out of the paper. The solenoid is surrounded by a barrier that
More informationSPLINE ESTIMATION OF SINGLE-INDEX MODELS
SPLINE ESIMAION OF SINGLE-INDEX MODELS Li Wang and Lijian Yang University of Georgia and Mihigan State University Supplementary Material his note ontains proofs for the main results he following two propositions
More informationCovariances for Bivariate Selection Model Second-Step Estimates
Eon 240B, Fall 2001 Daniel MFadden Covarianes for Bivariate Seletion Model Seond-Step Estimates Consider the bivariate latent variable model with normal disturbanes, (1) y * = + C, w * = zα + σν, where
More information22.54 Neutron Interactions and Applications (Spring 2004) Chapter 6 (2/24/04) Energy Transfer Kernel F(E E')
22.54 Neutron Interations and Appliations (Spring 2004) Chapter 6 (2/24/04) Energy Transfer Kernel F(E E') Referenes -- J. R. Lamarsh, Introdution to Nulear Reator Theory (Addison-Wesley, Reading, 1966),
More informationThe Corpuscular Structure of Matter, the Interaction of Material Particles, and Quantum Phenomena as a Consequence of Selfvariations.
The Corpusular Struture of Matter, the Interation of Material Partiles, and Quantum Phenomena as a Consequene of Selfvariations. Emmanuil Manousos APM Institute for the Advanement of Physis and Mathematis,
More informationNormative and descriptive approaches to multiattribute decision making
De. 009, Volume 8, No. (Serial No.78) China-USA Business Review, ISSN 57-54, USA Normative and desriptive approahes to multiattribute deision making Milan Terek (Department of Statistis, University of
More informationCase I: 2 users In case of 2 users, the probability of error for user 1 was earlier derived to be 2 A1
MUTLIUSER DETECTION (Letures 9 and 0) 6:33:546 Wireless Communiations Tehnologies Instrutor: Dr. Narayan Mandayam Summary By Shweta Shrivastava (shwetash@winlab.rutgers.edu) bstrat This artile ontinues
More informationAn iterative least-square method suitable for solving large sparse matrices
An iteratie least-square method suitable for soling large sparse matries By I. M. Khabaza The purpose of this paper is to report on the results of numerial experiments with an iteratie least-square method
More informationOn Component Order Edge Reliability and the Existence of Uniformly Most Reliable Unicycles
Daniel Gross, Lakshmi Iswara, L. William Kazmierzak, Kristi Luttrell, John T. Saoman, Charles Suffel On Component Order Edge Reliability and the Existene of Uniformly Most Reliable Uniyles DANIEL GROSS
More informationResearch Letter Distributed Source Localization Based on TOA Measurements in Wireless Sensor Networks
Researh Letters in Eletronis Volume 2009, Artile ID 573129, 4 pages doi:10.1155/2009/573129 Researh Letter Distributed Soure Loalization Based on TOA Measurements in Wireless Sensor Networks Wanzhi Qiu
More informationSupplementary Information. Infrared Transparent Visible Opaque Fabrics (ITVOF) for Personal Cooling
Supplementary Information Infrared Transparent Visible Opaque Fabris (ITVOF) for Personal Cooling Jonathan K. Tong 1,Ɨ, Xiaopeng Huang 1,Ɨ, Svetlana V. Boriskina 1, James Loomis 1, Yanfei Xu 1, and Gang
More informationCapacity-achieving Input Covariance for Correlated Multi-Antenna Channels
Capaity-ahieving Input Covariane for Correlated Multi-Antenna Channels Antonia M. Tulino Universita Degli Studi di Napoli Federio II 85 Napoli, Italy atulino@ee.prineton.edu Angel Lozano Bell Labs (Luent
More informationA Spatiotemporal Approach to Passive Sound Source Localization
A Spatiotemporal Approah Passive Sound Soure Loalization Pasi Pertilä, Mikko Parviainen, Teemu Korhonen and Ari Visa Institute of Signal Proessing Tampere University of Tehnology, P.O.Box 553, FIN-330,
More informationFour-dimensional equation of motion for viscous compressible substance with regard to the acceleration field, pressure field and dissipation field
Four-dimensional equation of motion for visous ompressible substane with regard to the aeleration field, pressure field and dissipation field Sergey G. Fedosin PO box 6488, Sviazeva str. -79, Perm, Russia
More informationGreen s function for the wave equation
Green s funtion for the wave equation Non-relativisti ase January 2019 1 The wave equations In the Lorentz Gauge, the wave equations for the potentials are (Notes 1 eqns 43 and 44): 1 2 A 2 2 2 A = µ 0
More informationTheory. Coupled Rooms
Theory of Coupled Rooms For: nternal only Report No.: R/50/TCR Prepared by:. N. taey B.., MO Otober 00 .00 Objet.. The objet of this doument is present the theory alulations to estimate the reverberant
More informationDetermination of the reaction order
5/7/07 A quote of the wee (or amel of the wee): Apply yourself. Get all the eduation you an, but then... do something. Don't just stand there, mae it happen. Lee Iaoa Physial Chemistry GTM/5 reation order
More informationDelivery guide for Environmental Noise Data:
Delivery guide for Environmental Noise Data: DF4_DF8: Strategi noise maps for major roads, major railways, major airports and agglomerations Type of Doument: Guidelines Annex DF4_DF8 Prepared by: Colin
More informationConcerning the Numbers 22p + 1, p Prime
Conerning the Numbers 22p + 1, p Prime By John Brillhart 1. Introdution. In a reent investigation [7] the problem of fatoring numbers of the form 22p + 1, p a, was enountered. Sine 22p + 1 = (2P - 2*
More informationStrauss PDEs 2e: Section Exercise 3 Page 1 of 13. u tt c 2 u xx = cos x. ( 2 t c 2 2 x)u = cos x. v = ( t c x )u
Strauss PDEs e: Setion 3.4 - Exerise 3 Page 1 of 13 Exerise 3 Solve u tt = u xx + os x, u(x, ) = sin x, u t (x, ) = 1 + x. Solution Solution by Operator Fatorization Bring u xx to the other side. Write
More information