[R] Confidence intervals in ANOVA

Hallo! I have the a model with 3 time points, 2 treatments and N subjects. I can calculate an ANOVA but I can not calculate the CI of the interaction term (time and treatment), which I need for a closer look at the effect of the treatment to the 3 time points. I do NOT want to use lme

[R] Add row to data frame

Hi all, is there an easy way to build up a data frame by sequentially adding individual rows? The data frame consists of numeric and character columns. I thought of rbind, but I ended up with numeric values for the character columns. Pascal __

[R] graphical parametres...

Hi...it is probably trivial, but I do not know how to do the following: I want the output of a xy plot to be plotted in different colors according to a given condition ... I want to plot temperature dependency of flow stress ; for some (rare) occurencies I've got a special condition (chemical

Re: [R] Add row to data frame

Pascal == Pascal A Niklaus [EMAIL PROTECTED] on Mon, 08 Dec 2003 11:47:02 +0100 writes: Pascal Hi all, is there an easy way to build up a data Pascal frame by sequentially adding individual rows? yes, pretty easy, but usually not recommended because quite inefficient. rbind()

[R] help

Hi, there, I have a question about multinom rountine. The response is a matrix, and the predicator is a data frame called dat. I would like multinom to get response and data from within this function test instead of from R Console session.

[R] tkrplot with grid lattice plots

Hello, i tried to use a lattice graphics with tkrplot it seems that it doen't work, here is the exemple used: library(tkrplot) tt - tktoplevel() tktitle(tt)-Exemple randdata-data.frame(x=rnorm(100), y=rnorm(100), idobs=rep(1:10, each=10)) plot.graph-function() { + plot(randdata$x,

Re: [R] graphical parametres...

Anne Piotet wrote: Hi...it is probably trivial, but I do not know how to do the following: I want the output of a xy plot to be plotted in different colors according to a given condition ... I want to plot temperature dependency of flow stress ; for some (rare) occurencies I've got a special

Re: [R] Help

Stephen Opiyo wrote: Hi, I have a data set (data frame) approx. 50 rows * 600 columns. The columns are separated by commas. I would like to know how to remove those commas between the columns. What should I do to remove those commas? I guess during the import of the data? Please read the

[R] Reading formated databases

Hi R-users: How can I read an ascii database that is controled by the column number? For example: 7349593Luis Miguel Ariza Gutierrez 32342123 9394583X XX 34234930 39483 CCC CC39203230 3484932YY

[R] aggregate and names of factors

Hello, I use the function 'aggregate' a lot. One small annoyance is that it is necessary to name the factors in the 'by' list to get the names in the resulting data.frame (else, they appear as Group.1, Group.2...etc). For example, I am forced to write: aggregate(y,list(f1=f1,f2=f2),mean)

Re: [R] extracting p value from GEE

Hi, Many thanks for your kind help. best regards, Yu-Kang From: Emmanuel Paradis [EMAIL PROTECTED] To: Tu Yu-Kang [EMAIL PROTECTED] CC: [EMAIL PROTECTED] Subject: Re: [R] extracting p value from GEE Date: Fri, 05 Dec 2003 15:48:41 +0100 At 11:53 04/12/2003 +, vous avez rit: Dear R users,

Re: [R] aggregate and names of factors

Christophe Pallier [EMAIL PROTECTED] writes: Hello, I use the function 'aggregate' a lot. One small annoyance is that it is necessary to name the factors in the 'by' list to get the names in the resulting data.frame (else, they appear as Group.1, Group.2...etc). For example, I am forced

Re: [R] Reading formated databases

Kenneth Cabrera wrote: How can I read an ascii database that is controled by the column number? ?read.fwf -- Chuck Cleland, Ph.D. NDRI, Inc. 71 West 23rd Street, 8th floor New York, NY 10010 tel: (212) 845-4495 (Tu, Th) tel: (732) 452-1424 (M, W, F) fax: (917) 438-0894

Re: [R] Reading formated databases

Kenneth Cabrera wrote: Hi R-users: How can I read an ascii database that is controled by the column number? For example: 7349593Luis Miguel Ariza Gutierrez 32342123 9394583X XX 34234930 39483 CCC

[R] thanks!! get mean of several rows

