N=2^q$ for a positive integer $q$, then the $\beta_j$’s can be chosen to be the design matrix for a saturated model of a $2^q$ factorial design in which the levels of the factors are set at $\pm1$ \[@BHH, Chapter 5\]. In addition, assume that $X_1,\ldots,X_m$ are random diagonal matrices of size $N$ and $Y_{jk}$, $j=1,\ldots,N; k=1,\ldots,m$ are random variables such that all the diagonal elements of the $X_j$’s and all the $Y_{jk}$’s are i.i.d. symmetric Bernoulli random variables. Then define $$\label{Uj}
U_j = \pmatrix{
Y_{j1}X_1
\cr
\vdots
\cr
Y_{jm}X_m}
\beta_j.$$ One can easily show that for any $Nm\times Nm$ matrix $M$, $E (\frac{1}{N}\sum_{j=1}^N U_j'MU_j ) = \operatorname{tr}(M)$. Thus, we can use this definition of the $U_j$’s in (\[ascore\]), and the resulting estimating equations are still unbiased.

This design is closely related to a class of designs introduced by @avron, who propose selecting the $U_j$’s as follows. Suppose $H$ is a Hadamard matrix, that is, an $n\times n$ orthogonal matrix with elements $\pm1$. @avron actually consider $H$ a multiple of a unitary matrix, but the special case $H$ Hadamard makes their proposal most similar to ours. Then, using simple random sampling (with replacement), they choose $N$ columns from this matrix and multiply this $n\times N$ matrix by an $n\times
n$ diagonal matrix with diagonal entries made up of independent symmetric Bernoulli random variables. The columns of this resulting matrix are the $U_j$’s. We are also multiplying a subset of the columns of a Hadamard matrix by a random diagonal matrix, but we do not select the columns by simple random sampling from some arbitrary Hadamard matrix.

The extra structure we impose yields beneficial results in terms of the variance of the randomized trace approximation, as the following calculations show.
N=2^q$ for a positive integer $ q$, then the $ \beta_j$ ’s can be chosen to be the design matrix for a saturated model of a $ 2^q$ factorial invention in which the degree of the factors are set at $ \pm1 $ \[@BHH, Chapter 5\ ]. In addition, simulate that $ X_1,\ldots, X_m$ are random diagonal matrices of size $ N$ and $ Y_{jk}$, $ j=1,\ldots, N; k=1,\ldots, m$ are random variable such that all the aslant elements of the $ X_j$ ’s and all the $ Y_{jk}$ ’s are i.i.d. symmetric Bernoulli random variables. Then define $ $ \label{Uj } 
 U_j = \pmatrix { 
 Y_{j1}X_1 
 \cr 
 \vdots 
 \cr 
 Y_{jm}X_m } 
 \beta_j.$$ One can easily show that for any $ Nm\times Nm$ matrix $ M$, $ east (\frac{1}{N}\sum_{j=1}^N U_j'MU_j) = \operatorname{tr}(M)$. Thus, we can use this definition of the $ U_j$ ’s in (\[ascore\ ]), and the result estimating equations are still unbiased. 

 This design is closely associate to a class of designs introduced by @avron, who nominate selecting the $ U_j$ ’s as follows. Suppose $ H$ is a Hadamard matrix, that is, an $ n\times n$ extraneous matrix with elements $ \pm1$. @avron actually regard $ H$ a multiple of a unitary matrix, but the special case $ H$ Hadamard make their proposal most similar to ours. Then, using simple random sample distribution (with surrogate), they choose $ N$ columns from this matrix and multiply this $ n\times N$ matrix by an $ n\times 
 n$ diagonal matrix with diagonal entries made up of independent symmetric Bernoulli random variables. The column of this resulting matrix are the $ U_j$ ’s. We are also reproduce a subset of the column of a Hadamard matrix by a random diagonal matrix, but we do not select the columns by childlike random sampling from some arbitrary Hadamard matrix. 

 The extra structure we impose yield beneficial results in terms of the variance of the randomized trace approximation, as the following calculation show.
N=2^q$ vor a positive integer $q$, then the $\beta_j$’s can bx choseh to be ghe design matrix for a satucatee modtj of a $2^q$ factorial design ij which rhe owvels of tis factovf ars set et $\pm1$ \[@BHH, Chaptgr 5\]. In addithon, assume thad $B_1,\lbots,X_m$ are random diagonal matrices jf size $N$ and $Y_{jk}$, $j=1,\ldotf,N; k=1,\kqots,j$ are random variables such that aml the viagonal elemenys of the $X_j$’s and all the $J_{jk}$’s are i.i.d. symmetric Bernoulli tzndjn variables. Ghen defint $$\nabel{Uj}
U_j = \pmatrix{
Y_{j1}X_1
\cr
\vdots
\cr
Y_{jm}X_m}
\beta_j.$$ Une ccn easily syoq tjdt for any $Im\timef Nm$ matrix $M$, $E (\frac{1}{T}\sum_{j=1}^N I_j'MU_j ) = \operatprneme{te}(M)$. Thus, we can use thms definition of the $U_j$’s in (\[avckre\]), and the resulribg esjimathng dwuagiohs ade stipl nnbiased.

Thia design is closely related to a bkass of desifns inerjduced by @avron, who propose selecting tve $H_j$’s as follows. Suppose $Y$ is a Hadamard matrid, that is, an $n\times n$ orthogonal matrix with elements $\pm1$. @aeron ecguaoln covwifer $H$ a multiple of a unitary matrix, but the flevisl case $H$ Hadaiard makes yhfit proposal mosj similcd fo ours. Then, using simple randim sampligg (woth replacement), they choose $N$ columns fgom rhis matrix and muptiply this $n\tikes N$ matrix by an $n\times
n$ dnagonam matrix wihh diagonzu entries made uo on itdependent symmetric Bernotlli randim vcriables. The columgs of this resulting matrix are thf $U_j$’s. Wa are also multiplying a subset of the columns of a Hadsmdrd matrix yy a rsndom diagonaj matrix, but wg do not felecg the colujns by vimple randjm sampling fslm some arbivrary Hadwmare marrix.

The dxtra structurr we impose yields veneficial results in ferms of the vaxnabce of the randpmixed tgacx appwmximation, as the foukowine calculatiivs snow.
N=2^q$ for a positive integer $q$, then can chosen to the design matrix a factorial design in the levels of factors are set at $\pm1$ \[@BHH, 5\]. In addition, assume that $X_1,\ldots,X_m$ are random diagonal matrices of size $N$ $Y_{jk}$, $j=1,\ldots,N; k=1,\ldots,m$ are random variables such that all the diagonal elements of $X_j$’s all $Y_{jk}$’s i.i.d. symmetric Bernoulli random variables. Then define $$\label{Uj} U_j = \pmatrix{ Y_{j1}X_1 \cr \vdots \cr Y_{jm}X_m} One can easily show that for any $Nm\times matrix $M$, $E (\frac{1}{N}\sum_{j=1}^N ) = \operatorname{tr}(M)$. Thus, we use definition of $U_j$’s (\[ascore\]), the resulting estimating are still unbiased. This design is closely related to a class of designs introduced by @avron, who selecting the follows. Suppose is Hadamard that is, an orthogonal matrix with elements $\pm1$. @avron a multiple of a unitary matrix, but the case $H$ makes their proposal most similar to Then, using simple random sampling (with replacement), they $N$ columns from this matrix and multiply this $n\times N$ matrix by an $n\times n$ with diagonal entries made of independent symmetric random The of resulting matrix the $U_j$’s. We are also multiplying a subset of the columns a Hadamard matrix by a random diagonal matrix, but we select columns by simple sampling from some arbitrary matrix. extra structure we impose results terms of randomized approximation, as the following show.
N=2^q$ for a positive integer $q$, theN the $\beta_j$’s Can be ChoSen To Be thE desIgn matrix for a sATuraTed model of a $2^q$ factorial dEsign In WHich THe LevelS of the fACtORS arE sEt At $\pM1$ \[@Bhh, CHapteR 5\]. In AdditioN, assume thaT $X_1,\lDoTs,X_m$ are randoM DiAgonal matrIceS of size $N$ and $Y_{Jk}$, $j=1,\Ldots,N; K=1,\lDotS,M$ are rAndOm varIables SUch thaT all the diAgONal eleMEnts of tHE $x_j$’S and All the $Y_{jk}$’s are i.i.d. sYMmETric Bernoulli rAndom vArIAbLES. ThEn dEfine $$\label{uj}
u_j = \pmaTRix{
Y_{j1}X_1
\cR
\VdOTS
\Cr
Y_{JM}X_m}
\beta_j.$$ One caN easily show THat For any $nm\TimES Nm$ matRix $M$, $E (\FrAC{1}{N}\sUm_{j=1}^N U_j'MU_j ) = \opEratOrname{tr}(M)$. thus, we CAn use thIS definiTion of The $u_j$’s In (\[asCOrE\]), aNd tHe REsuLTiNg eSTimAting equAtIoNs are StilL UNBIaseD.

ThIs deSign iS closely relatEd tO a clASs oF desiGns inTrodUcEd by @aVron, whO propOsE selecting the $U_j$’S as fOllows. SupPosE $H$ Is a haDamarD Matrix, ThaT is, An $n\timeS n$ orthoGOnaL mATRIx With elements $\pm1$. @avroN aCTUaLly consiDer $H$ a mULtIpLE of a unitArY maTrix, BUT the sPeciAL cAse $H$ HadaMard maKEs ThEir propOsAl most SiMilAr tO ours. tHen, uSing siMple randOm samPLing (with replacEMent), they choosE $n$ cOLUmNS froM thIs matrix and MultIPly tHis $n\TImEs N$ MAtrix By an $n\TiMEs
N$ Diagonal matrix with dIaGonal eNtrieS made up of indePendent symMETRic BernoUlli RAnDOm variables. The ColumNs of this reSUlting maTrix aRe the $U_j$’s. we are also MULtiplyinG a sUbsEt oF thE COlUmns of a HadamaRD MatrIx By a randOm dIagonal MatRix, But We dO nOt select tHe columnS bY sImPlE raNdom sAMpling frOm SomE aRbiTrary hAdamarD matrIx.

ThE eXtRA stRucture WE iMPOse yIeLdS benEfiCiAl resUlts IN teRms of thE variance Of tHE ranDoMiZed tracE approximatioN, aS the followInG caLculatIONs show.
N=2^q$ for a positive inte ger $q$, t hen t he$\b et a_j$ ’s c an be chosen t o  bethe design matrix fora sat ur a tedm od el of  a $2^q $  f a c tor ia ldes ig n  i n whi chthe lev els of the  fa ct ors are seta t$\pm1$ \[@ BHH , Chapter 5\ ].In add it ion ,  assu methat$X_1,\ l dots,X _m$ are r an d om dia g onal ma t r ic es o f size $N$ and $Y _ {j k }$, $j=1,\ldot s,N; k =1 , \l d o ts, m$are random  v ariab l es such th a t all the diagonalelements of the  $X_j$ ’s  an d  all t he $Y _{ j k}$ ’s are i.i. d. s ymmetricBernou l li rand o m varia bles.The n d efin e  $ $\ lab el { Uj} U_ j = \pm atrix{
Y _{ j1 }X_1\cr\ v d o ts
\ crY_{j m}X_m }
\beta_j.$$One  can eas ily s how t hatfo r any  $Nm\t imesNm $ matrix $M$, $ E (\ frac{1}{N }\s um _{j =1 }^N U _ j'MU_j  )= \ operato rname{t r }(M )$ . T hu s, we can use this  d e f in ition of  the $ U _j $’ s  in (\[a sc ore \]), a nd th e re s ul ting est imatin g  e qu ationsar e stil lunb ias ed.
T hisdesign  is clos ely r e lated to a cla s s of designsi nt r o du c ed b y @ avron, whoprop o se s elec t in g t h e $U_ j$’sas fo l lows. Suppose $H$ i sa Hada mardmatrix, thatis, an $n\ t i m es n$ or thog o na l  matrix with e lemen ts $\pm1$. @avron a ctual ly consi der $H$ a m ultipleofa u nit ary m at rix, but thes p ecia lcase $H $ H adamard  ma kes  th eir  p roposal m ost simi la rto  o urs . The n , usingsi mpl eran dom s a mpling  (wit h re pl ac e men t), the y  c h o ose$N $colu mns  f rom t hism atr ix andmultiplythi s  $n\ ti me s N$ ma trix by an $n \t imes
n$ di ag ona l matr i x  with di agonal entries made upo f indep end ent s ymme tric Bern oul li ran dom variab les. T he co lu mns o f thi s re sul ti ng matrixa r e t he $U _j $’s.  We are  also multiplyinga  su bset of the c olu mnso f  a  Ha d am a rdma t rix b y a random diag onal matri x, bu t we do no t  se le ct thecolumns  by s i mple ra ndom samp ling from  s omea r bit rary Hadam ard matr ix.

Thee xtras tr uctur e w e impo se  yi eldsbenefi c ial  resu lts in  t erms o f the  v arianceof the randomized trace  appro ximat ion , as thefol l owi ng calcul atio ns show.
N=2^q$ for_a positive_integer $q$, then the_$\beta_j$’s can_be_chosen to_be_the design matrix_for a saturated_model of a $2^q$_factorial design in_which_the levels of the factors are set at $\pm1$ \[@BHH, Chapter 5\]. In addition,_assume_that $X_1,\ldots,X_m$_are_random_diagonal matrices of size $N$_and $Y_{jk}$, $j=1,\ldots,N; k=1,\ldots,m$ are_random variables_such that all the diagonal elements of the_$X_j$’s_and all the_$Y_{jk}$’s are i.i.d. symmetric Bernoulli random variables. Then define_$$\label{Uj}
U_j = \pmatrix{
Y_{j1}X_1
\cr
\vdots
\cr
Y_{jm}X_m}
\beta_j.$$ One can easily_show that for_any_$Nm\times_Nm$ matrix $M$, $E_(\frac{1}{N}\sum_{j=1}^N U_j'MU_j ) = \operatorname{tr}(M)$. Thus,_we can use this definition of_the $U_j$’s in (\[ascore\]), and the resulting_estimating equations are still unbiased.

This design_is closely related to a_class of_designs introduced by @avron, who_propose selecting the_$U_j$’s as_follows. Suppose $H$_is a Hadamard matrix, that is,_an $n\times n$_orthogonal matrix with elements $\pm1$. @avron_actually_consider $H$ a_multiple_of_a unitary_matrix, but the_special_case $H$_Hadamard_makes their proposal most similar to_ours._Then, using simple random sampling (with replacement),_they choose $N$ columns_from_this matrix and multiply_this $n\times N$ matrix by_an $n\times
n$ diagonal matrix with diagonal_entries made_up of_independent symmetric Bernoulli random variables. The columns of this resulting matrix_are the $U_j$’s. We are also_multiplying a subset of_the columns_of_a Hadamard matrix_by_a random_diagonal matrix, but we do not select_the columns_by simple random sampling from some_arbitrary Hadamard matrix.

The extra_structure_we impose yields beneficial results in_terms of the variance of the_randomized trace approximation, as the_following_calculations_show.