4Misc_Start4Platform@ xHH@Rg(HHdh xHH@Rg(HHdh x HH@Rg(HHdh ^Graph*WDashSettings#  ! ,o6Normal@ Lucida Console<HHHH$$4 4 4 4 4 4 vhomeB:dC:Users:LaserMaven:Desktop:C:Users:LaserMaven:Desktop{k\LaserMaۨ\[YC:\Users\LaserMaven\C:\UsersDesktoperMan\HC:\Users\LaserMan\H\RecentWindows(Graph0:wVarFactor,wVarFactor,wVarFactor 4Misc_EndTXOPState_Start PeakFunctions2-64rMaven\4XOPState_End\•GraphWave() !y_ u ֫-g-g,wVarFactord????@prCs?]䷽@o>@L(D @Zw'@Vy.@yf 2@$z|4@)6@͗8@1A-:@N\;@ަMK<@6%=@GJ՟=@U\$>@dR>@%א>@1땝>@nD,q>@Xˊ>@*s_ ?@30?@E-%?@Mp,?@y51?@S-5?@!~DB8?@v;:?@-ݥ;?@0@?@uV>?@*>?@L3>?@Ѓ>?@+??@0шo ??@;G@,*VN@GR@N#U@f%2@Y@lbA\@R#3_@҅.3Qa@y崪b@>d@t::e@ 0Uf@YSQAWg@bBh@ i@2$V'i@^P~j@R5k@S;Ük@.2l@u}l@FS;l@D+m@Zsm@͑[m@Ϛnm@O([n@S'@n@ޝNdn@rDn@= n@+n@|n@_sIn@vUn@jzn@!o@[ o@Eo@!T\ao@]o@Ԋ㳪#o@7a'o@++o@"A,.o@II1o@3m3o@њim5o@4  7o@|4k8o@9o@y译:o@ܓg;o@;o@1nP>o@!x>o@"/>o@=->o@{ ?o@vW*B?o@rG^?o@7u?o@N߉?o@|?o@ ?o@ є?o@q?o@FUx?o@Atx?o@^?o@,?o@ Y?o@ןm?o@%fY?o@B %?o@=?o@L_?o@aA\?o@ ^9?o@%^B@^d?L@E4T@:Ge[@q4wa@8~if@Ąrk@;op@.(Ps@|Uv@@La@DMG@.y{@41@Sڝ@o>D@9 @ g"@ 6@|0RI@3BǼ[@Z_l@c~%{@5}A@r @:È@7ʦ@e d@52\Ğ@* ͞@s4֞@잵+ޞ@@h@@($@ldw@ȝ@m@-5, @p@1@ g;@!˻e@`jٙQ@feZ @a+"@G$@^]P&@u~(@HƼ*@X}-g-gwLambda???????*// Platform=WindowsNT, IGORVersion=7.050, architecture=Intel, systemTextEncoding="Windows-1252", historyTextEncoding="UTF-8", procwinTextEncoding="Windows-1252", recreationTextEncoding="UTF-8", build=29747 #pragma TextEncoding = "UTF-8" Silent 101 // use | as bitwise or -- not comment. DefaultFont "Arial" MoveWindow/P 681.75,44,1218,366.5 Graph0() MoveWindow/C 678.75,410.75,1216.5,491 KillStrings/Z root:gWMSetNextTextFilesTextEncoding Window Graph0() : Graph PauseUpdate; Silent 1 // building window... Display /W=(35.25,41.75,666,490.25) wVarFactor[*][0],wVarFactor[*][1],wVarFactor[*][2] ModifyGraph lSize=1.5 ModifyGraph rgb(wVarFactor)=(0,0,0),rgb(wVarFactor#1)=(1,12815,52428) ModifyGraph grid=2 ModifyGraph log(left)=1 ModifyGraph tick=2 ModifyGraph mirror=1 ModifyGraph minor=1 ModifyGraph lblMargin(left)=12 ModifyGraph standoff=0 ModifyGraph gridStyle=2 ModifyGraph lblLatPos(left)=2 Label left "Inverse Variance Reduction" Label bottom "Number of Fitting Points" SetAxis/A/N=1 left Legend/C/N=text0/J/A=MC/X=-22.27/Y=-19.23 " λ\r\\s(wVarFactor#2) \\{wLambda[2]}\r\\s(wVarFactor#1) \\{wLambda[1]}\r\n\\s(wVarFactor) \\{wLambda[0]}" EndMacro #pragma TextEncoding = "Windows-1252" #pragma rtGlobals=3 // Use modern global access method and strict wave access. // Cramer-Rao Lower Bound on exponential slope variance; from derivative of // log likelihood function. function fvar(n, lambda) // denominator in: Variance >= sigma^2/Denominator variable n, lambda // lambda is exponential slope factor; exp(-lambda) // per data point (data normalized to 1 at x (or t)=0 variable rho = exp(-2*lambda) variable vsum vsum = (-n^2+2*n-1)*rho^(n+2) vsum+= (2*n^2-2*n-1)*rho^(n+1) vsum+= -n^2*rho^n + rho^2 + rho vsum/= (1-rho)^3 // finite sum formula from Gradsheyn&Ryzhik 0.114 return vsum end function GraphWave() // 2D wave for Graph0 plot make/O/D/N=(100,3) wVarFactor make/O/D/N=3 wLambda = {0.2, 0.1, 0.05} setscale/P x, 2, 1,"", wVarFactor wVarFactor = fvar(x, wLambda[q]) end