search for: esoph

...est to add a panel function to levelplot (or perhaps to an other 3d lattice function) which is able to translate the z values into the size of the rectangles. It could be used to display categorical data. I append the proposed code and two examples: - panel.catlevelplot() - example1.catlevelplot.esoph() - example2.catlevelplot.esoph() Wolfram Fischer #------ CODE -------------------------------------------------------------- panel.catlevelplot <- function (x, y, z, wx, wy, zcol, col.regions, subscripts , ... , z.factor.min = 0.02 # factor for z range expansion...

Un texte encapsul? et encod? dans un jeu de caract?res inconnu a ?t? nettoy?... Nom : non disponible URL : <https://stat.ethz.ch/pipermail/r-help/attachments/20080801/a1a7c3e9/attachment.pl>

This arose in testing [.terms and has me confused. data(esoph)?? # use a standard data set t0x <- terms(model.frame( ~ tobgp, data=esoph)) t1 <-? terms(model.frame(ncases ~ agegp + tobgp, data=esoph)) t1x <- (delete.response(t1))[-1] > all.equal(t0x, t1x) [1] TRUE # the above is wrong, because they actually are not the same > all.equal(at...

...{with default par()s !}: get(getOption("device"))() par(mfrow = c(2,2)) ## part of example(interaction.plot) _improved_ using with() : data(OrchardSprays) with(OrchardSprays, { interaction.plot(treatment, rowpos, decrease) interaction.plot(rowpos, treatment, decrease) }) data(esoph) with(esoph, { interaction.plot(agegp, alcgp, ncases/ncontrols) interaction.plot(agegp, tobgp, ncases/ncontrols, trace.label="tobacco", fixed=TRUE) }) ------------------------------- Something that should be easy to fix, but needs some care (and time). Volu...

Full_Name: Axel Benz Version: 1.7.1 OS: Windows Submission from: (NULL) (137.251.33.43) Hello, I guess many people will answer me again that this is a S language feature, but I am only a stupid computer scientist and I simply do not understand this logic, despite of reading a lot about S: > test field tuckey 4 Kreis2 -1 5 Kreis5 -2 9 Metall -3 17

On 05/04/2019 9:03 a.m., Therneau, Terry M., Ph.D. via R-devel wrote: > This arose in testing [.terms and has me confused. > > data(esoph)?? # use a standard data set > > t0x <- terms(model.frame( ~ tobgp, data=esoph)) > t1 <-? terms(model.frame(ncases ~ agegp + tobgp, data=esoph)) > t1x <- (delete.response(t1))[-1] > > > all.equal(t0x, t1x) > [1] TRUE > > # the above is wrong, because they...

Hello, I need to write out a data matrix as a STATA 7 file and this happens perfectly with write.dta(), except I cannot seem to export the labelnames to Stata. So far I have tried the following: # X is the data matrix that is to be exported attributes(X)$var.labels <- c("apple", "banana", "cat") write(X, filename = "text.dta", version = 7) When i

...es". > > Thanks for the pointer, > > Terry > > > > On 4/5/19 9:00 AM, Duncan Murdoch wrote: >> On 05/04/2019 9:03 a.m., Therneau, Terry M., Ph.D. via R-devel wrote: >>> This arose in testing [.terms and has me confused. >>> >>> data(esoph)?? # use a standard data set >>> >>> t0x <- terms(model.frame( ~ tobgp, data=esoph)) >>> t1 <-? terms(model.frame(ncases ~ agegp + tobgp, data=esoph)) >>> t1x <- (delete.response(t1))[-1] >>> >>> ? > all.equal(t0x, t1x) >>>...

...ke: *** [check-all] Error 2 I looked at "tests/reg-tests-1.Rout.fail"; it's 1427 lines long. The error given is > ## Comments: > > > ## PR 796 (aic in binomial models is often wrong) > ## > a1 <- glm(cbind(ncases, ncontrols) ~ agegp + tobgp * alcgp, + data = esoph, family = binomial())$aic > a1 [1] 236.9645 > a2 <- glm(ncases/(ncases+ncontrols) ~ agegp + tobgp * alcgp, + data = esoph, family = binomial(), weights=ncases+ncontrols)$aic > a2 [1] 236.9645 > stopifnot(a1 == a2) Error: a1 == a2 is not TRUE Execution halted I am running on an A...

...>>> >>> >>> >>> On 4/5/19 9:00 AM, Duncan Murdoch wrote: >>>> On 05/04/2019 9:03 a.m., Therneau, Terry M., Ph.D. via R-devel wrote: >>>>> This arose in testing [.terms and has me confused. >>>>> >>>>> data(esoph)?? # use a standard data set >>>>> >>>>> t0x <- terms(model.frame( ~ tobgp, data=esoph)) >>>>> t1 <-? terms(model.frame(ncases ~ agegp + tobgp, data=esoph)) >>>>> t1x <- (delete.response(t1))[-1] >>>>> >>>&...

I have a dataframe with the following variables: idnum area gender race etc. I would like to make a table that looks like area gender race M F B W A 1 4 5 3 5 1 2 6 7 4 6 3 etc. Basically, I want to make a single broad table with a number of sub-set tables. I have tried: cbind(table(area, gender), table(area, race)) But, when I do this, I lose the labels gender / race. This makes it

...; https://stat.ethz.ch/pipermail/r-devel/2017-January/073621.html In the mean time, I've spent several hours on the issue and also committed changes to R-devel "in two iterations". In the case there is a *Left* hand side part to xtabs() formula, see the help page example using 'esoph', it uses tapply(..., FUN = sum) and I now think there is a missing feature in tapply() there, which I am proposing to change. Look at a small example: > D2 <- data.frame(n = gl(3,4), L = gl(6,2, labels=LETTERS[1:6]), N=3)[-c(1,5), ]; xtabs(~., D2) , , N = 3 L n A B C D E F...

Dear List: Most of the analysis I do involves contingency tables. I am migrating to R from Stata and I have a number of questions about using contingency tables in R. I suspect that most of the things I want to do are very short R scripts that people on this list probably have. I wonder if you would be willing to share them. First, the presentation of tables by table() is not

...; https://stat.ethz.ch/pipermail/r-devel/2017-January/073621.html In the mean time, I've spent several hours on the issue and also committed changes to R-devel "in two iterations". In the case there is a *Left* hand side part to xtabs() formula, see the help page example using 'esoph', it uses tapply(..., FUN = sum) and I now think there is a missing feature in tapply() there, which I am proposing to change. Look at a small example: > D2 <- data.frame(n = gl(3,4), L = gl(6,2, labels=LETTERS[1:6]), N=3)[-c(1,5), ]; xtabs(~., D2) , , N = 3 L n A B C D E F...

...rmail/r-devel/2017-January/073621.html > > In the mean time, I've spent several hours on the issue > and also committed changes to R-devel "in two iterations". > > In the case there is a *Left* hand side part to xtabs() formula, > see the help page example using 'esoph', > it uses tapply(..., FUN = sum) and > I now think there is a missing feature in tapply() there, which > I am proposing to change. > > Look at a small example: > >> D2 <- data.frame(n = gl(3,4), L = gl(6,2, labels=LETTERS[1:6]), N=3)[-c(1,5), ]; xtabs(~., D2) &g...

...r2 (beavers) Body Temperature Series of Two Beavers cars Speed and Stopping Distances of Cars chickwts Chicken Weights by Feed Type co2 Mauna Loa Atmospheric CO2 Concentration discoveries Yearly Numbers of Important Discoveries esoph Smoking, Alcohol and (O)esophageal Cancer euro Conversion Rates of Euro Currencies euro.cross (euro) Conversion Rates of Euro Currencies eurodist Distances Between European Cities faithful Old Faithful Geyser Data fdeaths (UKL...

...titute(Y > 0, list(Y = formula[[2]]))), ##or subset = as.expression(bquote(.(formula[[2]]) > 0)), ##or subset = bquote(.(formula[[2]]) > 0), na.action = na.action) mf } ## never works tst(ncases ~ agegp + alcgp, data = esoph) traceback() #--> shows that inside model.frame.default # eval(substitute(subset, ...)) is called as well ---- Happy quizzing.. Martin Maechler, ETH Zurich

Good day My name is Jos? Oscar, I'm from Mexico and I have some questions about foreign in your write.foreig( ) function. We know that this function generates the inputs to be able to run them or execute them in another program like SPSS, SAS or Stata. In these cases, when creating an example file, a matrix and using the function to execute the .sps file directly from SPSS, I don't care or

...rmail/r-devel/2017-January/073621.html > > In the mean time, I've spent several hours on the issue > and also committed changes to R-devel "in two iterations". > > In the case there is a *Left* hand side part to xtabs() formula, > see the help page example using 'esoph', > it uses tapply(..., FUN = sum) and > I now think there is a missing feature in tapply() there, which > I am proposing to change. > > Look at a small example: > >> D2 <- data.frame(n = gl(3,4), L = gl(6,2, labels=LETTERS[1:6]), N=3)[-c(1,5), ]; xtabs(~., D2) &g...

...; >> In the mean time, I've spent several hours on the issue >> and also committed changes to R-devel "in two iterations". >> >> In the case there is a *Left* hand side part to xtabs() formula, >> see the help page example using 'esoph', >> it uses tapply(..., FUN = sum) and >> I now think there is a missing feature in tapply() there, which >> I am proposing to change. >> >> Look at a small example: >> >>> D2 <- data.frame(n = gl(3,4), L = gl(6,...