From 00c0aea10ddc9c7b510e51af5aa68a639873202f Mon Sep 17 00:00:00 2001
From: Viktoria Petrova <vipet103@hhu.de>
Date: Mon, 25 Nov 2024 15:00:39 +0100
Subject: [PATCH] add statistical analysis assay and field experiment protocol

---
 assays/StatisticalAnalysis/README.md          |   0
 assays/StatisticalAnalysis/dataset/.gitkeep   |   0
 assays/StatisticalAnalysis/isa.assay.xlsx     | Bin 0 -> 6900 bytes
 assays/StatisticalAnalysis/protocols/.gitkeep |   0
 .../protocols/FieldExperiment.md              |  19 ++++++++++++++++++
 5 files changed, 19 insertions(+)
 create mode 100644 assays/StatisticalAnalysis/README.md
 create mode 100644 assays/StatisticalAnalysis/dataset/.gitkeep
 create mode 100644 assays/StatisticalAnalysis/isa.assay.xlsx
 create mode 100644 assays/StatisticalAnalysis/protocols/.gitkeep
 create mode 100644 assays/StatisticalAnalysis/protocols/FieldExperiment.md

diff --git a/assays/StatisticalAnalysis/README.md b/assays/StatisticalAnalysis/README.md
new file mode 100644
index 0000000..e69de29
diff --git a/assays/StatisticalAnalysis/dataset/.gitkeep b/assays/StatisticalAnalysis/dataset/.gitkeep
new file mode 100644
index 0000000..e69de29
diff --git a/assays/StatisticalAnalysis/isa.assay.xlsx b/assays/StatisticalAnalysis/isa.assay.xlsx
new file mode 100644
index 0000000000000000000000000000000000000000..3d433531242932d3bd5d748c8c8ca306e897d1af
GIT binary patch
literal 6900
zcmai3by(DE(+25YdZk+$K}tzMbXmGva_MF1kd*GFM3F{9kdl-X38lMHk&s4uY5f*?
z(Q`cK{odz~z4qFjd*-?BnVp&6y{Zam=(s2-DA*{q_Fj7NGc7oJs3<7g=qM<p$WIMq
zVGgbk2UlZFFGq-r5w9oIz9dP>q4g#p<WNFxkq3-#q^Tv>np4$7QMivX;wi&TGqd~b
z?TiNrdl4r?kqIeRipZ9CM*DS@ni%<UAV@4VN==vyXL^_eN?CW5_BJ2Rl#Q+Z*ongB
zwJb^Eyo30xS8i}B&)1yTAt`xzNDP7MlDK}UCa)miSeQFE%oztS>9B1sh&w&T+MLRO
zNG-ZR_g#1UBBTd?c=s9=N!Q-P$qHxG&(uaTi}i*E2Vn_x*w;PWsfwx#Q1ODr<*m7v
zQNx-&G86}K6R!s!VG_Jr9Pf<5>)~5Yo`NKjau{i?l04zkA*i)jL<`KJ;umIHA1_9Q
zVx6BUD^?iX7a?@s-|8Yh_Fee?$Tzq@9}DDU73Uy<=&yns%1R#)qtQpFVjXSz&LtWd
z%YXfD6F{KR4^7Ie{35<XOmOtW=zH3`8nFZWV8*qc7ZM~>Cu^_G#JED%&<)(Bisz|S
z|C0%npP4|C|Lb5P-x)hY>|OY-Y%a64n<+oNiu_LD8VU;DFB{$q%aVk5%569CrMG<K
z7Ehl$fUvXk8W?gvu?b4HusBdx$Gm6{`nDAjCygItgB9zvv24||-X|Wj7Nd@!!RLNO
zMTXJSAJ*T}^WDaYh!YPpj2TP7F&95|$?{y08C=-a1lJ=6bz(k>Kgx2ekx+WZ6-(Bw
z(qkjTF)+Fs?%oo+UsSr~deRlJ6kJ%-FFjq)`tqEASjgVN0Ja*xX_~`erGM7VZlOpU
z@z`s0D=k!KDVtMuaM?(EqTx)UAOkd*omG=&O{`I<Jl~`L){SPy2xoa9T0y4vM;*=c
zkU(Y+Jsk5hf#Dk>gFYx%jQo=|G$bRQ_P;3mrvc8T!M`yZ%<@|Mg4sGGvs8aI^ME<q
znZaOo7r%rnqHWu50`6^{<Gw!Bch?7V#Jg}jHZ(Xy3$ChBw3oLQfu0{~QG;l(v#~PU
zizib)fY$pCMZ;ATN3qJXIKhAwrShgNG9SU2Q!7CGw1EUV9G8t{a@1p8ghp_g*4dXO
zUP+w8v3ExfEI{}`NiKY4Pl<Z?<?Wu)u9tleJ=k?1$1Ecs(bn)(URUA{_iHme=a~0z
zl$}9f8VB8uUO5H_vt)$!@=>MVvA}p>WI-P{v7n>y-sN+L%Som<;VT1tqA2Yzt}>_n
z98%kj5ywXqy$cC9jTtw)7<aEP3ZBUIeJH%{wX?8B|KEuwyd>y?8W(E_#P!cIt&LQx
z4@MH%iX<HP58Iy*{tG^zB&)$vZvp~l54qxIsYOLt7_Sw%t5wUV3PRHn8$i<HYxips
zYfcXuL=_fGv&M+xMX4nTPn*;BQkP3(2xX`XPXb%G8RJ*#iY(r8bq<7`z9Nr^hWK{q
zk3|JGV#_*o6CEcG!}Ez-?_;kh-KDY%bEBdD93GpJLUj**t9qEW61Zk$kIAF^m2PeN
z8Cw_Aw`Ut&5#{4srZ4N04eiNhJ#uPqgj4WTK5MpFH0iI3eTy~0fBuHP^vxRPJ$ne(
zW>iqZ%Z2rUhKDva(4ERm)BdBB%|MYOW#^LUqD_1!e^TGmKE1RK5a%9vo{~mKVmODl
zlc#u#J?>5Kr@lkpF4&!-dwqzK^U@4SoH{%1Z)#ws>e`Zx4oW8gk1{!nmjpts+MyRx
zU;0qE)zPplFn9GUIklx;iroR(TsPXc-Imx>Y&)PAZ^{r^;?f6i*bXL*0Fb_$1x^H^
zb2(#NYkk2ezP4kM=0V%qfbeup6Cs-GllXif%11DlM+cMaHwfOQpyQ)m=8^SOfY$<-
zOJqE`;I%Z%5IXpG_*1XOv5@7ZZ)bJ(XGU@8<lgf-%lGPJpGh@K-&ETH<N_L<t#eiw
zJDIkH0we0f7~)$*aM-}~fri(N*`B@h!yd-*$uyBc!4|;v$unt1#ZJKWd1WH>xHmAa
zS^TCc5kW=R31YC&lGuXiU^Dy2G^uwJ)I@CO$!hwB(rD!LV@KxwLU_I1D@>0NttzeY
zqHi=plQkid_g~&ENs_?Yi8z&=JO^;_@Za?~VeibU8dKEw0e<l3O0yOTYPXcF^s3Pv
z#;@}ntw?|?0%^!TNZl}etv`)*prBpjls*;m6ye{?&`k5TTk5LHFB>fMJrv6fISz%9
z<B;sKa$T%Vogo$)uFf_NRxUrQHc>?-NE*<w<rmrad2^gDnqOKLpM|YPSe1zG!Lwu~
zgvuv9#bX3-3_oXZXMD=+jBn*$__$TcB-WjT@%A&z*NgZ{5IdV#e*dZE36HdRERCQF
zO!Lo-J8sV8#9kD+p=O^w{p%?Mi!g--TQdfO``_Imko2j4k(DU>vheP{e#Drw7ndNY
zO%GzQVBW0nG){H=dA`kqF#^OEuQFVrXv{p92-*)RV<L9750WmCjd@?`>(Tl{%^t2?
zSydEwl08*323uLOR+y>B)r<eSP$Z!p^SYa^0mr~U>vrk=yjzN!WJi7y3nske6+A1C
zVJkB%V4{Z!Psv}72CRV+;?ScMyq&aN>C$G(;x;Ro{CS1-)M95F@EGxfr0#;oClB7f
zG_3z&LakT%fVl}>F42Ty_a$3|5$UefCQfMLkBtJi%~tn<l(xH^*Xbso%a@kp7O2z%
zT-1d|W8}cW^is#@S2_8w9o*F#0)_l)^p~ki=oh*XjU1fq440b>auQsc{e1n|Z#MN*
zU^BvGKR5%wm0B3=W@b5<_PCD-+XTMf9Asu5p)dyvf;c18$bi!$yl9xNZM<miGp0ao
z6AbN*3L?_1j&a{`3Hk5L)0(L@HbbUPD5<%ATRR^P>FAByUrngU+R~WwVn^v)9!El{
zM0cm+?kyyH!3Oz+6lFnV1jG<xhm>2~z3UUDIXu_D0H_&gBX~vy&^QVVl39I#_vdj_
z=dm+lfJE`Zyz~aLvC7hRppr=FL@^+~Q;>-)aJXqG{Ad!$!r7O%OBYfGinM4cIBQX{
z>Q%^{oXbr+3s2c*7l|c?rraEeOUK_X6eNfY4#wgcQI^Vg+bpUk%c+nSzA>R9<x*GI
z;;1hF_2|fvV@0%7!;bGjN8mj5&BG?0^|9m=b(CFST^qd@l@E8)ZVO{|Kd1Wi=~1n|
z_V(@bF#?|%rSD`uZXq7R+vm--6{dtqLZ1j4XKdW-i~Kx7%AjA?-5b_N$$fBf*(2Df
zt!GNKnmTNs#u2}m(f1s2J2~j4X_nEeTj^e8)+4iHykpi#6fjx`O6ZKM<e+1uccjZ%
zj3YG<i$=q4<cNG0VA^1IxXI~r!(VwCED>3Avs>Ouws{f=J_5EbTW}gI@eJB7G4L5T
zih7@&ZR~7Wsbr){e0P&P%>O7H`QEE;Ldth<fXJB2{|91Wy<m22hY9X~(4$lK{b)~V
zEl*~>!Ky4n%KP|y!vFkC6Ak9#`YenH6dvM)l?n9F!^Rs`3?MpgkjKAvj|_(lMICLU
zOB(fa3S&{5Te)rPDUCn~p(nGX`WrfQ@#K4GjkLkIj!*fgO${5;X1WZIzqSwQO5KL<
zy=lOvU0*_s(|n970PNl?7{)7?amfs<9RfbkrRusR{7#n{5di2YZD?OsERPqY1H1?d
zA7?>xW<QXbVL^Y}>x%wvEE4LzMNNC(ZLQ%*SFW{2TC;#eqa#{Gm0Rs4BwFVAgKrGi
zS3q>VZ_0p%pzVF2hUI>!E_Qk^T_%SsOg40=!S`vu40#u8%1WGA_-pZs9=3-8t%b?A
zA_Bz``Ol=onanfgd2sT-#pM$DNvmZyY>{%$LJekv@N<dc7Ug^046faBlj>i=Hcj6Z
zCVH|ly>yS7Hr=y07(w`+EIo_y{rjfWJL0JV5~~KN7EhqXm3=(9?`?<S<pBa5g?56m
ztTH%LXyTR_KJHd_5muo&fE}hj4>sG*{Vz!)`NI-xDoL7P&kCCZ7R0^q{yaA+X|OJ<
zZ#I15I{dR6{BiG)cjp(6<j3uP4m0?Gj+21v4D2^uS1In9HGwM9)nDs!iKvTgb`y;2
zH<2c+lynT}Pc8U8uULLTE>3B@MwXX!aIFCQke=JFo;?=9pT_c3s}G{J=U;3aiP(kk
zxajp5jtgWKC+j`o)Bet*QpC#b^Ys=}5?$x@4lyPkAM}lLQ+mu1WzH&}HbH`{+nLJm
zVs3+HFs*j+n1X^(luVyH-)YnB)q;ocL#}~uRH}ODKkAs{rM7(Wb$33G4V$5kDs0$}
zpD^Tk?|r+2v6+UvlPVDkuWJrP`)4f9iG?ghR3ry$$DpLmD6^r_a_gko_SSAM48j{|
ztgtqt8|LV9n2ZvNXGvXrjoxz^yaM_48<Cb?O3SU{C2^f0sbD6`P*fVa^)||lbxVz{
zwaDcAO*{NPQ0NjSRiVTJI>K-aLrC&@O)v4;2zs9*?#aD`LICXpo{^_6vTs!^Q=Aza
zK<L-GeUnP0SV&xg5&&irTbNWOZ=P_uo$j$tR7pwo=&dWie?g3Zyz4B&8xUZfE}+L&
zzQc_<jP-82T{M^fZB&b@LsJ^i5fp5w?*j7VeWRmosGze~^3t%TxQ0Dv1f@?U$wy$|
zi;cJZjoR3x)srIbbR41&=TV|=v%FQwf_Q{+k4z~}J;QurwX_Dh$Cm0X+x+%8cIL2T
zV)bQJN{R}!UEK{$?#88-YMa!8vqhw~1$E*2Tsk`WO)&;gcTe>A*$`Wye02-H56anR
z2#UjaJvPt-%BioKRQxx1_Odr?Xv^jwu{O*hz{+E<6$omLc&d$R^}p3Y8fFT`uf6DG
z;BJ`u=zC<9*}wD@v3m5r!}kq)b=eVjm99RPL80R#C*x8w<B!zn+<WEPt0J%ZluJg7
zYff1UW^9t4KP=H$8*hA6JA3f)e&vr&!B~y*k<v%XB@^sp-5EO0w(c@-hq%Ym!Fijb
zGcg$>euA|^#XdU?H(}lD({Hb4+~x4X5DR$dj@+2hkOPe1KYNF(mp$a?!jYF`gj|!!
zLbndNxMwxc2)H8xY;7Vn3sX6T65ZDW<cm+xu8+hkoo<O?l90Gf=Oem5l)*FiPHxq$
z=1TMR3L3s5!_twwuSj0MX!~@&EcIim`KWFYn6t`^L$Vxm)`C9aJsrCD^eszbfM()+
zQ}sbYT>j)O`BW~WTNunwgcI~vp3~acwpmZTkyno!JTzzG(D7v8hU`cI<ojF<yhD<O
zs(pGeYQr)@;<=L##qtO~;W|!Z+|tkD4Oj9k9#klc1U&?PTIqg>nzT|hkKw)Pc9$ZU
zf)uA5#*9I{(`GuSr)~qG)?CrnoAk{R3+u7Ud;97m!EL9`?gSY;Xk3VTZx^~NXPsyk
zC8s*goota-%)8E$kq^)KvF}Pf7&33=A@{L+L)xX94Ws(v%&1%{BaW?ijoKka?Y8S6
z5Jo8UOwCaUQ!(=pzLB0C#t_}lwofqn$yTFPD-E&04v8>%5`q)Dg($7AXLA_TpIqVQ
zxc$jSE<~2jDjL-UkF7nip=`DY&s~RctlMz<RG`Cg=hGUT5lGdBddrpON1AZ?Y>b?d
zTbES!sSM?QV&q5p{G%vuYBQNEXiv!3HKNeRdhUr6exvmtob@n(tgo!OHkJgR&2xP3
z1%K||PCDo`uA;5XKYEw0322}yt|W@tG#I7In9uZh1&=x~qI=}xFLY=*!M#f{WV<O>
zIWzVso$qY$?YSGmc6P_|Ll|ZpgqfiyI*PVWN=ZL*;7-uF-;WkV-Ju-tG&85=aMw<1
z+mb)sd++GZ?k+p_#E!ii<=3H<!<Cq~!770!0`xpQPZb?IcsF7FTB0l~-yy_{o+K70
z&6TrP1^Vwxm<7yS)fwja=aXO8*EUEna(`w+K8pj8b^OcB)Y0)`cTP{Z3!Ax#AF_2y
z9{MdN&4E2sTAM$Y`jPHy6#1EjmM&Lr&1Cd=p@j33F2eac>&guh2}A>#97Qv-%aL@X
zQ(w4V8z$vQ!(^c6%_Q&!m`<kLUHoQDN^$M^a2AIHtZI2iGqWhoMRb)X)+q~MumF#{
zbd)&4YdClmby;8&VohbY&HHSEP@-qOqA}jnWM61p#|f9Lpj~U=2em%dt?i+aZ$Os}
zZDTXnag^oEQKIo&>Fl}r^><ks&XSt9t)rHN^d#5RySxf(FlkC%${kvWIgiIx53$?&
zg-u9p#`cOHv@w!f_%D|jMa@$J2t}W)FzRh{E7xf7+fqg(`Vp>%=-*>wEm&e-{#;j+
zx$i7!;bfCU3L1YJF4jZ1`k?Z0`wxbP{Sgl^`OT})r5<C)ZUF}9Pr3YP>#-%bRF&6I
z#W$R`mDd^$`1XMwbA-SFg=YIAC?3xG8Se<QB7s|`O(ON#I>i~?dgb!DiD*Q@$kh7R
zzMM1N<OS9=Z*a5Go9}2>iM{MR2fF)7#>k8<AiEFoUo&P7bB6p(nR?HK{MUHM#WWiv
zpOqIW6E_{H5zqOgNDEflg#NU>W+RW$-+1{h$`cKkm{{e4uc#4CH`5X7dakw$Nz8i|
zl#c2+&dkKx1A)dHk7e`GxvBzh7ru*^b5C?mbcD`z*rHLezG<Q6c<OlThQYl=)wLrE
z@Y4N<+2VCWH5Mz`eTAkzADpvxJbfN0NV5S2rh%cBjx0@iUxM<cW^#l8X)>s+?5=`T
za{IY*<#`SIg929F>SoX{Mf)JawV=eH@}eR&A)CDBAK`w*8OAeevW3sJoEY>tJzIHh
zjtYBgh~;rT8Ll*Jd(nQaLikgdr1NPR?8{^x))#ZsSOc5E*+OPMlJ7KPFS(#0dm9+G
zyd3~%=z8oOMk3rL$FBoT@eFip8fE^4RkKhLwD6WKmQ`U-i_bUUD(kFW^G@OfFW94D
z6GJR1y$P!um8-`&-z{_z@Uxj6>%jq!9llEAzh}deN7l~}nGIKDHi$0gz+VNZm;@yi
zq{tV#br#8fSi%#8IyNnj9p`?RCV*>K-6JZQuVYpBhlIm}cmO4@(?&boSWK@(Wjvhu
z{_tn5if7B}SsKkl4mw@669rQxWJxVTCXH=i1jyOQxLDFx+`dT)yWc~z%NLy8E1>RO
zCGA2#xFf&*BBb5nBZQ1eM%KICLdiX4xu-YKeOHlv!5`QRt!~9CyA#Em&PYtaqOIah
zYn`KaK%W-kxzeiWMGbsCcE7+c7;LIlH+7rd)b7o?hLc3R*C`gXIx@p2bt(FU;<#Xz
zX&`d`qjQsO^}EBpQH<bieD{}Txi_U(Uk~A=f}NEXq;_-|+~=A{Lem@Mn&b4&_el8E
z1Se&zhPhrnUpdR~dDY0C`y~NqL1FZPszg}WM>Ses1OC!fOHon>&0HMG0PK?;QkGa6
zAsadzJPx7R=+q@!6bPyOW3^(7o6M{EoqVwO)QSBkc&ZAh*T_)*yE2Z<#$V3A=;Xg0
z{Ig|V2)dW&qANxsGkhgN{ypfQ0`C<H1*P;C=BEnz{{_3yK`+lmuoh$~{ZvJ-BK`Ar
z7rN->xd6;v0{;Dif9j-HK;-`k$gs#bm$1JU?f<;*g`9YKF2Zvlr5B{1zln=iSHEmo
zf69wjXcQE!e_(%$jlUVY^zx_F^V<s%*3TRKL-_ga;eS#10{f?jE1BZ&*Z)(f_#L?W
zudv^;#qYTPzAXNB5Qgj-NbX#It&YE;|Le-^B!8LKzhRdV{~YVTBdYur_MajDC;HD3
w^Bawd)YyKl2v-B=Pxzlb@i!b!{;O23Bzjc^jGw<X0CN3!kE}=_<%RkG0Z#93O#lD@

literal 0
HcmV?d00001

diff --git a/assays/StatisticalAnalysis/protocols/.gitkeep b/assays/StatisticalAnalysis/protocols/.gitkeep
new file mode 100644
index 0000000..e69de29
diff --git a/assays/StatisticalAnalysis/protocols/FieldExperiment.md b/assays/StatisticalAnalysis/protocols/FieldExperiment.md
new file mode 100644
index 0000000..43010df
--- /dev/null
+++ b/assays/StatisticalAnalysis/protocols/FieldExperiment.md
@@ -0,0 +1,19 @@
+## Field experiment
+
+Due to the strong dependence of photosynthesis on light intensity (Ogren, 1993), we considered three light intensity clusters when analysing field measurements: LL, ML, and HL conditions. These light intensity clusters were identified by K-means clustering of PAR and LEF. In addition, we also compared three main developmental phases of barley, i.e. slow expansion phase (SEP) (ZS<30), rapid expansion phase (REP) (30≤ZS<60), and anthesis and senescence phase (ASP) (ZS≥60). These two factors, light intensity (L) and developmental phase (S), each with three levels, were considered when analysing the MultispeQ parameters from the field experiments based on the following linear model with the quantitative covariates light intensity (PAR) and developmental stage (ZS):
+**(5)**
+$$y_(p)ijklmnopqr = \mu + G_i + E_j + (G:E)_ij + ZS_k + L_l + S_m + (G:L)_il + (G:S)_im + M_n + D_o + PAR_ijklmnor + T_ijklmnor + (E:R)_jp + (E:R:B)_jpq + \epsilon_(p)ijklmnopqr$$
+
+where $y_(p)ijklmnopqr$ was the observed MultispeQ parameter across all light conditions and all developmental stages, μ the general mean, $G_i$ the effect of the *i*th inbred, $E_j$ the effect of the *j*th environment, $(G:E)_ij$ the interaction between the *i*th inbred and the *j*th environment, $ZS_k$ the effect of the *k*th Zadoks score of barley development, $L_l$ the effect of the *l*th light intensity cluster, $S_m$ the effect of the *m*th barely developmental phase, $(G:L)_il$ the interaction between the *i*th inbred and the *l*th light intensity cluster, $(G:S)_im$ the interaction between the *i*th inbred and the *m*th barley developmental phase, $M_n$ the effect of the *n*th MultispeQ device, $D_o$ the effect of measurement date, $(E:R)_jp$ the effect of the *p*th replicate nested within the *j*th environment, $(E:R:B)_jpq$ the effect of the *q*th block nested within the *p*th replicate in the *j*th environment, $PAR_ijklmnor$ the light intensity of each measurement, $T_ijklmnor$ the ambient temperature of each measurement, and $\epsilon_(p)ijklmnopqr$ the random error.
+
+To estimate adjusted entry means for MultispeQ parameters of all inbreds, $G_i, E_j, (G:E)_ij, ZS_k, L_l, S_m, (G:L)_il, and (G:S)_im$ were treated as fixed effects, and $M_n, D_o, (E:R)_jp, (E:R:B)_jpq$ as random effects, $PAR_ijklmnor and T_ijklmnor$ were covariates. Furthermore, we calculated adjusted entry means for all inbreds for each light intensity cluster as well as each developmental phase.
+
+In addition, to evaluate the effect of each fixed factor and covariate, analysis of variance (ANOVA) was conducted.
+
+To assess the heritability of each photosynthesis-related parameter at each developmental stage, which was considerably shorter than the above-mentioned three developmental phases, data were separated into eight stages from Zadoks principal growth stages. The adjusted entry means were calculated based on the following model:
+
+*(6)*
+$$y_(pd)ijlnopqr = \mu + G_i + E_j + M_n + D_o + PAR_ijklmnor + T_ijklmnor + (E:R)_jp + (E:R:B)_jpq + \epsilon_(p)ijklmnopqr$$
+where, $y_(pd)ijlnopqr$ was the photosynthesis-related parameter for each developmental stage across all other factors. Due to convergence problems, the interaction between $G_i and E_j$ was removed from this model.
+
+To assess the similarities among the barley genotypes with respect to their photosynthesis parameters, we performed hierarchical clustering by Ward’s minimum variance theory (Ward, 1963) using the adjusted entry means of PSII parameters and SPAD at three different developmental phases. Furthermore, principal component analysis (PCA) was conducted by using the adjusted entry means calculated for each inbred in each of the developmental phases described before. In addition, a PCA was conducted by using the environmental factors (temperature and precipitation) of the inbreds at the country of origin. The relationship between photosynthesis-related parameters and morphological or growth-related parameters of the inbreds was evaluated by Pearson’s correlation coefficient among adjusted entry means.
\ No newline at end of file
-- 
GitLab