similar to: interpreting one-way anova tables

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

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

Hello, I'm new to R and have been following many guides including the two-way anova (http://www.personality-project.org/r/r.anova.html). Using that walkthrough including the supplied data I do str(data.ex2) and receive the appropriate types of data as follows: > str(data.ex2) 'data.frame': 16 obs. of 4 variables: $ Observation: int 1 2 3 4 5 6 7 8 9 10 ... $ Gender :

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 =

I have a vector of 2,1,0 I want to change to 0,1,2 respectively (the data is allele dosages) I have tried multiple nested if/else statements and looked at the ?if help and cannot work out what is wrong, other people have posted code which is identical and they state works. Any help would be greatly appreciated. > A[1:20] [1] 1 1 0 0 1 0 1 0 1 0 0 0 1 1 0 1 1 1 0 0 > B <-

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

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

I note that "current implementations of R use 32-bit integers for integer vectors," but I am working with large arrays that contain integers from 0 to 3, so they could be stored as unsigned 8-bit integers. Can R do this? (FYI -- This is for storing minor-allele counts for genetic studies. There are 0, 1 or 2 minor alleles and 3 would represent missing.) It is theoretically possible

> From: Thomas Baier > > > I am confused as to what I need and must do. > > > > Your message said > > ftp://cran.r-project.org/pub/R/src/contrib/Devel/rcom_0.97.tar.gz > > > > The DESCRIPTION file inside that download says > > Package: rcom > > Version: 0.95 > > Date: 2004-07-14 > > That's too bad. Sorry for that. > >

Everyone is back at work and starting to use computers on our smallest network which has Internet access through a rather old Linksys router.? Two systems were left on and screen-locked over the extra long weekend.? There does not appear to have been any Internet access interruption in our absence. A Firefox browser on one system was left pointing to a commonly used web site:

I am trying to use the SeasonalMannKendall function in the Kendall package. My dataset (alb_data) is in the same format as the example dataset (manaus) in the package. > class(manaus) [1] "ts" > is.ts(manaus) [1] TRUE > typeof(manaus) [1] "double" > alb_data=read.table("R:/albemarle_manken.txt", header=T) > head(alb_data) year Jan Feb

Dear Rxperts, Running the following code: ======================================================= twlo=10; twhi=20; wt=154; vd=0.5; cl=0.046; tau=6; t=3; F=1; wtkg <- wt/2.2 # convert lbs to kg vd.pt <- wtkg * vd # compute weight-based vd (L) cl.pt <- wtkg * cl # compute CL (L/hr) k <- cl.pt/vd.pt # compute k (hr^-1) cave <-

Dear prof. Harrell, I'm not able to use the force option with fastbw, here an example of the error I've got (dataset stagec rpart package): > fitstc <- cph(Surv(stagec$pgtime,stagec$pgstat) ~ age + eet + g2 + grade + gleason + ploidy, data=stagec) > fbwstc <- fastbw(fitstc,rule="aic",type="individual") > fbwstc Deleted Chi-Sq d.f. P Residual d.f.

I have a data set for different time intervals. The data has three comment lines before data for each time interval. For each time interval there are 500 data points. I want to change the dataset such that I have the following format: t1 t2 t3 ................ 0.00208 0.00417 0.00625 ................. a1 a2 a3 ...................

Hi there, I have a question regarding a data set that looks like the following: 8 repeated measurements for 3 dosage levels. However, the measurements for each dosage levels are taken on a different group of subjects. example: subject time dosage measure 11         1          1          154 11         2          1          150 :          :           :             : 11        

Hi, Does anybody know how to obtain the imputed SNP genotype probabilities from the snpStats package? I am interested in using an imputation method implemented in R to be further used in a simulation study context. I have found the snpStats package that seems to contain suitable functions to do so. As far as I could find out from the package vignette examples and its help, it gives the

Dear all, I often import data sets from other programs that contain user defined missing values. As an example consider a variable MXPLOCO with a user defined missing value of as -9. > summary(MXPLOCO) Min. 1st Qu. Median Mean 3rd Qu. Max. NA's -9.000 4.750 5.000 5.349 6.000 8.500 2.000 The following command successfully converts the value -9 to NA. However,

Dear All: Could you please tell me how I can draw figure formatted like figure 4.18 of MASS (4th) with the attached data set? Thanks Zhongming Yang --------------------------------- -------------- next part -------------- An embedded and charset-unspecified text was scrubbed... Name: sample.txt Url: https://stat.ethz.ch/pipermail/r-help/attachments/20050316/abfdb85e/sample.txt

Dear all: I m a newer on R.? I have some problem when I use?anova function.? I use anova function to get Type I SS results, but I also need to get Type III SS results.? However, in my code, there is some different between the result of Type I SS and Type III SS.? I don?t know why the ?seqe? factor disappeared in the result of Type III SS.? How can I do?? Here is my example and result.