similar to: Do we have similar function to Minitab's INDICATOR?

Dear R-Friends, Sorry to bother everyone with my earlier question: "Do we have similar R function to work like Minitab's INDICATOR?" I found the way to make things work: e.g., > x <- c(2,2,5,3,6,5,NA) > model.matrix(~ factor(x) - 1) factor(x)2 factor(x)3 factor(x)5 factor(x)6 1 1 0 0 0 2 1 0 0 0

Hi I am using the following script to average a set of data 0f 62 columns into 31 colums. The data consists of values of ln(0.01) or -4.60517 instead of NA's. These need to be averaged for each row (i.e 2 values being averaged). What I would I need to change for me to meet the conditions: 1. If each run of the sample has a value, the average is given 2. If only one run of the sample has a

Hi All I am having 2 very long character strings (550chars) and I want to put them as expressions together with c(). The problem is that I also get these double-quotes, as seen below in 'fct'. How can I remove these double-quotes? I tried as.name() but it did not work (because of size?). These are creating trouble with subsequent programs, which I tested with strings that for some

Dear all, I found that the plotmath seems won't work on xlab in plot (and other plot functions including matplot, contour, ...): > plot(1:10, 1:10, xlab=expression(alpha)) the plot shows nothing for x-label. My R-version is R.1.1.0 > version _ platform Windows arch x86 os Win32 system x86, Win32 status major 1

Dear R-friends, For the following data: > xy x y i 1 731 0.313 2 2 739 0.340 2 3 790 0.373 2 4 855 0.451 2 5 980 0.608 2 6 575 0.156 1 7 608 0.207 1 8 630 0.249 1 9 670 0.332 1 10 838 0.377 1 11 964 0.466 1 > coplot(y ~ x|i, data=xy) coplot gives 3 panels, rather than 2, namely one for i=1 and one for i=2. Futhermore, when I extand data fram xy to have i=3 as follows:

Dear R and ESS users, I got the following error when trying to load a source file (smart.R in directory E:\R\) into XEmacs through the menu iESS and "Load source file": Error in file(file, "r") : unable to open connection In addition: Warning message: cannot open file `E:Rsmart.R' It seems the slash "/" or "\" being missing, but I don't know

Dear R-Friends, How can I calculate a double integral like \int_a^b \int_c^y g(x, y) dx dy where a, b, c are constants, g(x, y), e.g., g(x, y) = tan(x + y). I tried to nested integrate() and adapt(), but none of them working, seemingly due to the limits can not be specified constants. Best regards, C. Joseph Lu Department of Statistics National Cheng-Kung University Voice:

Hi, >> The patch in PR is wrong since it will make both darwin and windows >> unhappy. I will comment on gentoo PR about possible proper solution. > > I'll give it a try as you comment on the gentoo PR. > But my ppc box is in my office, so I'll post a result on next Monday (JST). I posted my result on the gentoo PR. https://bugs.gentoo.org/show_bug.cgi?id=403519#c12

Dear Friends: I like so much to work with R program. Congratulations for your work. I need R for work with multivariate data. My question is: With the pairs(X) command my output is a pairwise scatterplot symmetric matrix. Like: | X1 |X1 vs X2|X1 vs X3|X1 vs X4| |X2 vs X1| X2 |X2 vs X3|X2 vs X4| |X3 vs X1|X3 vs X2| X3 |X3 vs X4| |X4 vs X1|X4 vs X2|X4 vs X3| X4 | It is

I just thought of a useful metaphore for the problem I face. I am dealing with a problem in business finance, with two kinds of related events. However, imagine you have a known amount of carbon (so many kilograms), but you do not know what fraction is C14 (and thus radioactive). Only the C14 will give decay events (and once that event has occurred, the atom that decayed will never decay

Hello > I posted my result on the gentoo PR. > https://bugs.gentoo.org/show_bug.cgi?id=403519#c12 > https://bugs.gentoo.org/show_bug.cgi?id=403519#c13 > https://bugs.gentoo.org/show_bug.cgi?id=403519#c14 I commented there. You need to declare function like extern "C" { static void* LLVM_ATTRIBUTE_USED PPCCompilationCallbackC <everything else goes here> }; -- With

Dear R-friends, When I selected "Save as postscript" via mouse left-button on the active "R Graphics" window, those colored lines I made in the plot changed in the saved postscript file, e.g., color blue becomes green, red becomes blue, ... I have not changed any setting/options. What's going on? My version of R is > version _ platform i386-pc-mingw32 arch

Dear R-friends, As the binary windows version of R being distributed, what is the correct way to install contributed packages without using winst.exe as before? The answer in "R for Windows FAQ" says: "You can install pre-compiled packages either from CRAN or from a local .zip file by using install.packages: see its help page." But the information available in

Hello all, here's a collection of answers I got on my question concerning partial correlation coefficients: Some people gave a simple formula for the three-variable-case, as did Dave Lucy: pcor <- function(v1, v2, v3) { c12 <- cor(v1, v2) c23 <- cor(v2, v3) c13 <- cor(v1, v3) partial <- (c12-(c13*c23))/(sqrt(1-(c13^2)) * sqrt(1-(c23^2)))

Dear userR, With the following results, are they correct or acceptable? > x <- c(1.4, 1.2, 2.8) > sum(x) [1] 5.4 > sum(x) == 5.4 [1] FALSE > (1.4 + 1.2 + 2.8) - 5.4 [1] -8.881784e-16 > (1.4 + 1.2) - 2.6 [1] -4.440892e-16 > 2.6 - 1.5 - 1.1 [1] 0 > version _ platform i386-pc-mingw32 arch i386 os mingw32 system i386, mingw32 status major 2 minor

I am using R to create colour plots as bitmap files. I have created a palette of 256 colours that I want to use, but the resultant output files contain only 12 colours used in dither patterns. The result is the same whether I use bmp() or png(). How can I get the true colours that I specified instead of dither patterns? (making the driver's job easier!) I am using R 1.1.0 for Windows

Hi R-help. I am trying to run a linear mixed model with nested factors with either lme or lmer and I am having no luck obtaining the same results as Minitab. Here is Minitab's code: MTB > GLM 'count' = site year replicate(site year) site*year; SUBC> Random 'year' 'replicate'; Can you tell me how to code this in R? The settings are typeII, Tukey,

HI, Using c11<- 0.01 c12<- 0.01 c1<- 0.10 c2<- 0.10 One possible problem is that: dim(res5) #[1] 513? 20 res6<-aggregate(.~m1+n1+m+n,data=res5[,c(1:6,9:12,21:24)] ,max) #Error in `[.data.frame`(res5, , c(1:6, 9:12, 21:24)) : ?# undefined columns selected A.K. ________________________________ From: Joanna Zhang <zjoanna2013 at gmail.com> To: arun <smartpink111 at

Dear all, Does R or S-plus or any of their packages provide any command to form any of the following procedures to find Best Regression Equation - 1. 'All Possible Regressions Procedures' (is there any automated command to perform 2^p regressions and ordering according to criteria R2(adj), mallows Cp, s2- by not setting all the regression models manually), 2. 'Backward

Dear R-friends, After change (and saved) font style to "bold" in GUI preferences from menu-bar Edit, R 1.2.0, Win98, the "buflines" in the file Rconsole changed to "bufiles" and that causes error in invoking R next time. Is this a bug? Merry Christmas & Best regards, C. Joseph Lu Department of Statistics National Cheng-Kung University Tainan, Taiwan