Sobol-Hoeffding Decomposition with Application to Global Sensitivity Analysis
|
|
- Maud Sims
- 5 years ago
- Views:
Transcription
1 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI Sobol-Hoeffding Decomposition with Application to Global Sensitivity Analysis Olivier Le Maître with Colleague & Friend Omar Knio LIMSI, CNRS UPR-325, Orsay, France PhD course on UQ - DTU
2 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI Table of content Sobol-Hoeffding decomposition Notations S-H decomposition 2 Application to Global SA Analysis of the variance (ANOVA) Sensitivity indices Example 3 Computation of the SI Monte-Carlo estimation SI from PC expansion
3 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI Notations L 2 functions over unit-hypercubes Let L 2 (U d ) be the space of real-valued squared-integrable functions over the d-dimensional hypercube U: ˆ f : x U d f (x) R, f L 2 (U d ) f (x) 2 dx < U d L 2 (U d ) is equipped with the inner product,, f, g L 2 (U d ), u, v := ˆ f (x)g(x)dx, U d and norm 2, f L 2 (U d ), f 2 := f, f /2
4 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI Notations L 2 functions over unit-hypercubes Let L 2 (U d ) be the space of real-valued squared-integrable functions over the d-dimensional hypercube U: ˆ f : x U d f (x) R, f L 2 (U d ) f (x) 2 dx < U d NB: all subsequent developments immediately extend to product-type situations, where x A = A A d R d, and weighted spaces L 2 (A, ρ), ρ : x A ρ(x) 0, ρ(x) = ρ (x ) ρ d (x d ) (eg: ρ is a pdf of a random vector x with mutually independent components)
5 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI Notations Ensemble notations Let D = {, 2,, d} Given i D, we denote i := D \ i its complement set in D, such that For instance i = {, 2} and i = {3,, d}, i = D and i = i i = D, i i =
6 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI Notations Ensemble notations Let D = {, 2,, d} Given i D, we denote i := D \ i its complement set in D, such that i i = D, i i = Given x = (x,, x d ), we denote x i the vector having for components the x i i, that is D i = {i,, i i } x i = (x i,, x i i ), where i := Card(i) For instance ˆ ˆ f (x)dx i = f (x,, x d ) dx i, U i U i i i and ˆ ˆ i / i f (x)dx i = f (x,, x d ) dx i, U d i U d i i D
7 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI S-H decomposition Sobol-Hoeffding decomposition Any f L 2 (U d ) has a unique hierarchical orthogonal decomposition of the form d d f (x) = f (x,, x d ) =f 0 + f i (x i ) + d i= d d i= j=i+ k=j+ i= j=i+ d f i,j (x i, x j )+ f i,j,k (x i, x j, x k ) + + f,,d (x,, x d ) Hierarchical: st order functionals (f i ) 2nd order functionals (f i,j ) 3rd order functionals (f i,j,l ) d-th order functional (f,,d ) Decomposition in a sum of 2 k functionals Using ensemble notations: f (x) = f i (x i ) i D
8 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI S-H decomposition Sobol-Hoeffding decomposition Any f L 2 (U d ) has a unique hierarchical orthogonal decomposition of the form f (x) = f i (x i ) i D Orthogonal: the functionals if the S-H decomposition verify the following orthogonality relations: ˆ f i (x i )dx j = 0, i D, j i, U ˆ f i (x i )f j (x j )dx = f i, f j = 0, i, j D, i j U d It follows the hierarchical construction ˆ f = f (x)dx = f =D U ˆ d f {i} = f (x)dx {i} f = f D\{i} f i D U ˆ d f i = f (x)dx i f j = f i f j i D, i 2 U i j i j i
9 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI S-H decomposition Table of content Sobol-Hoeffding decomposition Notations S-H decomposition 2 Application to Global SA Analysis of the variance (ANOVA) Sensitivity indices Example 3 Computation of the SI Monte-Carlo estimation SI from PC expansion
10 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI Analysis of the variance (ANOVA) Parametric sensitivity analysis Consider x as a set of d independent random parameters uniformly distributed on U d, and f (x) a model-output depending on these random parameters It is assumes that f is a 2nd order random variable: f L 2 (U d ) Thus, f has a unique S-H decomposition f (x) = f i (x i ) i D Further, the integrals of f with respect to i are in this context conditional expectations, ˆ E [f x i ] = f (x)dx i = g(x i ) U i i D, so the S-H decomposition follows the hierarchical structure f = E [f ] f {i} = E [ ] f x {i} E [f ] i D f i = E [f x i ] f j j i i D, i 2
11 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI Analysis of the variance (ANOVA) Variance decomposition Because of the orthogonality of the S-H decomposition the variance V [f ] of the model-output can be decomposed as i V [f ] = V [f i ], V [f i ] = f i, f i i D V [f i ] is interpreted as the contribution to the total variance V [f ] of the interaction between parameters x i i The S-H decomposition thus provide a rich mean of analyzing the respective contributions of individual or sets of parameters to model-output variability However, as there are 2 d contributions, so one needs more "abstract" characterizations
12 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI Analysis of the variance (ANOVA) Sensitivity indices To facilitate the hierarchization of the respective influence of each parameter x i, the partial variances V [f i ] are normalized by V [f ] to obtain the sensitivity indices: S i (f ) = V [f i] V [f ], i S i (f ) = i D The order of the sensitivity indices S i is equal to i = Card(i) st order sensitivity indices The d first order indices S {i} D characterize the fraction of the variance due the parameter x i only, ie without any interaction with others Therefore, d S {i} (f ) 0, i= measures globally the effect on the variability of all interactions between parameters If d i= S {i} =, the model is said additive, because its S-H decomposition is d f (x i,, x d ) = f 0 + f i (x i ), i= and the impact of the parameters can be studied separately
13 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI Sensitivity indices Sensitivity indices To facilitate the hierarchization of the respective influence of each parameter x i, the partial variances V [f i ] are normalized by V [f ] to obtain the sensitivity indices: S i (f ) = V [f i] V [f ], i S i (f ) = i D The order of the sensitivity indices S i is equal to i = Card(i) Total sensitivity indices The first order SI S {i} measures the variability due to parameter x i alone The total SI T {i} measures the variability due to the parameter x i, including all its interactions with other parameters: T {i} := i i S i S {i} Important point: for x i to be deemed non-important or non-influent on the model-output, S {i} and T {i} have to be negligible Observe that i D T {i}, the excess from characterizes the presence of interactions in the model-output
14 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI Sensitivity indices
15 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI Sensitivity indices Sensitivity indices In many uncertainty problem, the set of uncertain parameters can be naturally grouped into subsets depending on the process each parameter accounts for For instance, boundary conditions BC, material property ϕ, external forcing F, and D is the union of these distinct subsets: D = D BC D ϕ D F The notion of first order and total sensitivity indices can be extended to characterize the influence of the subsets of parameters For instance, S Dϕ = i D ϕ S i, measures the fraction of variance induced by the material uncertainty alone, while T DF = S i i D F measures the fraction of variance due to the external forcing uncertainty and all its interactions
16 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI Example Let (ξ, ξ 2 ) be two independent centered, normalized random variables ξ i N(0, ), i =, 2 Consider the model-output f : (ξ, ξ 2 ) R 2 R given by f (ξ, ξ 2 ) = (µ + σ ξ ) + (µ 2 + σ 2 ξ 2 ) Determine the S-H decomposition of f 2 Compute the st order and total sensitivity indices of f 3 Comment 4 Repeat for f (ξ, ξ 2 ) = (µ + σ ξ ) (µ 2 + σ 2 ξ 2 )
17 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI E [f ξ ] = (µ + µ 2 ) + σ ξ f (ξ ) = E [f ξ ] E [f ] = σ ξ E [f ξ 2 ] = (µ + µ 2 ) + σ 2 ξ 2 f 2 (ξ 2 ) = E [f ξ 2 ] E [f ] = σ 2 ξ 2 Example f (ξ, ξ 2 ) = µ + µ 2 + σ ξ + σ 2 ξ 2 E [f ] = (µ + µ 2 ) E [f ξ, ξ 2 ] = f (ξ, ξ 2 ) f,2 (ξ, ξ ) = E [f ξ, ξ 2 ] E [f ] f (ξ ) f 2 (ξ 2 ) = 0 Then, V [f ] = σ 2 + σ2 2, so σ 2 σ2 2 S = T = σ 2 + and S 2 = T 2 = σ2 2 σ 2 + σ2 2 Comment: obvious case, as f is a linear (additive) model
18 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI E [f ξ ] = µ µ 2 + µ 2 σ ξ f (ξ ) = E [f ξ ] E [f ] = µ 2 σ ξ E [f ξ 2 ] = µ µ 2 + µ σ 2 ξ 2 f 2 (ξ 2 ) = E [f ξ 2 ] E [f ] = µ σ 2 ξ 2 Example E [f ] = µ µ 2 f (ξ, ξ 2 ) = µ µ 2 + µ 2 σ ξ + µ σ 2 ξ 2 + σ σ 2 ξ ξ 2 E [f ξ, ξ 2 ] = f (ξ, ξ 2 ) f,2 (ξ, ξ ) = E [f ξ, ξ 2 ] E [f ] f (ξ ) f 2 (ξ 2 ) = σ σ 2 ξ ξ 2 Then, V [f ] = µ 2 σ2 2 + µ2 2 σ2 + σ2 σ2 2, so µ 2 2 S = σ2 µ 2 σ2 2 + µ2 2 σ2 + σ2 σ2 2 µ 2 and S 2 = σ2 2 µ 2 σ2 2 + µ2 2 σ2 +, σ2 σ2 2 µ 2 2 T = σ2 + σ σ 2 µ 2 µ 2 σ2 2 + µ2 2 σ2 + and T 2 = σ2 2 + σ σ 2 σ2 σ2 2 µ 2 σ2 2 + µ2 2 σ2 +, σ2 σ2 2 Comment: fraction of variance due to interactions is σ 2 σ2 2 /(µ2 σ2 2 + µ2 2 σ2 + σ2 σ2 2 )
19 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI Example f (ξ, ξ 2 ) = µ µ 2 + µ 2 σ ξ + µ σ 2 ξ 2 + σ σ 2 ξ ξ 2 Then, V [f ] = µ 2 σ2 2 + µ2 2 σ2 + σ2 σ2 2, so µ 2 2 S = σ2 µ 2 σ2 2 + µ2 2 σ2 + σ2 σ2 2 µ 2 and S 2 = σ2 2 µ 2 σ2 2 + µ2 2 σ2 +, σ2 σ2 2 µ 2 2 T = σ2 + σ σ 2 µ 2 µ 2 σ2 2 + µ2 2 σ2 + and T 2 = σ2 2 + σ σ 2 σ2 σ2 2 µ 2 σ2 2 + µ2 2 σ2 +, σ2 σ2 2 Comment: fraction of variance due to interactions is σ 2 σ2 2 /(µ2 σ2 2 + µ2 2 σ2 + σ2 σ2 2 ) Example: (µ, σ ) = (, 3) and (µ 2, σ 2 ) = (2, 2), so S = 9/9, S 2 = /9, S,2 = 9/9 One can draw the conclusions: ξ is the most influential variable as S > S 2 and T > T 2 interactions are important as S S 2 = 9/9 05, especially for ξ 2 for which (T 2 S 2 )/T 2 = 9/0
20 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI Example Table of content Sobol-Hoeffding decomposition Notations S-H decomposition 2 Application to Global SA Analysis of the variance (ANOVA) Sensitivity indices Example 3 Computation of the SI Monte-Carlo estimation SI from PC expansion
21 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI Monte-Carlo estimation -st order sensitivity indices by Monte-Carlo sampling The S i can be computed by MC sampling as follow Recall that S {i} (f ) = V [ f {i} ] V [f ] Observe: E [E [f x i ]] = E [f ] = V [E [f x i ])] V [f ] E = [ E [f x i ] 2] E [E [f x i ]] 2 E [f ] and V [f ] can be estimated using MC sampling Let χ M = {x (),, x (M) } be a set of independent samples drawn uniformly in U d, the mean and variance estimators are: V [f ] Ê [f ] = M M l= f (x (l) ), V [f ] = M M f (x (l) ) 2 Ê [f ]2 l= It now remains to compute the variance of conditional expectations, V [ E [ f x {i} ]] Any idea?
22 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI Monte-Carlo estimation Monte-Carlo estimator for variance of conditional expectation Observe: ˆ f (x i, x i )f (x i, x U d+ i i )dx idx i dx i ˆ ˆ ˆ = dx U i i f (x i, x i )dx i f (x i, x U i U i i )dx i ˆ [ˆ ] 2 = dx U i i f (x i, x i )dx i U i
23 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI Monte-Carlo estimation Monte-Carlo estimator for variance of conditional expectation Observe: ˆ f (x i, x i )f (x i, x U d+ i i )dx idx i dx i ˆ ˆ ˆ = dx U i i f (x i, x i )dx i f (x i, x U i U i i )dx i ˆ [ˆ ] 2 = dx U i i f (x i, x i )dx i U i V [ E [ ]] [ f x {i} = E E [ ] ] 2 f x {i} E [ E [ ]] 2 f x {i} = E [E [ ] ] 2 f x {i} E [f ] 2
24 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI Monte-Carlo estimation Monte-Carlo estimator for variance of conditional expectation Observe: ˆ f (x i, x i )f (x i, x U d+ i i )dx idx i dx i ˆ ˆ ˆ = dx U i i f (x i, x i )dx i f (x i, x U i U i i )dx i ˆ [ˆ ] 2 = dx U i i f (x i, x i )dx i U i V [ E [ ]] [ f x {i} = E E [ ] ] 2 f x {i} E [ E [ ]] 2 f x {i} = E [E [ ] ] 2 f x {i} E [f ] 2 = lim M M M ( ) ( ) f x (l) {i}, x (l) f x (l) {i} {i}, x (l) E [f ] 2 {i} l= ( independent samples x (l) (l), x {i} {i} ) (k) and x {i}
25 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI Monte-Carlo estimation Monte-Carlo estimators for st order SI S {i} Draw 2 independent sample sets, with size M, χ M and χ M 2 Compute estimators Ê [f ] and V [f ] from χ M (or χ M ) [M model evaluations] 3 For i =, 2,, d: Estimate variance of conditional expectation through V [ E [ ]] f x {i} = M Estimator of the -st order SI: M l= ( ) ( ) f x (l) {i}, x (l) f x (l) {i} {i}, x (l) Ê [f ]2 {i} V [ E [ ]] f ξ {i} Ŝ {i} (f ) = V [f ] [M new model evaluations at ( x (l) {i}, x (l) {i} )] Requires (d + ) M model evaluations
26 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI Monte-Carlo estimation χ M : x () x () i x () d x (l) x (l) i x (l) d x (M) x (M) i x (M) d χ M : x () x () i x () d x (l) x (l) i x (l) d x (M) x (M) i x (M) d χ {i} M : x () x () i x () d x (l) x (l) i x (l) d x (M) x (M) i x (M) d V [ E [ f x {i} ]] M M l= f (x {i}, x {i} )f ( x {i}, x {i} ) ( M M l= f (x {i}, x {i} )) 2
27 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI Monte-Carlo estimation Total sensitivity indices by Monte-Carlo sampling The T {i} can be computed by MC sampling as follow Recall that T {i} (f ) = S i (f ) = i {i} i D\{i} S i (f ) = V [ E [ f x {i} ]] V [f ] As for V [ E [ f x {i} ]], we can derive the following Monte-Carlo estimator for V [ E [ f x {i} ]], using the two independent sample sets χm and χ M V [ E [ ]] f x {i} = M M f l= ( ) ( ) x (l), x (l) f x (l), x (l) E [f ] 2, {i} {i} {i} {i} so finally ( T {i} (f ) = V [f ] M M f l= ) ( ) ( ) x (l), x (l) f x (l), x (l) E [f ] 2 {i} {i} {i} {i} The MC estimation of the T {i} needs d M additional model evaluations
28 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI Monte-Carlo estimation χ M : x () x () i x () d x (l) x (l) i x (l) d x (M) x (M) i x (M) d χ M : x () x () i x () d x (l) x (l) i x (l) d x (M) x (M) i x (M) d χ {i} M : x () x () i x () d x (l) x (l) i x (l) d x (M) x (M) i x (M) d V [ E [ f x {i} ]] M M l= f (x {i}, x {i} )f (x {i}, x {i} ) ( M M l= f (x {i}, x {i} )) 2
29 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI Monte-Carlo estimation MC estimation of st order and total sensitivity indices: requires M (2d + ) simulations convergence of estimators in O(/ M) slow convergence, but d-independent convergence not related to smoothness of f works also (in practice) with advanced sampling schemes (QMC, LHS, ) to proceed with QMC/LHS? How
30 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI SI from PC expansion S-H decomposition from PC expansions Consider the model-output f : ξ Ξ R d R, where ξ = (ξ,, ξ d ) are independent real-valued rv with joint-probability density function p ξ (x,, x d ) = d p i (x i ) Let { Ψα} be the set of d-variate orthogonal polynomials, Ψα(ξ) = d i= i= ψ (i) α i (ξ i ), with ψ (i) l 0 π l the uni-variate polynomials mutually orthogonal with respect to the density p i If f L 2 (Ξ, p ξ ), it has a convergent PC expansion d f (ξ) = lim Ψα(ξ)fα, α = α i No α No i=
31 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI SI from PC expansion S-H decomposition from PC expansions Given the truncated PC expansion of f, ˆf (ξ) = Ψα(ξ)fα, α A one can readily obtain the PC approximation of the S-H functionals through ˆfi (ξ i ) = Ψα(ξ i )fα, where the multi-index set A(i) is given by α A(i) A(i) = {α A; α i > 0 for i i, α i = 0 for i / i} A For the sensitivity indices it comes α A(i) S i (ˆf ) = f α 2 Ψα, Ψα α A f α 2 Ψα, Ψα, T α T (i) {i}(ˆf ) = f α 2 Ψα, Ψα α A f α 2, Ψα, Ψα where T (i) = {α A; α i > 0 for i i}
32 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI SI from PC expansion Subsets of PC multi-indices for S-H functionals
33 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI SI from PC expansion Example: Let (ξ, ξ 2 ) be two independent random variables with uniform distributions on the unit interval Consider the model-output f : (ξ, ξ 2 ) [0, ] 2 R given by f (ξ, ξ 2 ) = tanh(3(ξ 02ξ 2 5)( + ξ 2 )) f(ξ, ξ 2 ) ξ ξ 08
34 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI SI from PC expansion S-H decomposition at No = 3 (P = 0): ξ ξ ξ ξ ξ ξ ξ ξ 08
35 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI SI from PC expansion S-H decomposition at No = 5: (P = 2) ξ ξ ξ ξ ξ ξ ξ ξ 08
36 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI SI from PC expansion S-H decomposition at No = 7: (P = 36) ξ ξ ξ ξ ξ ξ ξ ξ 08
37 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI SI from PC expansion S-H decomposition at No = 9: (P = 55) ξ ξ ξ ξ ξ ξ ξ ξ 08
38 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI SI from PC expansion Convergence of sensitivity indices: No Card(A) S {} (ˆf ) S {2} (ˆf ) S {,2} (ˆf ) ˆf f ( ) 2567 ( 2) 390 ( 2) 7 ( 2) ( ) 2545 ( 2) 676 ( 2) 26 ( 2) ( ) 2542 ( 2) 734 ( 2) 96 ( 3) ( ) 2542 ( 2) 745 ( 2) 36 ( 3) ( ) 2542 ( 2) 748 ( 2) 9 ( 4)
39 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI SI from PC expansion Convergence of sensitivity indices for naive MC sampling: S 0956 S e+06 e+07 M e+06 e+07 M
40 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI SI from PC expansion Homework: Homma-Saltelli function f (x, x 2, x 3 ) = sin(x ) + a sin 2 (x 2 ) + bx 4 3 sin(x ), where x, x 2 and x 3 are iid random variables uniformly distributed on [ π, π] Letting be the expectation (over [ π, π] 3 ), we have sin = 0, sin 2 = /2, sin 4 = 3/8, x n = { π n n+, n even, 0, n odd Compute the SH decomposition of f and the associated sensitivity indices Solution: f (x, x 2, x 3 ) = a/2+sin(x )(+bπ 2 /4)+a(sin 2 (x 2 ) /2)+b sin(x )(x 4 3 π4 /5)) 2 Compute PC expansion of f by non-intrusive spectral projection for polynomial degrees to 6 and extract st order and total sensitivity indices for a = 7 and b = 0 3 Compute -st order and total sensitivity indices by Monte-Carlo simulation
41 Sobol-Hoeffding decomposition Application to Global SA Computation of the SI SI from PC expansion Further readings: Hoeffding, W (948) A class of statistics with asymptotically normal distribution The annals of Mathematical Statistics, 9, pp Sobol, I M (993) Sensitivity estimates for nonlinear mathematical models Wiley,, pp Sobol, I M (200) Global sensitivity indices for nonlinear mathematical models and their Monte Carlo estimates Mathematics and Computers in Simulations, 55, pp Homma T, Saltelli A, (996) Importance measures in globla sensitivity analysis of nonlinear models Reliab Eng Syst Safety, 52:, pp -7 Crestaux T, Le Maître O and Martinez JM, (2009) Polynomial Chaos expansion for sensitivity analysis Reliabg Eng Syst Safety, 94:7, pp 6-72
Local and Global Sensitivity Analysis
Omar 1,2 1 Duke University Department of Mechanical Engineering & Materials Science omar.knio@duke.edu 2 KAUST Division of Computer, Electrical, Mathematical Science & Engineering omar.knio@kaust.edu.sa
More informationPC EXPANSION FOR GLOBAL SENSITIVITY ANALYSIS OF NON-SMOOTH FUNCTIONALS OF UNCERTAIN STOCHASTIC DIFFERENTIAL EQUATIONS SOLUTIONS
PC EXPANSION FOR GLOBAL SENSITIVITY ANALYSIS OF NON-SMOOTH FUNCTIONALS OF UNCERTAIN STOCHASTIC DIFFERENTIAL EQUATIONS SOLUTIONS M. Navarro, O.P. Le Maître,2, O.M. Knio,3 mariaisabel.navarrojimenez@kaust.edu.sa
More informationPolynomial chaos expansions for sensitivity analysis
c DEPARTMENT OF CIVIL, ENVIRONMENTAL AND GEOMATIC ENGINEERING CHAIR OF RISK, SAFETY & UNCERTAINTY QUANTIFICATION Polynomial chaos expansions for sensitivity analysis B. Sudret Chair of Risk, Safety & Uncertainty
More informationSENSITIVITY ANALYSIS IN NUMERICAL SIMULATION OF MULTIPHASE FLOW FOR CO 2 STORAGE IN SALINE AQUIFERS USING THE PROBABILISTIC COLLOCATION APPROACH
XIX International Conference on Water Resources CMWR 2012 University of Illinois at Urbana-Champaign June 17-22,2012 SENSITIVITY ANALYSIS IN NUMERICAL SIMULATION OF MULTIPHASE FLOW FOR CO 2 STORAGE IN
More informationAnalysis of covariance (ANCOVA) using polynomial chaos expansions
Research Collection Conference Paper Analysis of covariance (ANCOVA) using polynomial chaos expansions Author(s): Sudret, Bruno; Caniou, Yves Publication Date: 013 Permanent Link: https://doi.org/10.399/ethz-a-0100633
More informationEstimating functional uncertainty using polynomial chaos and adjoint equations
0. Estimating functional uncertainty using polynomial chaos and adjoint equations February 24, 2011 1 Florida State University, Tallahassee, Florida, Usa 2 Moscow Institute of Physics and Technology, Moscow,
More informationAlgorithms for Uncertainty Quantification
Algorithms for Uncertainty Quantification Lecture 9: Sensitivity Analysis ST 2018 Tobias Neckel Scientific Computing in Computer Science TUM Repetition of Previous Lecture Sparse grids in Uncertainty Quantification
More informationGaussian Process Regression and Emulation
Gaussian Process Regression and Emulation STAT8810, Fall 2017 M.T. Pratola September 22, 2017 Today Experimental Design; Sensitivity Analysis Designing Your Experiment If you will run a simulator model,
More informationPolynomial Chaos and Karhunen-Loeve Expansion
Polynomial Chaos and Karhunen-Loeve Expansion 1) Random Variables Consider a system that is modeled by R = M(x, t, X) where X is a random variable. We are interested in determining the probability of the
More informationDerivative-based global sensitivity measures for interactions
Derivative-based global sensitivity measures for interactions Olivier Roustant École Nationale Supérieure des Mines de Saint-Étienne Joint work with J. Fruth, B. Iooss and S. Kuhnt SAMO 13 Olivier Roustant
More informationA Unified Framework for Uncertainty and Sensitivity Analysis of Computational Models with Many Input Parameters
A Unified Framework for Uncertainty and Sensitivity Analysis of Computational Models with Many Input Parameters C. F. Jeff Wu H. Milton Stewart School of Industrial and Systems Engineering Georgia Institute
More informationConcepts in Global Sensitivity Analysis IMA UQ Short Course, June 23, 2015
Concepts in Global Sensitivity Analysis IMA UQ Short Course, June 23, 2015 A good reference is Global Sensitivity Analysis: The Primer. Saltelli, et al. (2008) Paul Constantine Colorado School of Mines
More informationGaussian process regression for Sensitivity analysis
Gaussian process regression for Sensitivity analysis GPSS Workshop on UQ, Sheffield, September 2016 Nicolas Durrande, Mines St-Étienne, durrande@emse.fr GPSS workshop on UQ GPs for sensitivity analysis
More informationwhere r n = dn+1 x(t)
Random Variables Overview Probability Random variables Transforms of pdfs Moments and cumulants Useful distributions Random vectors Linear transformations of random vectors The multivariate normal distribution
More informationStochastic Spectral Approaches to Bayesian Inference
Stochastic Spectral Approaches to Bayesian Inference Prof. Nathan L. Gibson Department of Mathematics Applied Mathematics and Computation Seminar March 4, 2011 Prof. Gibson (OSU) Spectral Approaches to
More informationIntroduction to Smoothing spline ANOVA models (metamodelling)
Introduction to Smoothing spline ANOVA models (metamodelling) M. Ratto DYNARE Summer School, Paris, June 215. Joint Research Centre www.jrc.ec.europa.eu Serving society Stimulating innovation Supporting
More informationUncertainty Quantification in Computational Science
DTU 2010 - Lecture I Uncertainty Quantification in Computational Science Jan S Hesthaven Brown University Jan.Hesthaven@Brown.edu Objective of lectures The main objective of these lectures are To offer
More informationA reduced-order stochastic finite element analysis for structures with uncertainties
A reduced-order stochastic finite element analysis for structures with uncertainties Ji Yang 1, Béatrice Faverjon 1,2, Herwig Peters 1, icole Kessissoglou 1 1 School of Mechanical and Manufacturing Engineering,
More informationIntroduction to Uncertainty Quantification in Computational Science Handout #3
Introduction to Uncertainty Quantification in Computational Science Handout #3 Gianluca Iaccarino Department of Mechanical Engineering Stanford University June 29 - July 1, 2009 Scuola di Dottorato di
More informationPolynomial chaos expansions for structural reliability analysis
DEPARTMENT OF CIVIL, ENVIRONMENTAL AND GEOMATIC ENGINEERING CHAIR OF RISK, SAFETY & UNCERTAINTY QUANTIFICATION Polynomial chaos expansions for structural reliability analysis B. Sudret & S. Marelli Incl.
More informationDynamic System Identification using HDMR-Bayesian Technique
Dynamic System Identification using HDMR-Bayesian Technique *Shereena O A 1) and Dr. B N Rao 2) 1), 2) Department of Civil Engineering, IIT Madras, Chennai 600036, Tamil Nadu, India 1) ce14d020@smail.iitm.ac.in
More informationEXPOSITORY NOTES ON DISTRIBUTION THEORY, FALL 2018
EXPOSITORY NOTES ON DISTRIBUTION THEORY, FALL 2018 While these notes are under construction, I expect there will be many typos. The main reference for this is volume 1 of Hörmander, The analysis of liner
More informationSparse polynomial chaos expansions in engineering applications
DEPARTMENT OF CIVIL, ENVIRONMENTAL AND GEOMATIC ENGINEERING CHAIR OF RISK, SAFETY & UNCERTAINTY QUANTIFICATION Sparse polynomial chaos expansions in engineering applications B. Sudret G. Blatman (EDF R&D,
More informationGlobal Sensitivity Analysis in Structural Optimization
Global Sensitivity Analysis in Structural Optimization Uwe Reuter Department of Civil Engineering, TU Dresden, Germany Martin Liebscher and Heiner Müllerschön DYNAmore GmbH, Stuttgart, Germany Summary:
More informationSparse Pseudo Spectral Projection Methods with Directional Adaptation for Uncertainty Quantification
Noname manuscript No. (will be inserted by the editor) Sparse Pseudo Spectral Projection Methods with Directional Adaptation for Uncertainty Quantification J. Winokur D. Kim F. Bisetti O.P. Le Maître O.M.
More informationLectures. Variance-based sensitivity analysis in the presence of correlated input variables. Thomas Most. Source:
Lectures Variance-based sensitivity analysis in the presence of correlated input variables Thomas Most Source: www.dynardo.de/en/library Variance-based sensitivity analysis in the presence of correlated
More information18Ï È² 7( &: ÄuANOVAp.O`û5 571 Based on this ANOVA model representation, Sobol (1993) proposed global sensitivity index, S i1...i s = D i1...i s /D, w
A^VÇÚO 1 Êò 18Ï 2013c12 Chinese Journal of Applied Probability and Statistics Vol.29 No.6 Dec. 2013 Optimal Properties of Orthogonal Arrays Based on ANOVA High-Dimensional Model Representation Chen Xueping
More informationPerformance Evaluation of Generalized Polynomial Chaos
Performance Evaluation of Generalized Polynomial Chaos Dongbin Xiu, Didier Lucor, C.-H. Su, and George Em Karniadakis 1 Division of Applied Mathematics, Brown University, Providence, RI 02912, USA, gk@dam.brown.edu
More informationFourier Sin and Cos Series and Least Squares Convergence
Fourier and east Squares Convergence James K. Peterson Department of Biological Sciences and Department of Mathematical Sciences Clemson University May 7, 28 Outline et s look at the original Fourier sin
More informationKeywords: Sonic boom analysis, Atmospheric uncertainties, Uncertainty quantification, Monte Carlo method, Polynomial chaos method.
Blucher Mechanical Engineering Proceedings May 2014, vol. 1, num. 1 www.proceedings.blucher.com.br/evento/10wccm SONIC BOOM ANALYSIS UNDER ATMOSPHERIC UNCERTAINTIES BY A NON-INTRUSIVE POLYNOMIAL CHAOS
More informationLearning gradients: prescriptive models
Department of Statistical Science Institute for Genome Sciences & Policy Department of Computer Science Duke University May 11, 2007 Relevant papers Learning Coordinate Covariances via Gradients. Sayan
More informationDinesh Kumar, Mehrdad Raisee and Chris Lacor
Dinesh Kumar, Mehrdad Raisee and Chris Lacor Fluid Mechanics and Thermodynamics Research Group Vrije Universiteit Brussel, BELGIUM dkumar@vub.ac.be; m_raisee@yahoo.com; chris.lacor@vub.ac.be October, 2014
More informationPractical unbiased Monte Carlo for Uncertainty Quantification
Practical unbiased Monte Carlo for Uncertainty Quantification Sergios Agapiou Department of Statistics, University of Warwick MiR@W day: Uncertainty in Complex Computer Models, 2nd February 2015, University
More informationCollocation based high dimensional model representation for stochastic partial differential equations
Collocation based high dimensional model representation for stochastic partial differential equations S Adhikari 1 1 Swansea University, UK ECCM 2010: IV European Conference on Computational Mechanics,
More informationMath Real Analysis II
Math 4 - Real Analysis II Solutions to Homework due May Recall that a function f is called even if f( x) = f(x) and called odd if f( x) = f(x) for all x. We saw that these classes of functions had a particularly
More informationSingular Integrals. 1 Calderon-Zygmund decomposition
Singular Integrals Analysis III Calderon-Zygmund decomposition Let f be an integrable function f dx 0, f = g + b with g Cα almost everywhere, with b
More informationThe following definition is fundamental.
1. Some Basics from Linear Algebra With these notes, I will try and clarify certain topics that I only quickly mention in class. First and foremost, I will assume that you are familiar with many basic
More informationAdvanced Algorithms for Multidimensional Sensitivity Studies of Large-scale Air Pollution Models Based on Sobol Sequences
Advanced Algorithms for Multidimensional Sensitivity Studies of Large-scale Air Pollution Models Based on Sobol Sequences I. Dimov a,, R. Georgieva a, Tz. Ostromsky a, Z. Zlatev b a Department of Parallel
More informationInformation Science 2
Information Science 2 Probability Theory: An Overview Week 12 College of Information Science and Engineering Ritsumeikan University Agenda Terms and concepts from Week 11 Basic concepts of probability
More informationStochastic Collocation Methods for Polynomial Chaos: Analysis and Applications
Stochastic Collocation Methods for Polynomial Chaos: Analysis and Applications Dongbin Xiu Department of Mathematics, Purdue University Support: AFOSR FA955-8-1-353 (Computational Math) SF CAREER DMS-64535
More informationVectors in Function Spaces
Jim Lambers MAT 66 Spring Semester 15-16 Lecture 18 Notes These notes correspond to Section 6.3 in the text. Vectors in Function Spaces We begin with some necessary terminology. A vector space V, also
More informationMetric Spaces. Exercises Fall 2017 Lecturer: Viveka Erlandsson. Written by M.van den Berg
Metric Spaces Exercises Fall 2017 Lecturer: Viveka Erlandsson Written by M.van den Berg School of Mathematics University of Bristol BS8 1TW Bristol, UK 1 Exercises. 1. Let X be a non-empty set, and suppose
More informationBeyond Wiener Askey Expansions: Handling Arbitrary PDFs
Journal of Scientific Computing, Vol. 27, Nos. 1 3, June 2006 ( 2005) DOI: 10.1007/s10915-005-9038-8 Beyond Wiener Askey Expansions: Handling Arbitrary PDFs Xiaoliang Wan 1 and George Em Karniadakis 1
More informationEFFICIENT SHAPE OPTIMIZATION USING POLYNOMIAL CHAOS EXPANSION AND LOCAL SENSITIVITIES
9 th ASCE Specialty Conference on Probabilistic Mechanics and Structural Reliability EFFICIENT SHAPE OPTIMIZATION USING POLYNOMIAL CHAOS EXPANSION AND LOCAL SENSITIVITIES Nam H. Kim and Haoyu Wang University
More informationChapter 4: Interpolation and Approximation. October 28, 2005
Chapter 4: Interpolation and Approximation October 28, 2005 Outline 1 2.4 Linear Interpolation 2 4.1 Lagrange Interpolation 3 4.2 Newton Interpolation and Divided Differences 4 4.3 Interpolation Error
More informationMonte Carlo Methods in Statistical Mechanics
Monte Carlo Methods in Statistical Mechanics Mario G. Del Pópolo Atomistic Simulation Centre School of Mathematics and Physics Queen s University Belfast Belfast Mario G. Del Pópolo Statistical Mechanics
More informationSafety Envelope for Load Tolerance and Its Application to Fatigue Reliability Design
Safety Envelope for Load Tolerance and Its Application to Fatigue Reliability Design Haoyu Wang * and Nam H. Kim University of Florida, Gainesville, FL 32611 Yoon-Jun Kim Caterpillar Inc., Peoria, IL 61656
More informationStochastic Spectral Methods for Uncertainty Quantification
Stochastic Spectral Methods for Uncertainty Quantification Olivier Le Maître 1,2,3, Omar Knio 1,2 1- Duke University, Durham, North Carolina 2- KAUST, Saudi-Arabia 3- LIMSI CNRS, Orsay, France 38th Conf.
More informationApplication and validation of polynomial chaos methods to quantify uncertainties in simulating the Gulf of Mexico circulation using HYCOM.
Application and validation of polynomial chaos methods to quantify uncertainties in simulating the Gulf of Mexico circulation using HYCOM. Mohamed Iskandarani Matthieu Le Hénaff Carlisle Thacker University
More informationSensitivity analysis of parametric uncertainties and modeling errors. in generalized probabilistic modeling
ECCOMAS European Congress on Computational Methods in Applied Sciences and Engineering Sensitivity analysis of parametric uncertainties and modeling errors in generalized probabilistic modeling Maarten
More informationChaospy: A modular implementation of Polynomial Chaos expansions and Monte Carlo methods
Chaospy: A modular implementation of Polynomial Chaos expansions and Monte Carlo methods Simen Tennøe Supervisors: Jonathan Feinberg Hans Petter Langtangen Gaute Einevoll Geir Halnes University of Oslo,
More informationarxiv: v2 [math.pr] 27 Oct 2015
A brief note on the Karhunen-Loève expansion Alen Alexanderian arxiv:1509.07526v2 [math.pr] 27 Oct 2015 October 28, 2015 Abstract We provide a detailed derivation of the Karhunen Loève expansion of a stochastic
More informationAsymptotic normality of global clustering coefficient of uniform random intersection graph
Asymptotic normality of global clustering coefficient of uniform random intersection graph Mindaugas Bloznelis joint work with Jerzy Jaworski Vilnius University, Lithuania http://wwwmifvult/ bloznelis
More informationIntroduction to uncertainty quantification An example of application in medicine
Introduction to uncertainty quantification An example of application in medicine Laurent Dumas Laboratoire de Mathématiques de Versailles (LMV) Versailles University Short course, University of Mauritius,
More informationRandom Eigenvalue Problems Revisited
Random Eigenvalue Problems Revisited S Adhikari Department of Aerospace Engineering, University of Bristol, Bristol, U.K. Email: S.Adhikari@bristol.ac.uk URL: http://www.aer.bris.ac.uk/contact/academic/adhikari/home.html
More informationAn Introduction to Wavelets and some Applications
An Introduction to Wavelets and some Applications Milan, May 2003 Anestis Antoniadis Laboratoire IMAG-LMC University Joseph Fourier Grenoble, France An Introduction to Wavelets and some Applications p.1/54
More informationRandom Matrix Eigenvalue Problems in Probabilistic Structural Mechanics
Random Matrix Eigenvalue Problems in Probabilistic Structural Mechanics S Adhikari Department of Aerospace Engineering, University of Bristol, Bristol, U.K. URL: http://www.aer.bris.ac.uk/contact/academic/adhikari/home.html
More informationReal Analysis Problems
Real Analysis Problems Cristian E. Gutiérrez September 14, 29 1 1 CONTINUITY 1 Continuity Problem 1.1 Let r n be the sequence of rational numbers and Prove that f(x) = 1. f is continuous on the irrationals.
More information8 STOCHASTIC SIMULATION
8 STOCHASTIC SIMULATIO 59 8 STOCHASTIC SIMULATIO Whereas in optimization we seek a set of parameters x to minimize a cost, or to maximize a reward function J( x), here we pose a related but different question.
More informationELEG 3143 Probability & Stochastic Process Ch. 6 Stochastic Process
Department of Electrical Engineering University of Arkansas ELEG 3143 Probability & Stochastic Process Ch. 6 Stochastic Process Dr. Jingxian Wu wuj@uark.edu OUTLINE 2 Definition of stochastic process (random
More informationRECOMMENDED TOOLS FOR SENSITIVITY ANALYSIS ASSOCIATED TO THE EVALUATION OF MEASUREMENT UNCERTAINTY
Advanced Mathematical and Computational Tools in Metrology and Testing IX Edited by F. Pavese, M. Bär, J.-R. Filtz, A. B. Forbes, L. Pendrill and K. Shirono c 2012 World Scientific Publishing Company (pp.
More informationUncertainty Propagation and Global Sensitivity Analysis in Hybrid Simulation using Polynomial Chaos Expansion
Uncertainty Propagation and Global Sensitivity Analysis in Hybrid Simulation using Polynomial Chaos Expansion EU-US-Asia workshop on hybrid testing Ispra, 5-6 October 2015 G. Abbiati, S. Marelli, O.S.
More informationMeasure and Integration: Solutions of CW2
Measure and Integration: s of CW2 Fall 206 [G. Holzegel] December 9, 206 Problem of Sheet 5 a) Left (f n ) and (g n ) be sequences of integrable functions with f n (x) f (x) and g n (x) g (x) for almost
More informationProbability Theory and Statistics. Peter Jochumzen
Probability Theory and Statistics Peter Jochumzen April 18, 2016 Contents 1 Probability Theory And Statistics 3 1.1 Experiment, Outcome and Event................................ 3 1.2 Probability............................................
More informationQuasi-Monte Carlo integration over the Euclidean space and applications
Quasi-Monte Carlo integration over the Euclidean space and applications f.kuo@unsw.edu.au University of New South Wales, Sydney, Australia joint work with James Nichols (UNSW) Journal of Complexity 30
More informationHowever, reliability analysis is not limited to calculation of the probability of failure.
Probabilistic Analysis probabilistic analysis methods, including the first and second-order reliability methods, Monte Carlo simulation, Importance sampling, Latin Hypercube sampling, and stochastic expansions
More information4th Preparation Sheet - Solutions
Prof. Dr. Rainer Dahlhaus Probability Theory Summer term 017 4th Preparation Sheet - Solutions Remark: Throughout the exercise sheet we use the two equivalent definitions of separability of a metric space
More informationCourse: ESO-209 Home Work: 1 Instructor: Debasis Kundu
Home Work: 1 1. Describe the sample space when a coin is tossed (a) once, (b) three times, (c) n times, (d) an infinite number of times. 2. A coin is tossed until for the first time the same result appear
More informationON MATRIX VALUED SQUARE INTEGRABLE POSITIVE DEFINITE FUNCTIONS
1 2 3 ON MATRIX VALUED SQUARE INTERABLE POSITIVE DEFINITE FUNCTIONS HONYU HE Abstract. In this paper, we study matrix valued positive definite functions on a unimodular group. We generalize two important
More informationAn introduction to Birkhoff normal form
An introduction to Birkhoff normal form Dario Bambusi Dipartimento di Matematica, Universitá di Milano via Saldini 50, 0133 Milano (Italy) 19.11.14 1 Introduction The aim of this note is to present an
More informationA framework for global reliability sensitivity analysis in the presence of multi-uncertainty
A framework for global reliability sensitivity analysis in the presence of multi-uncertainty Max Ehre, Iason Papaioannou, Daniel Straub Engineering Risk Analysis Group, Technische Universität München,
More informationGeneralized Sobol indices for dependent variables
Generalized Sobol indices for dependent variables Gaelle Chastaing Fabrice Gamboa Clémentine Prieur July 1st, 2013 1/28 Gaelle Chastaing Sensitivity analysis and dependent variables Contents 1 Context
More informationComputational methods for uncertainty quantification and sensitivity analysis of complex systems
DEPARTMENT OF CIVIL, ENVIRONMENTAL AND GEOMATIC ENGINEERING CHAIR OF RISK, SAFETY & UNCERTAINTY QUANTIFICATION Computational methods for uncertainty quantification and sensitivity analysis of complex systems
More informationUncertainty Quantification of an ORC turbine blade under a low quantile constrain
Available online at www.sciencedirect.com ScienceDirect Energy Procedia 129 (2017) 1149 1155 www.elsevier.com/locate/procedia IV International Seminar on ORC Power Systems, ORC2017 13-15 September 2017,
More informationUncertainty representation and propagation in chemical networks
Uncertainty representation and propagation in chemical networks P. Pernot Laboratoire de Chimie Physique, Orsay Plan 1 Introduction 2 Uncertainty propagation 3 Treatment of uncertainties for chemical rates
More informationPredictive Engineering and Computational Sciences. Local Sensitivity Derivative Enhanced Monte Carlo Methods. Roy H. Stogner, Vikram Garg
PECOS Predictive Engineering and Computational Sciences Local Sensitivity Derivative Enhanced Monte Carlo Methods Roy H. Stogner, Vikram Garg Institute for Computational Engineering and Sciences The University
More informationSensitivity Analysis of Compact Models in Nanodevice Modeling
Sensitivity Analysis of Compact Models in Nanodevice Modeling 1 (in collaboration with Asen Asenov 2 Binjie Cheng 2 Ivan Dimov 1 Urban Kovac 2 Campbell Millar 2 ) 1 Institute of Information and Communication
More informationSome methods for sensitivity analysis of systems / networks
Article Some methods for sensitivity analysis of systems / networks WenJun Zhang School of Life Sciences, Sun Yat-sen University, Guangzhou 510275, China; International Academy of Ecology and Environmental
More informationAsymptotics of minimax stochastic programs
Asymptotics of minimax stochastic programs Alexander Shapiro Abstract. We discuss in this paper asymptotics of the sample average approximation (SAA) of the optimal value of a minimax stochastic programming
More informationIf we want to analyze experimental or simulated data we might encounter the following tasks:
Chapter 1 Introduction If we want to analyze experimental or simulated data we might encounter the following tasks: Characterization of the source of the signal and diagnosis Studying dependencies Prediction
More informationUncertainty Quantification in Computational Models
Uncertainty Quantification in Computational Models Habib N. Najm Sandia National Laboratories, Livermore, CA, USA Workshop on Understanding Climate Change from Data (UCC11) University of Minnesota, Minneapolis,
More informationOutline of Fourier Series: Math 201B
Outline of Fourier Series: Math 201B February 24, 2011 1 Functions and convolutions 1.1 Periodic functions Periodic functions. Let = R/(2πZ) denote the circle, or onedimensional torus. A function f : C
More information7: FOURIER SERIES STEVEN HEILMAN
7: FOURIER SERIES STEVE HEILMA Contents 1. Review 1 2. Introduction 1 3. Periodic Functions 2 4. Inner Products on Periodic Functions 3 5. Trigonometric Polynomials 5 6. Periodic Convolutions 7 7. Fourier
More information5. Some theorems on continuous functions
5. Some theorems on continuous functions The results of section 3 were largely concerned with continuity of functions at a single point (usually called x 0 ). In this section, we present some consequences
More informationSample Quantum Chemistry Exam 2 Solutions
Chemistry 46 Fall 7 Dr. Jean M. Standard Name SAMPE EXAM Sample Quantum Chemistry Exam Solutions.) ( points) Answer the following questions by selecting the correct answer from the choices provided. a.)
More informationHochdimensionale Integration
Oliver Ernst Institut für Numerische Mathematik und Optimierung Hochdimensionale Integration 14-tägige Vorlesung im Wintersemester 2010/11 im Rahmen des Moduls Ausgewählte Kapitel der Numerik Contents
More informationFast Numerical Methods for Stochastic Computations
Fast AreviewbyDongbinXiu May 16 th,2013 Outline Motivation 1 Motivation 2 3 4 5 Example: Burgers Equation Let us consider the Burger s equation: u t + uu x = νu xx, x [ 1, 1] u( 1) =1 u(1) = 1 Example:
More informationThe stochastic modeling
February 211 Gdansk A schedule of the lecture Generation of randomly arranged numbers Monte Carlo Strong Law of Large Numbers SLLN We use SAS 9 For convenience programs will be given in an appendix Programs
More informationStochastic representation of random positive-definite tensor-valued properties: application to 3D anisotropic permeability random fields
Sous la co-tutelle de : LABORATOIRE DE MODÉLISATION ET SIMULATION MULTI ÉCHELLE CNRS UPEC UNIVERSITÉ PARIS-EST CRÉTEIL UPEM UNIVERSITÉ PARIS-EST MARNE-LA-VALLÉE Stochastic representation of random positive-definite
More informationELEMENTS OF PROBABILITY THEORY
ELEMENTS OF PROBABILITY THEORY Elements of Probability Theory A collection of subsets of a set Ω is called a σ algebra if it contains Ω and is closed under the operations of taking complements and countable
More informationMultiscale stochastic preconditioners in non-intrusive spectral projection
Multiscale stochastic preconditioners in non-intrusive spectral projection Alen Alenxanderian a, Olivier P. Le Maître b, Habib N. Najm c, Mohamed Iskandarani d, Omar M. Knio a, a Department of Mechanical
More informationVector spaces. DS-GA 1013 / MATH-GA 2824 Optimization-based Data Analysis.
Vector spaces DS-GA 1013 / MATH-GA 2824 Optimization-based Data Analysis http://www.cims.nyu.edu/~cfgranda/pages/obda_fall17/index.html Carlos Fernandez-Granda Vector space Consists of: A set V A scalar
More informationSensitivity analysis using the Metamodel of Optimal Prognosis. Lectures. Thomas Most & Johannes Will
Lectures Sensitivity analysis using the Metamodel of Optimal Prognosis Thomas Most & Johannes Will presented at the Weimar Optimization and Stochastic Days 2011 Source: www.dynardo.de/en/library Sensitivity
More informationFunctional Analysis Review
Outline 9.520: Statistical Learning Theory and Applications February 8, 2010 Outline 1 2 3 4 Vector Space Outline A vector space is a set V with binary operations +: V V V and : R V V such that for all
More informationWrapped Gaussian processes: a short review and some new results
Wrapped Gaussian processes: a short review and some new results Giovanna Jona Lasinio 1, Gianluca Mastrantonio 2 and Alan Gelfand 3 1-Università Sapienza di Roma 2- Università RomaTRE 3- Duke University
More informationTD 1: Hilbert Spaces and Applications
Université Paris-Dauphine Functional Analysis and PDEs Master MMD-MA 2017/2018 Generalities TD 1: Hilbert Spaces and Applications Exercise 1 (Generalized Parallelogram law). Let (H,, ) be a Hilbert space.
More informationReduction of Random Variables in Structural Reliability Analysis
Reduction of Random Variables in Structural Reliability Analysis S. Adhikari and R. S. Langley Department of Engineering University of Cambridge Trumpington Street Cambridge CB2 1PZ (U.K.) February 21,
More informationA6523 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Spring 2011
A6523 Signal Modeling, Statistical Inference and Data Mining in Astrophysics Spring 2011 Reading Chapter 5 (continued) Lecture 8 Key points in probability CLT CLT examples Prior vs Likelihood Box & Tiao
More informationarxiv: v1 [stat.co] 24 Jul 2015
The generalization of Latin hypercube sampling Michael D. Shields and Jiaxin Zhang Department of Civil Engineering, Johns Hopkins University Corresponding Author: michael.shields@jhu.edu arxiv:1507.06716v1
More informationEfficient Solvers for Stochastic Finite Element Saddle Point Problems
Efficient Solvers for Stochastic Finite Element Saddle Point Problems Catherine E. Powell c.powell@manchester.ac.uk School of Mathematics University of Manchester, UK Efficient Solvers for Stochastic Finite
More information