similar to: Stacking matrices

I read somewhere that t.test is equivalent to a hypothesis testing for one way ANOVA. But I'm wondering how they are equivalent. In the following code, the p-value by t.test() is not the same from the value in the last command. Could somebody let me know where I am wrong? > set.seed(0) > N1=10 > N2=10 > x=rnorm(N1) > y=rnorm(N2) > t.test(x,y) Welch Two Sample t-test data:

On 9/29/2017 3:37 PM, Rolf Turner wrote: > On 30/09/17 07:45, JLucke at ria.buffalo.edu wrote: > > <SNIP> > >> >> The conceptual paradigm for R is only marginally commensurate with >> that of >> standard statistical software. >> You must immerse yourself in R to become proficient. > > Fortune nomination. For newer list members wondering what

On 30/09/17 07:45, JLucke at ria.buffalo.edu wrote: <SNIP> > > The conceptual paradigm for R is only marginally commensurate with that of > standard statistical software. > You must immerse yourself in R to become proficient. Fortune nomination. cheers, Rolf -- Technical Editor ANZJS Department of Statistics University of Auckland Phone: +64-9-373-7599 ext. 88276

I wish to second this approach to learning R. I tried for several years to translate other stat programs or provide parallel analyses with R. This dabbling-in-R approach did not work . When a transferred to a research unit that could ill afford commercial software, I devoted my entire time to doing everything in R. This was a difficult learning process, but I eventually became proficient

I am trying to install rjags. Although I can install the library, I cannot load the package. See below. Here is my session > sessionInfo() R version 2.9.2 (2009-08-24) i386-pc-mingw32 locale: LC_COLLATE=English_United States.1252;LC_CTYPE=English_United States.1252;LC_MONETARY=English_United States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252 attached base packages: [1] stats

I am unable to install the rggobi package. I do not have trouble installing other packages. Any advice or suggestions? >sessionInfo() R version 2.8.1 (2008-12-22) i386-pc-mingw32 locale: LC_COLLATE=English_United States.1252;LC_CTYPE=English_United States.1252;LC_MONETARY=English_United States.1252;LC_NUMERIC=C;LC_TIME=English_United States.1252 attached base packages: [1] stats

guRus: I have a variable "beta" as an argument to R's beta function. So essentially I have a case of beta(alpha, beta). What surprises me is that R doesn't barf on this stupid programming practice. R gets the right answer. How does R know "beta the variable" from "beta the function"? Joseph F. Lucke Senior Statistician Research Institute on Addictions

The Lattice auto.key argument has a bug in R.12.2. R version 2.12.2 (2011-02-25) Platform: i386-pc-mingw32/i386 (32-bit) .... other attached packages: [1] lattice_0.19-17 loaded via a namespace (and not attached): [1] grid_2.12.2 If I set up my plot parameters as require(lattice) superpose.line.settings <- trellis.par.get("superpose.line") str(superpose.line.settings)

Although I have installed R many times on many machines without problems, I have recently encountered a problem with accessing packages. R has been installed successfully with the "--internet2" option. The session info is below. > sessionInfo() R version 2.8.1 (2008-12-22) i386-pc-mingw32 locale: LC_COLLATE=English_United States.1252;LC_CTYPE=English_United

windows XP R 2.10 As pointed out by Prof. Venables and Ripley (MASS 4th edition, p275), the results obtained from lme using method="ML" and method="REML" are often different, especially for small datasets. Is there any way to determine which method is preferable for a given set of data? Thanks, john John David Sorkin M.D., Ph.D. Chief, Biostatistics and Informatics

Hi R-users, Can Maple function be exported to R? I have a jacobian matrix (4X4)?from maple in algebraic form which involve modified Bessel function of the first kind. I just wonder whether we can use algebraic form into R before the value of the parameters can be estimated. Thank you so much for your attention and help.

Hi Yves lavaan looks like a very nice package. From the tutorial introduction I see you create path diagrams for some of the models you describe. How did you do this? I don't see a function for this in the package. I know there is a path.diagram function in the sem package that uses dot to draw the diagram, but I've always found the layouts from dot somewhat strange for path diagrams

k=lm(y~x) summary(k) returns R^2=0.9994 lm(y~x) is supposed to find coef. a anb b in y=a*x+b l=lm(y~x+0) summary(l) returns R^2=0.9998 lm(y~x+0) is supposed to find coef. a in y=a*x+b while setting b=0 The question is why do I get better R^2, when it should be otherwise? Im sorry to use the word "MS exel" here, but I verified it in exel and it gives: R^2=0.9994 when y=a*x+b is used

Hi, How can I make an ANOVA if I haven't got all data set but I know the numbers of subjects for each group, the mean and di standard deviation for each group? Thanks Angelo Scozzarella

Dear R users, i try to integrate lgamma from 0 to Inf. But here i get the message "roundoff error is detected in the extrapolation table", if i use 1.0e120 instead of Inf the computation works, but this is against the suggestion of integrates help information to use Inf explicitly. Using stirlings approximation doesnt bring the solution too. ## Stirlings approximation lgammaApprox

I am working on a project that requires me to do very large factorial evaluations. On R the built in factorial function and the one I created both are not able to do factorials over 170. The first gives an error and mine return Inf. Is there a way to have R do these larger calculations (the calculator in accessories can do 10000 factorial and Maple can do even larger) -- View this message in

Would like to try my luck to see if I can catch your eyes. I was trying to do some contrasts within ANOVA. I searched the archive and found a clue posted by Steffen Katzner ( http://tolstoy.newcastle.edu.au/R/help/06/01/19385.html) I have three levels for a factor names "StdLot" and want to make three comparisons, 1 vs 2, 1 vs 3 and 2 vs 3. First,

Hi, How can I get the values of mean and median (not only points but values too) on the boxplot. I am using boxplot function from graphics package. Following is my data set > df [1] 5 1 1 0 0 10 38 47 2 5 0 28 5 8 81 21 12 9 1 12 2 4 22 3 > mean.val <- sapply(df,mean) > boxplot(df,las = 1,col = "light blue") > points(seq(df), mean.val,

Dear All,   #I have the following code   Dose<-1000 Tinf <-0.5 INTERVAL <-8 TIME8 <-matrix(c((0*INTERVAL):(1*INTERVAL))) TIME7 <-matrix(c((0*INTERVAL):(2*INTERVAL))) TIME6 <-matrix(c((0*INTERVAL):(3*INTERVAL))) TIME5 <-matrix(c((0*INTERVAL):(4*INTERVAL))) TIME4 <-matrix(c((0*INTERVAL):(5*INTERVAL))) TIME3 <-matrix(c((0*INTERVAL):(6*INTERVAL))) TIME2

Hello, I''m a newbie here... please... can anyone help me about FULL AJAX ? for example if i''m scaffolding a simple model and controller for product which only have 1 field for NAME:STRING with the command : script/generate scaffold Product name:string and for the generated code I like to turn it into a full ajax app, for example if i''m click the NEW or EDIT link it