search for: x1i

Displaying 9 results from an estimated 9 matches for "x1i".

Did you mean: l1i
2010 Dec 10
0
locfit weights not working as expected
...d it to. I would gratefully appreciate any advice or help that anyone can give on this issue and I appreciate your time very much. library(MASS) library(splines) library(sm) library(quantreg) library(locfit) set.seed(20) Zi = rnorm(5000,0, 2) ei = mvrnorm(5000, c(.6,2), matrix(c(.7,0,0,.9),2,2)) X1i = exp(-.5*Zi + ei[,1]) X2i = exp(-.5*(Zi*log(T1i) -log(T1i) + Zi) + ei[,2]) X1i = log(X1i) boot.weight = rmultinom(1,5000, prob=rep(1/5000, 5000)) boot.subset = rep(1:5000, boot.weight) loc.model.1 = locfit(1*(X2i[boot.subset]< 10) ~ lp(X1i[boot.subset], Zi[boot.subset], style = c ? ? ?("...
2008 Jul 17
1
combining lists of pairs
Hi everybody, This has been causing me some trouble: I want to combine several lists of pairs into a single?vector but keep the pairs together. The lists are of the form: (1)?x1 x2 x3... N? (2) y1 y2 y3..? N I would like to keep it this way, simply appending one list to another, e.g. (1) x1i x2i... x1j x2j... N (2) etc. A?method like c(...) does not achieve this. It creates a single list, disregarding the pairs. Is there a more efficient way of doing it? Thanks! Oliver Marshall [[alternative HTML version deleted]]
2012 Sep 27
1
What to use for ti in back-transforming summary statistics from F-T double square-root transformation in 'metafor'
Hi Dr. Viechtbauer, I'm doing meta-analysis using your package 'metafor'. I used the 'IRFT' to transform the incident rate. But when I tried to back-transform the summary estimates from function rma, I don't know what's the appropriate ti to feed in function transf.iirft. I searched and found your post about using harmonic mean for ni to back-transform the double
2010 Dec 26
0
GLS with corAR(1) correlation structure residual/standard error calculation
...approach does not yield the reported results. Is anyone familiar with the inner workings of this function and can either explain the calculation of the standard errors or provide code that explains the inner workings of this function. Thanks! Example of the model I am running: model1<- gls(Y~ X1I + X2 + X3 + X4, data=Dat1, correlation = corAR1(), method = "ML") My understanding of model errors: Y = b_0 + X1 b_1+ ...Xk b_k + Z Z_t =phi Z_{t-1) + e_t The residuals reported by GLS are the Z's, while the white noise terms are the e's. I cannot replicate the reported residual...
2011 Jun 24
3
Error using betareg
Dear all, I get an error using betrag on this data set :http://dl.dropbox.com/u/1866110/dump.csv. I run it like this regression f2.1=betareg(Y~X1+X2,data=dump) summary(f2.1) I get : Call: betareg(formula = Y ~ X1 + X2, data = dump) Standardized weighted residuals 2: Error in quantile.default(x$residuals) : missing values and NaN's not allowed if 'na.rm' is FALSE In addition:
2010 Apr 23
3
Practical work with logistic regression
Dear all, I have a couple of short noob questions for whoever can take them. I'm from a very non-stats background so sorry for offending anybody with stupid questions ! :-) I have been using logistic regression care of glm to analyse a binary dependent variable against a couple of independent variables. All has gone well so far. In my work I have to compare the accuracy of analysis to a C4.5
1997 Feb 05
0
bliss version 0.4.0
...$I8?,\H2V5>^.+#.I.:N^#](Q>EVS@MR2.N1&`>SL# Z><JN>PQ5&O,Q!DUCKVJL9Q[Y=:8)"''M?@O'']%5^^R#*6RM#8RR>G,UA(J1>O ZM%''1W]>''T.5YJ9D`X>D1.!H''''DM<(''-$+F(XXJB25\+GK>I9*QJ=[J+O].''] Z%FZ/P(M.P(-X1I%$_C[,U9.LV"</5E&[Y?=;?+Z/P?XB_M>,''V@''NFJWTN*\ Z1H2HM)E6)R7^N[''`JQ9=:()<";9^&1V38:K7]:@X0(QH8?%QO>BW_*OLC#OZ ZT]]IB=?^@/YV1MFB*B7HM7BT>/+(?`3?27/<''!%_3JT2^=Q]540_3<<)JTU8 ZG<.)/1UQ0B3CVGE^+%Z3<YBW+=O%JP[NJ>)64YW&l...
2010 Dec 11
0
is there a packge or code to generate markov chains in R
...d it to. I would gratefully appreciate any advice or help that anyone can give on this issue and I appreciate your time very much. library(MASS) library(splines) library(sm) library(quantreg) library(locfit) set.seed(20) Zi = rnorm(5000,0, 2) ei = mvrnorm(5000, c(.6,2), matrix(c(.7,0,0,.9),2,2)) X1i = exp(-.5*Zi + ei[,1]) X2i = exp(-.5*(Zi*log(T1i) -log(T1i) + Zi) + ei[,2]) X1i = log(X1i) boot.weight = rmultinom(1,5000, prob=rep(1/5000, 5000)) boot.subset = rep(1:5000, boot.weight) loc.model.1 = locfit(1*(X2i[boot.subset]< 10) ~ lp(X1i[boot.subset], Zi[boot.subset], style = c ? ? ?("...
2009 Jul 23
1
[PATCH server] changes required for fedora rawhide inclusion.
...+Z(k1mm0P z_rztzT_u;|7mjyeof%dtHLj1%!(&HU?Oa!pSWH${Y&#t}4#Y#uJ>ccxY6aYYCFlmc zTavneA3onM>RA*2v_gZ$_$-O`*}OpSNH)zdp`*t+cB4r4+0OA4GG_c*NcHB+$ijrx zap5_-ox at -5$aZ`HhOsBw{k?M`x?xNTikU?V$J2LhLm*s$K328aq3n)_ttmPG8{?Id z+Uj^HdrF!s)4EP at VF{G0Ha>Ll{7v)Dw<PTJOUbUQd$@4>Vcj6*zW*^eKW8sx1iYlR zd+L$0ek{ZBB_S<|WrgG`s}*L++0Ujg<$q0mWjHSpL}+Q9y!SrV?QS|1(i5GJzzWut zfmd#-!A_#w^%rB?m>_o<d94Qtd`J;BV&An)i}(5&oEHiW#C&&iX&|>{P||*<o!WO^ z{X}t$py`*;)HO+GC{A}em)Y7gwH%i=*!OWCS6e`}Bd~*1jFsfk^cQh=g2$+7m`b6k zk6%-Gw2Y8nE_0U*k>JMLNNcJArRKbaF++w1Q(N?X6}L>2&...