similar to: bcanon bug???

I am trying to implement an ANOVA on a pair of quantile regression models in R. The anova.rq() function performs a basic check to see whether the models are nested, but I think this check is failing in my case. I think my models are nested despite the anova.rqlist() function saying otherwise. Here is an example where the GLM ANOVA regards the models as nested, but the quantile regression ANOVA

I'm trying to figure out how to use lmer to fit models with factors that have some nesting and some non-nested groupings. For example, in this paper: http://www.stat.columbia.edu/~gelman/research/published/parkgelmanbafumi.pdf we have a logistic regression of survey respondents' political preferences (1=Republican, 0=Democrat), regressing on sex, ethnicity, state (51 states within 5

Dear experienced R-users, I am having some probably trivial trouble estimating the confidence interval for the difference of two group means, with groups been of unequal sample size. I am using the "Bootstrap" package and the function "bcanon"(bcanon(x, nboot, theta, ...,alpha=c(0.025, 0.05, 0.1, 0.16, 0.84, 0.9, 0.95, 0.975)) for Nonparametric BCa confidence limits. The

I have a question regarding bootstrap coverage. I am trying to understand the benefits of using the bootstrap for small sample sets. To do this I created a normal population and then picked 10 from the populations and applied both traditional statistical methods and the Bootstrap (bcanon, 5000 bootstrap samples) to calculate a 95% confidence interval of on the mean. I saved the width of the

Hi everyone, I want to plot some pattern. For example, I have some coding as below, plot(0:11,type="n") polygon(c(6,8,10), c(7,3,5),bor=6) polygon(c(2,4,6), c(5,3,7),bor=6) polygon(c(6,10,8), c(7,9,11),bor=6) polygon(c(2,6,4), c(9,7,11),bor=6) This coding will give me four triangles, but how to generate more triangles in sequences so that not need to mention the coordinate for each

I am looking for a function like my.blockwisesum(vector, n) that computes sums of disjoint subsequences of length n from vector and can work with vector lengths that are not a multiple of n. It should give me for instance my.blockwisesum(1:10, 3) == c(6, 15, 24, 10) Is there a builtin function that can do this? One could do it by coercing the vector into a matrix of width n, and then use

I am having some problems adding packages under FreeBSD 2.2.8 and R-0.62.2 I have included the output of "uname -a" below first: == FreeBSD stiegl.niksun.com 2.2.8-RELEASE FreeBSD 2.2.8-RELEASE #0: Thu Dec 10 10:29:07 EST 1998 ath at stiegl.niksun.com:/scratch/obj/compile/STIEGL i386 == The CC in use is gcc (not egcs). The install failed as follows: == stiegl> R INSTALL

Hello everybody, I have a bunch of strings like this: "IBM POWER4+ 1.9GHz" "IBM RS64-III 500MHz" "IBM RS64-IV 600 MHz" "IBM RS64 IV 750MHz" "Intel Itanium 2 Processor 6M 1.5GHz" "Intel Itanium2 1 Ghz" "Intel Itanium2 1.5GHz" "Intel

Dear R users, I am new to R. I would like to calculate bootstrap confidence intervals using the BCa method for a parameter of interest. My situation is this: I already have a set of 1000 bootstrap replicates created from my original data set. I have already calculated the statistic of interest for each bootstrap replicate, and have also calculated the mean for this statistic across all the

Hi list, i noticed that when doing imap gssapi authentication with kerberos, dovecot (here 2.1.7) always searches /etc/krb5.keytab although i have auth_krb5_keytab = /etc/mail3.krb5.keytab in my etc/dovecot/dovecot.conf and doveconf -n also show this setting. If i combine the keytabs in krb5.keytab it works. Is there another location where i should put my configuration regarding gssapi/kerberos ?

Brian Ripley wrote (to somebody asking about "effect sizes"): > ... > Given that, I wonder if you are used to standard terminology. Good point. But I think for many of us there is more behind that. I personally belong to an (apparently fairly large) group of R users who may be enthusiastic, but are statistical laymen due to a lack of formal education in the area. The

Look at help(sink). I have figured out how to fix up draw.tree to do factor labels correctly and will be doing that in the next few days. Meanwhile I discovered that rpart() provides factor labels in a much more satisfactory way than tree() does. Because of this and because of the additional information rpart provides about the tree building, I recommend you switch to that. > I have created

The build of R-1.4.0 with Atlas 3.3.13, both compiled with gcc 2.95.3 on an Athlon Thunderbird (1.333 GHz) running otherwise stock Red Hat Linux 7.2, appears to have been successful. R passed all its built-in checks. This is the first time I've been able to get all of this to work together; the last few times I tried this (late June or early July IIRC), R did not pass its checks when linked

I'm trying to extend some work in the car and heplots packages that requires getting a table of multivariate means for one (or later, more) terms in an mlm object. I can do this for concrete examples, using aggregate(), but can't figure out how to generalize it. I want to return a result that has the factor-level combinations as rownames, and the means as the body of the table

Hi all. Belove is the example from the cluster-help page wtih the output. I simply cannot figure out how to relate the estimate and robust Std. Err to the p-value. I am aware this a marginal model applying the sandwich estimator using (here I guess) an emperical (unstructered/exchangeable?) ICC. Shouldent it be, at least to some extend, comparable to the robust z-test, for rx :

Hello, I did a cross-validation using cvlm from DAAG package but wasn't sure how to assess the result. Does this result means my model is a good model? I understand that the overall ms is the mean of sum of squares. But is 0.0987 a good number? The response (i.e. gailRel5yr) has min,1st Quantile, median, mean and 3rd Quantile, and max as follows: (0.462, 0.628, 0.806, 0.896, 1.000, 2.400) ?

I have created a tree and want to save some of the data so that I can create a html table from it. I would like to save the output from data.ltr (see example below) to a file, but haven't found a way to do that, keeping the nice format that typing data.ltr gives me (see output below). Is there a way to do this? Example: library (maptree) library (tree)

Hi All. While wandering in auth-pam.c I noticed that there's a few Portable-specific escapees from the xmalloc(foo * bar) cleanup. There's also a "probably can't happen" integer overflow in ssh-rand-helper.c with the memset: num_cmds = 64; - entcmd = xmalloc(num_cmds * sizeof(entropy_cmd_t)); + entcmd = xcalloc(num_cmds, sizeof(entropy_cmd_t));

Hello, I am trying to plot means and standard errors conditioned by a factor, using qplot. I am successful at getting the bar graph I want with a error bar, however I have tried many things and cannot get the bars to change colors. Currently showing as red and blue, but need it to be black and white for publication. Any suggestions please? Using a data set June, which is str:

R Help - I'm trying to run a correlation matrix with a covariate of "age" and will at some point will also want to covary other variables concurrently. I'm using the "psych" package and have tried other methods such as writing a loop to extract semi-partial correlations, but it does not seem to be working. How can I accomplish this? library(psych) > set.cor(y =