similar to: How to make a matrix of a number of factors?

I have the following code. I'm wondering why summary() doesn't show F value and Pr? Rscript multi_factor.R > a=3 > b=4 > c=5 > d=6 > e=7 > > A=1:a > B=1:b > C=1:c > D=1:d > E=1:e > > X=matrix(nr=a*b*c*d*e,nc=5) > colnames(X)=LETTERS[1:5] > > for(i_a in 1:a-1) { + for(i_b in 1:b-1) { + for(i_c in 1:c-1) { + for(i_d in 1:d-1) { +

I have the R code at the end. The last command gives me "1 observation deleted due to missingness". I don't understand what this error message. Could somebody help me understand it and how to fix the problem? > summary(afit) Df Sum Sq Mean Sq F value Pr(>F) A 2 0.328 0.16382 0.1899 0.82727 B 3 2.882 0.96057 1.1136 0.34644 C

When run from the console, the following code works, when run through the WEBrick server I get a "You have a nil object where you didn''t expect it!" class Competition < ActiveRecord::Base has_many :comp_dates, :order=>:position def initial_date comp_dates.first end def display_dates i_d = initial_date if i_d.nil? return "unknown"

Hi, I'm trying to set up R to run a simulation of two populations in which every 3.5 days, the initial value of one of the populations is reset to 1.5. I'm simulation an experiment we did in which we fed Daphnia populations twice a week with algae, so I want the initial value of the algal population to reset to 1.5 twice a week to simulate that feeding. I've use for loops and if/else

All- I am interested in estimating a parameter that is the starting value for an ODE model. That is, in the typical combined fitting procedure using nls and lsoda (alternatively rk4), I first defined the ODE model: minmod <- function(t, y, parms) { G <- y[1] X <- y[2] with(as.list(parms),{ I_t <- approx(time, I.input, t)$y dG <- -1*(p1 + X)*G +p1*G_b dX <-

I don't find a function to print a string to file. Would somebody let me know what function I should use?

Hi, I am wondering if some experienced users would help put the ready-to-run code of the examples in the manuals. It would help new users learn R faster by putting all the examples in an ready-to-run R script file. Can somebody help do so sometime and post the code along with the pdf manuals? http://cran.r-project.org/manuals.html Regards, Peng

I found this paper on ANOVA on unequal error variance. Has this be incorporated to any R package? Is there any textbook that discuss the problem of ANOVA on unequal error variance in general? http://www.jstor.org/stable/2532947?cookieSet=1

Some examples in the help page are too long to be copied from screen. Could somebody let me know some easy way on how to extract the example to a file so that I can play with them?

read.table terminates the program if the input file is empty. Is there way to let the program continue and return me a NULL instead of terminating the program? $ Rscript read_empty.R > read.table("empty_data.txt") Error in read.table("empty_data.txt") : no lines available in input Execution halted $ cat read_empty.R read.table("empty_data.txt") $ cat

Am 02.11.2013 um 11:38 schrieb Aris Adamantiadis <aris at 0xbadc0de.be>: > RFC4251 describes mpint to be multi-size and with positive values having > MSB clear, so it's clearly incompatible with raw string. > > Since you both agreed on the curve25519 implementation to use, I'll work > today on Markus' patch to make the changes Damien wanted. What do you want to

Hello . I have a little question .. I work with Anova (aov) and Manova (manova) and want to print some result in a text file. I see the results with the summary function or the summary.aov function but i do not know how to save the differents probabilities of error in a value. Sorry for my poor English Thanks a lot. -- Lamoral Julien

> > but doesn't emit optnone and nounwind attributes > s/nounwind/noinline On Mon, Jan 28, 2019 at 11:35 AM Jakub (Kuba) Kuderski < kubakuderski at gmail.com> wrote: > As far as I understand, gllvm doesn't run LTO pipeline or any cross-module > optimization, and the optimization level provided is only used to compile > each Translation Unit separately. > If you

Hi, I have llvm installed through MacPorts. But I don't see macports installed. Are there supposed to be manpages of commands like lli installed? -- Regards, Peng

Hi, I am wondering how to pop up the graphics window from Rscript. I run the following code, but I don't see the graphics window, even transiently. Regards, Peng $ Rscript plot.R > x=1:10 > y=1:10 > plot(x,y) >

I'm thinking of using external program 'grep' and pipe() to do so. But I'm wondering if there is a more efficient way to do so purely in R

I'm wondering where I can find the detailed descriptions on R memory management. Understanding this could help me understand the runtime of R program. For example, depending on how memory is allocated (either allocate a chuck of memory that is more than necessary for the current use, or allocate the memory that is just enough for the current use), the performance of the following program could

I can exclude columns by column number using '-'. But I wondering if there is an easy way to exclude some columns by column names. > x=cbind(c(1,2),c(3,4)) > x [,1] [,2] [1,] 1 3 [2,] 2 4 > colnames(x)=c('a','b') > x a b [1,] 1 3 [2,] 2 4 > x[,-'a'] Error in -"a" : invalid argument to unary operator > x[,-1] [1] 3 4

length(unique(c(1,2,2)))==length(c(1,2,2)) I use the above test to test if all the elements in a vector are unique. But I'm wondering if there is a convenient function to do so in R library.

Hi, x=matrix(1:60,nr=6) I can refer the last 2 rows by x[5:6,] If I don't know the total number of rows is 6, is there a way to refer the last 2 rows? Regards, Peng