similar to: Simple 2-Way Anova issue in R

Hi there, I checked google for aov. usually one uses summary to see whether the p-value is small. but I want to put aov in my script. how can I get the p-value, (0.1115, 0.6665, 0.6665 in the following example)? thanks YU > datafilename="http://personality-project.org/r/datasets/R.appendix2.data" > data.example2=read.table(datafilename,header=T) > aov.ex2 =

Hi, I am trying to reconcile anova table in R (summary(lm)) with individual t.test. datafilename="http://personality-project.org/R/datasets/R.appendix1.data" data.ex1=read.table(datafilename,header=T) #read the data into a table summary(lm(Alertness~Dosage,data=data.ex1)) gives: Call: lm(formula = Alertness ~ Dosage, data = data.ex1) Residuals:    Min     1Q Median     3Q    Max

Hi all, How does R know to regard a variable as a factor and not a character? For example, consider the following table: Observation Gender Dosage Alertness 1 m a 8 2 m a 12 3

Problem on flexmix when trying to apply signature developed in one model to a new sample. Dear R Users, R Core Team, I have a problem when trying to know the classification of the tested cases using two variables with the function of flexmix: After importing the database and creating a matrix: BM<-cbind(Data$var1,Data$var2) I see that the best model has 2 groups and use: ex2

Dear r heper, How can I change the strip text, for example (16,23] in the following example, to other more informative text such as "high level" on the fly? library(effects) Cowles$ex2 <- cut(Cowles$extraversion,3) mod.cowles <- glm(volunteer ~ sex+neuroticism*ex2,data=Cowles, family=binomial) eff.cowles <- allEffects(mod.cowles) plot(eff.cowles,

When investigating the legend(x,y, expression ( <math>, <math> ), ....) [as suggested by Jim Lindsey], I once again stumbled over the fact that expression objects behave somewhat ``wierdly'', IMHO. We can have expressions of length(.) > 1, and access elements of these like list elements. However, other natural things don't work with expressions; some example

I've got two tables.... first one(table1): ID chrom start end Ex1 2 152 180 Ex2 10 2000 2220 Ex3 15 3000 4000 second one ( table2): chrom location name 2 160 Alv 2 190 GNN 2 100

Hi everyone, I am beginner using R but I try to learn more. I need this graphic or similar but instead of tropical and temperate are three ontogenetic states Larva Met Juv CTMAX CTMIN SP SP SP This is my scrip: Ex = subset(Expr, Outlayer=="N") Ex2 = subset (Ex, S0 == 1) Ex3

Hello list, I want to make a large rulebased algorithm, to provide decision support for drug prescriptions. I have defined the algorithm in a function, with a for loop and many if statements. The structure should be as follows: 1. Iterate over a list of drug names. For each drug: 2. Get some drug related data (external dataset). Row of a dataframe. 3. Check if adaptions should be made to

Full_Name: Karsten Krug Version: 2.4.0 OS: Open Suse 10.0, Windows XP Submission from: (NULL) (88.134.13.50) I found a problem in the 'optimise' function for one dimensional optimisation. Example 1: Try to find a maximum of the function below with the use of 'optimise' in the interval [0,0.5]. The function follows a parabola and has two local maxima located at the margins of

I had an error message while running a macro from Yudi Pawitan's web site: > source("ex2-13.r") Error in parse(file, n, text, prompt) : syntax error on line 2 Inspecting ex2-13.r I found that the error was generated by another source() command. Clearly R does not like nested source()s, which is fair enough when you think about it. Still it's something that you might want

Using R2.12.1 on Ubuntu 10.04.1 I've tried to run the following code chunk in odfWeave <<fig1, echo=TRUE,fig=TRUE,width=7,height=4>>= x<-seq(1:100)/10 y<-sin(cos(x/pi)) imageDefs <- getImageDefs() imageDefs$dispWidth <- 4.5 imageDefs$dispHeight<- 4.5 setImageDefs(imageDefs) X11(type="cairo") plot(x,y) title(main="sin(cos(x/pi))")

Hello, I've been playing with the manova() function to do some pretty straightforward multivariate analyses, and I can't for the life of me figure out how to get at the discriminant functions used. When predicting several variables simultaneously, it's important to be able to gauge how much each variable is contributing to the analysis...a simple p-value isn't really enough. I

Dear all, I make habitat suitability models for animal species. The purpose of my research is to investigate the accuracy of different models. I clearly have a nested design: - accuracy_measure -> response variable - 2 model types (model_type) -> fixed effect - 230 species (species) -> random effect - 10 replicates/species (replicate) -> random effect - 10 subreplicate/replicate

Dear R-community, Using bwplot, how can I put the whiskers at percentile 5 and percentile 95, in place of the default position coef=1.5?? Using panel=panel.bwstrip, whiskerpos=0.05, from the package agsemisc gives satisfaction, but changes the appearance of my boxplot and works with an old version of R, what I don’t want, and I didn’t find the option in box.umbrella parameters Many thanks

Hi all. I'm working on an assignment for a psychology class and I am not sure how to adjust the y-axis so it displays the range: 0, 5, 10, 15 The code below is almost ideal: lsd = c(3, 5, 13) mar = c(1, 2, 3) g_range <- range(0, lsd, mar) plot(lsd, type="o", col="blue", ylim=g_range, axes=FALSE, ann=FALSE) axis(1, at=1:3,

Hello to everybody I have vector of drug dosage in string format 1.5 mg , 50 mg , 100 mg , 0.5 mg ... The task is to remove last free symbols: ' mg' Any suggestions are highly appreciated

In an Rweave tutorial written for possibly naive R users, I felt it necessary to explain why table() had to be used inside with(), whereas other tools like xtabs() had a data= argument. with() is quite nice for such cases, but it seems an unnecessary thing to learn right off. Before I turn this question into a request for R-devel, is there any inherent reason why it might be hard to add a

Hi list, I don't know if somebody has spent a lot of time debugging strange problems with if else positioning - the parser seems to recognize only the syntax bellow - this is the only way of making these pieces of code to work. As far as i'm concerned, no examples were available (it would be so awesome to have them in the introductory manual!) #Try to change the placement of the

Hi, I am using the: pcf.aov<-aov(meas~op+part, data=pcf.ex2), command to perform a two way ANOVA. When I save the: sumpcf<- summary.aov(pcf.aov), result of the summary.aov command in a variable I need to access the individual pieces of information in the summary. The summary appears to be a list and I am having a hard time finding a way to get at the information in the summary.