Dear all! Thanks to all who replied to my question on getting the means of several rows, and the one with the standard error + mean-plot! Many of them worked fine just as they were, others had to be adapted a bit. However, I can finally do my calculations, and find myself happy as a man could

Re: [R] Add row to data frame

Martin Maechler wrote: Pascal == Pascal A Niklaus [EMAIL PROTECTED] on Mon, 08 Dec 2003 11:47:02 +0100 writes: Pascal Hi all, is there an easy way to build up a data Pascal frame by sequentially adding individual rows? yes, pretty easy, but usually not recommended because

Re: [R] Add row to data frame

On Mon, 8 Dec 2003, Pascal A. Niklaus wrote: Martin Maechler wrote: Pascal == Pascal A Niklaus [EMAIL PROTECTED] on Mon, 08 Dec 2003 11:47:02 +0100 writes: Pascal Hi all, is there an easy way to build up a data Pascal frame by sequentially adding individual

[R] TukeyHSD changes if I create interaction term

Dear R community, I'm trying to understand this behavior of TukeyHSD. My goal is to obtain defensible, labelled multiple comparisons of an interaction term. Firstly, if I plot the TukeyHSD from the model that calculates its own interactions, then the y-axis labels appear to be reflected on

[R] stripchart problem

Hello, I am trying to plot age distribution data for a certain condition that runs in families. Below is a simplified view of the dataset, i.e. in this case there are four families, each line corresponding to one individual with age at diagnosis and sex. famdata family age sex 1fam1 2.1

[R] WinMenus - is there a way of knowing if a WinMenu or WinMenuItem already exists?

I am developing a menusystem using the functions WinMenuAdd and WinMenuAddItem etc. I want to be able to shift between different interfaces (ie. different sets of menu-trees). Therefore I would like to be able to ask whether a specific menu already exists in order to remove or add it Without

Re: [R] WinMenus - is there a way of knowing if a WinMenu or WinMenuItem already exists?

On Mon, 8 Dec 2003, Finn Sando wrote: I am developing a menusystem using the functions WinMenuAdd and WinMenuAddItem etc. I want to be able to shift between different interfaces (ie. different sets of menu-trees). Therefore I would like to be able to ask whether a specific menu already

Re: [R] stripchart problem

Hello, I am trying to plot age distribution data for a certain condition that runs in families. Below is a simplified view of the dataset, i.e. in this case there are four families, each line corresponding to one individual with age at diagnosis and sex. famdata family age sex 1

Re: [R] Windows Memory Issues

On Sat, 6 Dec 2003, Prof Brian Ripley wrote: I think you misunderstand how R uses memory. gc() does not free up all the memory used for the objects it frees, and repeated calls will free more. Don't speculate about how memory management works: do your homework! Are you saying that

[R] test for arima coef's significancy

