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&...