similar to: Contrasts within ANOVA frame (Repost)

Dear all, 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". I just want to make three comparisons, 1 vs 2, 1 vs 3 and 2 vs 3. First, contrasts(d3$StdLot,3)<-matrix(c(1,-1, 0,0,1,-1,1,0,-1),3,3) #d3 is the data

If you go to http://SingingFish.com you can see that this consumated Dec. 1. This gives the MP3 consortium the most comprehensive access to web distributed MP3 streaming and content, (and other formats)... Full article: http://www.thomson-multimedia.com/vus/04/042/f042.htm THOMSON MULTIMEDIA TO PURCHASE SINGINGFISH.COM TO BUILD STREAMING MEDIA CAPABILITIES Leading Technology Helps Web

Hi, everyone, I was trying to sort a data frame by two columns, one increasing, the other decreasing and got an error. "Error in FUN(left) : invalid argument to unary operator", The command is "BA[order(BA[1],-BA[2]),]". BA is the data frame. It was working if I used increasing on both columns. Why the decreasing symbol "-" is not working here? Thanks. -- Jun

Hi, I have run a groups of t test and obtained a list of the t-test results (about 30). How do I extract p-values from the 30 t-test results and make a table output? Thanks Jun Example t-test output: ============================ $Lambda_z Welch Two Sample t-test data: x[PA1.BS.DN$Onecode == 2] and x[PA1.BS.DN$Onecode == 3] t = -0.1128, df = 33.635, p-value = 0.9109 alternative

Hello, I'm trying to create a for loop for a data set, I have a list of results in this data set and I want to take the 1st two add them together and divide by the mean of the 1st to, then do the same for the 3rd and 4th values in the list and so on and each time return a value for the calculation on each pair. Any help would be appreciated. Thanks, Al [[alternative HTML version

Hello, I have this data: Time AMP 0 0.2000000 10 0.1958350 20 0.2914560 40 0.6763628 60 0.8494534 90 0.9874526 120 1.0477692 where AMP is the concentration of this metabolite with time. If you plot the data, you can see that it could be fitted using a logistic regression. For this purpose, I used this code: AMP.nls <- nls(AMP~SSlogis(Time,Asym, xmid, scal), data

Hello, I developed an application that I would like to run in linux. It is developed part in VB6 part in C++ (using VS 6 and VS 2005) and uses Codejock's components for the UI. The installer is built with InnoSetup. I managed to install it on Suse 9.3 with Wine 0.9.17 but I got errors (cannot load DLL) while registering msxml, capicom, and one of my VS 2005 dlls. I fixed the msxml

Hi, everyone, Since quantile calculation has nine different methods in R, I wonder how I specify a method when calling the bwplot() in lattice. I couldn't find any information in the documentation. Thanks. -- Jun Shen PhD PK/PD Scientist BioPharma Services Millipore Corporation 15 Research Park Dr. St Charles, MO 63304 Direct: 636-720-1589 [[alternative HTML version deleted]]

hi, all I have two questions on using R console effectively (this is on Mac, not sure if it applies to win platform): First, I'd like to make the console font bigger, the default is too small for my screen. There is a "Show Fonts" from Format menu where you can adjust it, but it seems only for current session. Next time I start R, I have to redo everything. My question is, is there

Hi, everyone, Try the following command to see if you get TRUE or FALSE. I get FALSE on a unix platform but TRUE on Windows. Any comment? all(paste(Theoph[1],Theoph[2])==paste(Theoph[[1]],Theoph[[2]])) -- Jun Shen PhD PK/PD Scientist BioPharma Services Millipore Corporation 15 Research Park Dr. St Charles, MO 63304 Direct: 636-720-1589 [[alternative HTML version deleted]]

I try to load the package lattice by library(lattice). Although I have installed it through CRAN by internet as well as from a local directory using the zip-file, it can't be loaded into the editor. Can anyone give advice, please. I am using R 2.9.1 on Widows XP. Thanks a lot, Luba [[alternative HTML version deleted]]

hi, Previously my email very simple and less information. Let say I have raw data as below:- user_id website time 20 google 0930 21 yahoo 0935 20 facebook 1000 25 facebook 1015 61 google 0940 60 yahoo 1050 40 friendster 1115 61 google 1200 40

Dear R-users, The following question bothered me for the whole afternoon: how can one pass a string as the conditioning argument to subset? I tried plain mystr, eval(mystr), expression(mystr), etc... I don't to be able to find the correct syntax > foo <- data.frame(a=1:10,b=10:1,c=rep(1:2,5)) > mystr<-"c==1" > subset(foo,c==1) a b c 1 1 10 1 3 3 8 1 5 5 6

I have problems getting a histogram with relative frequencies on the y-axis. Here is an example data set: > a <- c(4.626, 4.627, 4.627, 4.628, 4.629, 4.629, 4.630, 4.631, 4.632, 4.632) > d = hist(a,freq=F) > d$density [1] 299.9999 100.0000 200.0000 100.0000 100.0000 200.0000 The obtained densities are given by counts/(total n * bin width), with bin width being

Thank you Gabor Grothendieck for your answer !! i have charge the pachage but now i have this massage  read.table("c:/TOTAL.xls",h=T) [1] ÐÏ.à.. <0 rows> (or 0-length row.names) Warning messages: 1: In read.table("c:/TOTAL.xls", h = T) :   incomplete final line found by readTableHeader on 'c:/TOTAL.xls' 2: In if (!header) rlabp <- FALSE :   the condition has

HI,   I have searched the R-help achive, trying to find a way to sort columns in a data frame according to a specific list, but not be able to find a general way to do it.   My data frame contains columns, with each one named from a sample ID, and each sample was measured in a specific time, generating one column of data, now I would like to sort the columns according to the time course, so I have

Hello, I have a data set with many variables, and often I want to run a given function, like summary() or cor() or lmer() etc. on many combinations of one or more than one of these variables. For every combination of variables I want to analyze I have been writing out the code by hand, but given that I want to run many different functions over dozens and dozens of variables combinations it is

Hi, and my apologies for the following very naive question : I would like to read a column of numbers in R and plot a histogram. eg : x<-read.table("txSTART"); y<-as.numeric(x); and I do obtain the error : Error: (list) object cannot be coerced to type 'double'. Please could you let me know the way to fix it. thanks, bogdan <r-help@r-project.org> [[alternative

Hi everyone, I am trying to accomplish a small task that is giving me quite a headache. I would like to automatically generate a series of matrices and give them successive names. Here is what I thought at first: t1<-matrix(0, nrow=250, ncol=1) for(i in 1:10){ t1[i]<-rnorm(250) } What I intended was that the loop would create 10 different matrices with a single column of 250

Hi, R-listers, Failed, when I tried to merge df1 and df2 by "codetot" in df1 and "codetoto" in df2. I want to know the reason and how to merge them together. Data frames and codes I have used were listed as followed. Thanks a lot in advance. df1: popcode codetot p3need BCPy01-01 BCPy01-01-1 100.0000 BCPy01-01 BCPy01-01-2 100.0000 BCPy01-01 BCPy01-01-3 100.0000 BCPy01-02