similar to: Sweave in makefile under Win XP

Hi all. I think it would be nice to be able to combine levels of a factor on creation a la x <- rep(0:5,5) y <- factor(x,levels=0:5,labels=c('1','1',2:5)) ## (1) y [1] 1 1 2 3 4 5 1 1 2 3 4 5 1 1 2 3 4 5 1 1 2 3 4 5 1 1 2 3 4 5 Levels: 1 1 2 3 4 5 I thougt this would (should?) create a factor with 5 levels, ie combining 0 and 1 in x into one level in y. I find it

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 :

Hi helpeRs. When I do two-ways tables dimnames in columns are almost left-adjusted: <- factor(rbinom(30,1,.3),labels=c("first","second")) > b <- a > table(a,b) b a first second first 21 0 second 0 9 If I do a oneway table however, dimnames are right-adjusted: table(a) a first second 21 9 Actually I must admit I

Hi all. After upgrading to R2.0.1 i get > "??" [1] "??" > "??" [1] "\370" > "??" [1] "??" Whereas under R1.9.1 i get > "??" [1] "??" > "??" [1] "??" > "??" [1] "??" Any hints apreciated. Steen Steen Ladelund, statistician +4543233275

Hi, I need to produce a standard report for several variables in Sweave and thus would need the possibility to define a TeX-command which includes R-code like \newcommand{\meansd}[1]{The mean is \Sexpr{mean(#1)} and the standard deviation is \Sexpr{sd(#1)} . } and then just write \meansd{age} in the latex code to get the whole sentence. The above does not work, since

I would like to plot two data sets (frequency (y-axis) of mean values for 0-1(x=axis)) on a single histogram for comparison. The hist() only allow the overlay of two histograms, and although barplot() allows beside=TRUE, it does not show frequency values (like hist) but rather all of the values. Is there any way that I can use the hist() to plot two data sets similar to the barplot(). Any help or

Hi. I've done a standard barplotexample using R1-8.0, redhat8, probably using a GhostScript-integration. How can I improve the "sharpness" of the fonts. See for example the barplot (originally jpg coming from R) which is part of this screendump: http://www.zug.dk/Members/Nielssteenkrogh/1067525916309553166/element_103003161029 I'm not sure where to improve the this. Is it an

I've got a bunch of similar datasets, all of which I've fit to linear models. I'd like to easily create arrays of a specific parameter from each linear model (e.g., all of the intercepts in one array). I figured I'd put the model objects into an array, and then (somehow) I could easily create corresponding arrays of intercepts or residuals or whatever, but I can't the

We are constructing growth charts (age/weight and age/length) for children with diagnosis that impacts weight/length. But we we don't know how to use R for producing growth charts. We are collection data of Age, Weight and Length. The data are used to produce diagnosis-specific Growth charts (like the CDC Growth Charts:

I try to sort this dataframe: [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [1,] "CM" "BARBY" "INCREASED" " 0" " 2" " 0" " 1" " 1" [2,] "CM" "BARBY" "REDUCED" " 0" " 1" " 2" " 2" " 5" [3,] "CM"

Dear R-list. I'm doing af logistic analyses using gml. The model explaines variations in Adverse events infections (0 og 1) using age as explanatory variable. model2d<-glm(formula=AEorSAEInfecBac~Age,family=binomial("logit"),data=emrisk) I want to get predictions with 95% confidence limits for age 30 and age 60. I've been reading the "google" and "search

I administer a nano-setup on a dorm-network with a couple of hundred active users. The setup uses 2 x 2 2Mb/s DSLs, meaning two DSLs from each of two different ISPs. It works fine except for some minor glitches: https-sites often kicks users. This was solved by tying outbound https to a single DSL. Not the best solution but it works so far that users dont kicked from the sites anymore. Now

[image001.gif] Diploma's en titels zonder studie op basis van eerder verworven competenties Geachte heer of mevrouw, Graag willen wij u attenderen op onze unieke dienstverlening. Wellicht dat er voor u mogelijkheden zijn om zonder studie op basis van eerder verworven competenties via een External Programme in het bezit te komen van een academische Bachelor, Master

I''m not a programmer per se (including C), but... I''d like to be able to give the define in sch_sfq.c (of, say, 2.4.20), SFS_DEPTH other values than 128 as an argument on the tc commandline. It could be powers of two up to 2^7 (128) as it seems that 128 is the current maximum. I''m a little anxious to ask the question "How do I do that?" ;-) Instead

Hello, Just to let everyone know that this is still an issue. I am trying to install FreeBSD 9.0 amd64 on a Lenovo X121e and I can't get it to accept the geli passphrase during boot. I've confirmed using kern.geom.eli.visible_passphrase=1 that the passphrase is correct, and the same passphrase is accepted when the system is booted up. I've tried disabling kbdmux in

I am using R-2-15-1 with latest version of RTools under Windows 7 on a 64bit machine. When trying to build a dll from a simple test program matm.f I get the following error: > R CMD SHLIB matm.f gfortran -O3 -mtune=core2 -c matm.f -o matm.o f951.exe: out of memory allocating 312 bytes make: *** [matm.o] Error 1 Can someone please help NJ ________________________________

I am trying to make a barplot with a broken axis using gap.barplot (in the indispensable plotrix package). This works well when the data is a vector: > twogrp<-c(rnorm(10)+4,rnorm(10)+20) > gap.barplot(twogrp,gap=c(8,16),xlab="Index",ytics=c(3,6,17,20),ylab="Group values",main="Barplot with gap") But when the data is an array (for a bar plot with multiple

The following example, extracted from pp. 38-42 of Chambers (1998) Programming with Data (Springer), works for me in S-Plus 6.2 but not R 1.9.0pat, which returns the indicated error messages: > setClass('track', representation(x="numeric", y="numeric")) [1] "track" > tr1 <- new("track", x=1:3, y=4:6) > setMethod("plot", +

Hello! I want to thank everyone who helped me out with tips for load balancing asterisk machines in a cluster. I have encountered a new problem that is related to SIP diversion headers in the INVITE. I make calls through the manager interface and now want to add a SIP-Diversion header that changes the CallerID of a number that is not available on the trunk, the CallerID to be visible externally

Hi, I'd like to know if is it possible to make a 3D scatter plot (x y z) using R. In the on-line help and in the help mailing list I found only ways in which plot surfaces. Thanks for the help fulvio -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list -- Read http://www.ci.tuwien.ac.at/~hornik/R/R-FAQ.html Send