Dear sirs, I would like to know if there is a function to compute the pvalue for the significancy of arima coef in an arima object created by the arima function. I have written this one: pvalueArima-function(x,arima) { t-(arima$coef)/(diag(arima$var.coef)^0.5) df-length(x)-length(arima$coef)

[R] trouble with predict.l1ce

Dear R-help, I am having trouble with the predict function in lasso2. For example: data(Iowa) l1c.I - l1ce(Yield ~ ., Iowa, bound = 10, absolute.t=TRUE) predict (l1c.I) # this works is fine predict (l1c.I,Iowa) Error in eval(exper,envir, enclos) : couldn't find function Yield

[R] Recoding problem

Hello I have the following variables, all of which are logicals fmar15 fcoc15 fher15fcrk15fidu15 what I would like is a variable drug15 which equals idu if fidu15 is T; crk if fidu15 is F but fcrk is T, her if fher15 is T but fcrk15 and fidu15 are F and so on What's the best way

Re: [R] test for arima coef's significancy

On Mon, 8 Dec 2003, STOLIAROFF VINCENT wrote: Dear sirs, I would like to know if there is a function to compute the pvalue for the significancy of arima coef in an arima object created by the arima function. I have written this one: pvalueArima-function(x,arima) {

RE: [R] stripchart problem

stripchart always plots all the points in a given group with the same symbol, so you can't do what you want with it. Here are some alternatives: 1. Its not very nice but coplot can get you a chart somewhat in the vein you are looking for: with(famdata, coplot(age~family|family,

RE: [R] test for arima coef's significancy

Can you explain how you managed to derive a t distribution for this statistic, yet none of the references mentioned in the various help pages contain such a result? OK, let's say it was a naive and not well thought attempt. I tried it because I used to work with SAS and with the proc ARIMA, I

[R] Durbin Watson

the Durbin Watson function is in the car library. thanks, Erin __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help

RE: [R] Durbin Watson test

help.search(durbin) on my box gives: durbin.watson(car) Durbin-Watson Test for Autocorrelated Errors dwtest(lmtest) Durbin-Watson Test so you'll need either the `car' or `lmtest' package. HTH, Andy From: Of Erin Hodgess Hi R People: Where is the Durbin Watson test

[R] Durbin Watson thanks!

Hi all The Durbin Watson is also in the lmtest library as dwtest. Thanks to all of you who answered so promptly!! R Help rocks! Sincerely, Erin __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help

Re: [R] Recoding problem

I have the following variables, all of which are logicals fmar15 fcoc15 fher15fcrk15fidu15 what I would like is a variable drug15 which equals idu if fidu15 is T; crk if fidu15 is F but fcrk is T, her if fher15 is T but fcrk15 and fidu15 are F and so on What's the best way

[R] Matrix to Dates

Hello again R People: If I have a matrix with 2 columns z1 1960 1 1960 9 1961 6 __ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-help

[R] Matrix to dates

Let's try again! I have a matrix in which the first column is a four digit year, and the second column is a 2 digit month. How do I convert the matrix to a date function, please? Thanks, Erin Version 1.8.0 mailto: [EMAIL PROTECTED] __ [EMAIL

RE: [R] Matrix to dates

See ?DateTimeClasses, ?strptime, and ?as.character This example from strptime should get you going: ## read in date/time info in format 'm/d/y h:m:s' dates - c(02/27/92, 02/27/92, 01/14/92, 02/28/92, 02/01/92) times - c(23:03:20, 22:29:56, 01:03:30,

[R] R^2 analogue in polr() and prerequisites for polr()

Hi (1)In polr(), is there any way to calculate a pseudo analogue to the R^2. Just for use as a purely descriptive statistic of the goodness of fit? (2) And: what are the assumptions which must be fulfilled, so that the results of polr() (t-values, etc.) are valid? How can I test these

[R] Frequent crash printing graphics windows and wavethresh

I often have R (1.8.1) crash after I generate several graphics windows using windows() or X11(), print a graphics window and then rerun the same script. A windows message comes up saying Program error Rterm.exe .. (I can get the same problem using Rgui.) Xemacs tells me Process R trace trap at

[R] p-value from chisq.test working strangely on 1.8.1

Hello everybody, I'm seeing some strange behavior on R 1.8.1 on Intel/Linux compiled with gcc 3.2.2. The p-value calculated from the chisq.test function is incorrect for some input values: chisq.test(matrix(c(0, 1, 1, 12555), 2, 2), simulate.p.value=TRUE) Pearson's Chi-squared

Re: [R] Font

Use 'pointsize' argument, see ?postscript. This will proportionately scale all plot elements: title, axes labels and annotations, and plotting symbols (pch). If you want to control them separately, you may have to use the cex* parameters directly: cex, cex.axis, cex.lab, etc, see ?par. Best

[R] How to append to a data.frame?

Hi, I have a data.frame that I need to construct iteratively. At the moment, I'm doing: d-data.frame(x=c(),y=c(),z=()); # {and, within some loop} d-rbind(d,data.frame(x=newx,y=newy,z=newz); While this works, it is horribly verbose and probably not efficient, either. My real data.frame

Re: [R] R^2 analogue in polr() and prerequisites for polr()

many thanks! I was just asking for a r-square analogue, since the students I will present the results to, might like to know, how the measure of fit in an ordinal regression (e.g. the residual deviance) compare to measures they know (from introductory courses to linear regression) (such as the

Re: [R] How to append to a data.frame?

On Tue, 9 Dec 2003, David Kreil wrote: Hi, I have a data.frame that I need to construct iteratively. At the moment, I'm doing: d-data.frame(x=c(),y=c(),z=()); # {and, within some loop} d-rbind(d,data.frame(x=newx,y=newy,z=newz); While this works, it is horribly verbose