search for: tmpdat

Dear All, I have a problem in understanding how the interactions of 2 ways ANOVA work, because I get conflicting results from a t-test and an anova. For most of you my problem is very simple I am sure. I need an help with an example, looking at one table I am analyzing. The table is in attachment and can be imported in R by means of this command: scrd<-

...t$temp = as.factor(wheat$temp) ### convert to factor ### First, try to use nlsList wheat.list <- nlsList(Prop ~ SSlogis(end,phi1,phi2,phi3)| temp,pool=FALSE,data=wheat) ### ### next, try to use lm to estimate starting parameters. wheat.list = list() for (i in 1:length(levels(wheat$temp))){ tmpDat = wheat[wheat$temp == levels(wheat$temp)[i],] tmp.lm <- lm(Prop ~ end,data = tmpDat) tmp.nls <- nls(Prop ~ theta1 / (1 + exp(-(theta2 + theta3*end))), start = list(theta1 = .5,theta2=coef(tmp.lm)[1],theta3 = coef(tmp.lm)[2]), data = tmpDat,trace=TRUE) tmp2.nls <- nls(Prop...

Hi, I am running an analysis with an ordinal outcome and I need to run a test of the parallel regression assumption to determine if ordinal logistic regression is appropriate. I cannot find a function to conduct such a test. >From searching various message boards I have seen a few useRs ask this same question without a definitive answer - and I came across a thread that indicated there is no

...T) but it does not seem to work. trellis.par.set() does not have a component called "new" so I cannot use that to imitate what I would do with par. Does anyone know how to overlay two or more plots produced by plot.augPred? Thank you! Ramzi # Here is code that reproduces my problem: tmpdat = structure(list(ptno = c(4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 4L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 8L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 10L, 20L, 20L, 20L, 20L, 20L, 20L, 20L, 20L, 20L, 20L, 20L, 20L, 20L, 23L, 23L, 23L, 23L, 23L, 23L, 23L, 23L, 23L, 23L, 24L, 24L, 24L, 24L, 24L, 24L, 24L, 2...

Hello :) I'm just starting out with R and would appreciate your help with a couple of problems I am running into. I have used Sys.glob to get a list of all filenames having a particular file extension (in my case, *.txt) I would now like to use this list in the following manner: I would like to use each filename from the list and open that file into a tab separated matrix and proceed. How can

Hi, I am new in R and I have problem with reading this data file 0 TITLE Title 0 XLEGEND Legend -1 LABEL x 1 1 12 1 2 30 1 3 34 I want to read only lines starting with 1 (it indicates 1st plotting line) and create data set from second and third value on this row. Thank for advice Jakub Zlamal -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-help mailing list

...lsList > > wheat.list <- nlsList(Prop ~ SSlogis(end,phi1,phi2,phi3)| > > temp,pool=FALSE,data=wheat) ### > > > > ### next, try to use lm to estimate starting parameters. > > wheat.list = list() > > > > for (i in 1:length(levels(wheat$temp))){ > > tmpDat = wheat[wheat$temp == levels(wheat$temp)[i],] > > tmp.lm <- lm(Prop ~ end,data = tmpDat) > > tmp.nls <- nls(Prop ~ theta1 / (1 + exp(-(theta2 + theta3*end))), > > start = list(theta1 = .5,theta2=coef(tmp.lm)[1],theta3 = coef(tmp.lm)[2]), > > data = tmpD...