search for: rhohat

...oung, M Gail Woodbury and Karen Fryday-Field (1994): Statistical Methodology for the Concurrent Assessment of Intrarater and Intrarater Reliability: Using Goniometric Measurements as an Example. Physical Therapy 74 (8); 777-788 The function gives the estimated inter- and intrarater reliabilities (rhohat) for fixed and random rater effects, partial intrarater reliabilities for each rater, F test statistics with corresponding p-values, lower bounds for one-sided confiodence intervals and standard errors of measurement. The defaults are set up for use in a current project. Script for testing in the...

Dear R-users I need to calibrate kappa, rho, eta, theta, v0 in the following code, see below. However when I run it, I get: y <- function(kappahat, rhohat, etahat, thetahat, v0hat) {sum(difference(k, t, S0, X, r, implvol, q, kappahat, rhohat, etahat, thetahat, v0hat)^2)} > nlminb(start=list(kappa, rho, eta, theta, v0), objective = y, lower =lb, > upper =ub) Error in dots[[1L]][[1L]] : object of type 'closure' is not subsettable And I d...

I am working on a two stage threshold quantile regression model in R, and my aim is to estimate the threshold of the reduced-form equation (call it rhohat), and the threshold of the structural equation (call it qhat), in two stages. On the first stage, i estimate rhohat by quantile regression and obtain the fitted values. I use these fitted values to estimate qhat on the second stage. The code is as follows (thanks to Prof. Bruce Hansen, whose code i...

...quot;, "SPW", "SS", "SSB", "SSS", "SW", "SWF", "WB", "WLH", "WW") > names(data) <- species > cor.mat <- vcv.phylo(treeX, cor=TRUE) > regr <- compar.cheverud(data, cor.mat) > regr$rhohat [1] 5.541462 > 1 - var(regr$residuals)/var(data) [,1] [1,] -1.333095 I don't understand why the autoregression coefficient falls outside the interval -1 to 1, or why the calculation for obtaining an R-squared produces a value that doesn't make sense. Have I made a mistake i...

...ed, see table 2 (p. 780) and 3 (p. 281) >  sighat2e<-MSE >  sighat2Sfixed<-(MSS-MSE)/(mm*tt) >  sighat2Rfixed<-(MSR-MSSR)/(mm*nn) >  sighat2SRfixed<-(MSSR-MSE)/mm >  # the same for random and fixed, see table 2 (p. 780) and 3 (p. 281) >  sighat2e.part<-MSEpart >  rhohat.inter.random<-sighat2Srandom/ >  (sighat2Srandom+sighat2Rrandom+sighat2SRrandom+sighat2e) >  rhohat.inter.fixed<-(sighat2Sfixed-sighat2SRfixed/tt)/ >  (sighat2Sfixed+(tt-1)*sighat2SRfixed/tt+sighat2e) >  rhohat.intra.random<-(sighat2Srandom+sighat2Rrandom+sighat2SRrandom)/ >...