similar to: Plotting multiple curves on a plot

I have a simple X, Y data frame that I am trying to run regression analysis on. The linear regression looks great, but when I use lm(formula = y ~ poly(x, degree = 5)) I get the same coeffecients. So for example if I use degree =3 my formula would look like y = 4.2 x^3 + 3.2x^2 + 2.1x + 1.0 and my degree 5 would look like y = 6.5x^5+ 5.4x^4 + 4.2 x^3 + 3.2x^2 + 2.1x + 1.0, which doesn't make

PREFERED WAY OF DOING IT I have a data set of observations every second for a month long period, I want to extract the observations according to the date & time of another data frame ( the other data frame is in the same format). I want to do this to match these observations to my test observations (in the other data frame) which are done every 6 minutes. So basically im shrinking the

Hi R Users, I created interaction plots in ggplot2 and was trying to add output of two way ANOVA models, especially only interaction ( example treatment*control F(XX, XX) = xxx, p = xxx) into figures, but i was not able to add. Would you mind to help on how I can add information into each figure? I have attached the example data and the code that I used for this. dat<-structure(list(Sites

Dear R-helpers: Suppose I have a datafram called "test_frame" like this col1 col2 col3 col4 r1 x x x x r2 x x x x r3 x x x x .. x x x x rn x x x x I know I can get data of col3 by using test_frame[,3]. My question is, if I want to do this by specifying "col3" rather

Hello, I define a dataframe named "Price" first and there are other dataframes as well. Then I built an array named "code" containing all the dataframe name. So If I call code[1], it is "Price". Basicly, I would like to call the datafram with reference to the array"code".e,g, If I want to call dataframe "Prices", I use sth containing code[1]. It

Hello everyone, I would like to construct a datafram with the following command. eval(parse(text=paste("df=data.frame(", cmd, ")", sep=""))) But it comes out " Error in parse(file, n, text, prompt, srcfile, encoding) : syntax error, unexpected $undefined, expecting ',' in "df=data.frame(cbcDummy10to12 = 1,cbcForeWrld_Ret =

Hi All, Newbie question that i'm sure is easy, but i can't seem to apply properly I read in a datafram from a CSV file and i want to tell R that from coloum "n_0" to "n_32" the value "-1" is missing data i was looking at the is.na(xx) <- c(..,...,) idea but i can't seem to apply it properly, can anyone offer advice? On a side issue while i'm

Dear R Users,i am very new to R. I want your help on an issue regarding distance matrix and cluster analysis i had discharge data of 4 rivers(a,b,c,d) in 4 vectors each having 364 values > dput(qmu)structure(list(a = c(0.26, 0.25, 0.25, 0.25, 0.24, 0.23, 0.22, 0.21, 0.21, 0.21, 0.2, 0.19, 0.19, 0.19, 0.19, 0.18, 0.18, 0.18, 0.17, 0.17, 0.17, 0.17, 0.17, 0.17, 0.17, 0.17, 0.17, 0.17, 0.17, 0.17,

Hi! I am looking to choose a condom based on its pleasure score. I received some summarised data from 10 individuals: structure(list(Ramses = c(4, 4, 5, 5, 6, 3, 4, 4, 3, 4), Sheiks = c(5, 5, 6, 4, 7, 6, 4, 5, 6, 3), Trojans = c(7, 8, 7, 9, 6, 3, 2, 2, 2, 3), Unnamed = c(2, 1, 1, 3, 3, 4, 5, 4, 4, 3)), .Names = c("Ramses", "Sheiks", "Trojans", "Unnamed"),

Hello. I am trying to do a non-linear fit using the 'nls' command. The data that I'm using is as follows pH k 1 3.79 34.21 2 4.14 25.85 3 4.38 20.45 4 4.57 15.61 5 4.74 12.42 6 4.92 9.64 7 5.11 7.30 8 5.35 5.15 9 5.67 3.24 with a transformation of pH to H <- 10^-pH When using the nls command for a set of parameters - a, b and c, I receive two sets of errors: >

Hi i'm a french medical student, i have some data that i import from excel. My colomn of the datafram are the localisations of metastasis. If there is a metatsasis there is the symbol "_". i want to exclude the row without metastasis wich represent the NA data. so, i wrote this mela is the data fram mela1=ifelse(mela[,c(11:12,14:21,23,24)]=="_",1,0) # selection of the

Hi All I have a datafram which looks like this : CustomerID DietType 1 a 1 c 1 b 2 f 2 a 3 j 4 c 4 c 4 f And I would like to reshape this so I can

Hi list, I have the following entries, below, in today's log file (for yesterday, 10th May). I don't run the automated yum-updated and didn't run a yum update yesterday, and no packages were installed. Obviously the entries are old. I was wondering if anyone could offer an explanation? Thanks, Ned --------------------- yum Begin ------------------------ Packages

Hi, I am trying to do a lasso regression using the lars package with the following data (see attached): FastestTime WinPercentage PlacePercentage ShowPercentage BreakAverage FinishAverage Time7Average Time3Average Finish 116.90 0.14 0.14 0.29 4.43 3.29 117.56 117.77 5.00 116.23 0.29 0.43 0.14 6.14 2.14 116.84 116.80 2.00 116.41 0.00 0.14 0.29 5.71 3.71 117.24

Hi Allaisone, If you want a data frame as the output you will have to put up with a few NA values unless each Customer has the same number of diet types: a1df<-read.table(text="CustomerID DietType 1 a 1 c 1 b 2 f 2 a 3 j 4

Hi I am trying to build a regression model. My data looks like this: A B C D E 1 1 1 2.57 5 2 0 0 1.64 3 0 5 1 4.8 1 1 3 0 3.56 168 1 1 1 2.13 1 0 3 1 5 168 2 0 0 7.75 28 4 0 2 2.85 168 3 0 1 1.89 6 1 1 1 2.33 3 3 2 2 1.77 168 1 0 0 1.38 0.04 0 6 1 4.57 168 2 2 2 2.86 1 3 1 1 4.11 168 3 0 2 3 84 1 1 1 2.5 56 where E is a response variable and A, B, C,

Hi, I'm trying to use sapply to compute the min of several variables, each of them stored in data.frames, grouped as a list: Is it normal that mean() and min() produce different objects dimensions? > str(dats) List of 5 $ log20:'data.frame': 83 obs. of 5 variables: ..$ DATE : int [1:83] 2001081500 2001081512 2001081600 2001081612 2001081700 2001081712

Hi, I could try latest 7zip (4.57) on a triple boot machine (Debian SID i386, x86_64 and native Windows XP SP2), and found something surprising in the results : amd64 at 2400MHz linux32 -> 1814 linux64 -> 1946 win32 -> 2004 So performance is much better in Linux 64 bit than in 32 bit, while both are slower than native windows system. Is it normal, I mean does anyone else obtain

I believe you need to spend time with an R tutorial or two: a data frame (presumably the "table" data structure you describe) can *not* contain "blanks" -- all columns must be the same length, which means NA's are filled in as needed. Also, 8e^5 * 7e^4 = 5.6e^10, which almost certainly will not fit into any local version of R (maybe it would in some server version --

Ben, I would like to test the sulfur on the clover field, nitrogen on the clover field and then test for the presence of interaction. Sorry about the last email, seems it really screwed itself over, here it is again, hopefully nicer: Nitrogen(0) Nitrogen(20) Sulfur(0) 4.54 5.73 Sulfur(3) 4.64