similar to: bug

R is not for Relational, but could be... I was musing on data.frames when I decided what I needed was a SQL like 'select' statement. Before I get involved in writing something like that, has anyone else written anything that allows a relational like algebra on data.frames? Thanks, John Carter EMail: ece at dwaf-hri.pwv.gov.za Telephone : 27-12-808-0374x194 Fax:-

Hi, I have something like this: > x <- 1:10 > y2 <- 30+5*x+rnorm(x,sd=3) > y <- c(y1,y2) > x <- c(x,x) > plot(x,y) > x <- 1:10 > y1 <- 1+5*x+rnorm(x,sd=2) > y2 <- 30+5*x+rnorm(x,sd=5) > y <- c(y1,y2) > x <- c(x,x) > f <- factor(rep(c("a","b"),c(10,10))) > m <- lm(y~x+f) > anova(m) Analysis of Variance

Dear all, I have a question regarding the output of the coxph function. What does the 'max possible' exactly mean in the output below? Many thanks. coef exp(coef) se(coef) robust se z Pr(>|z|) smocc_zyban -0.4384 0.6451 0.8667 0.9473 -0.463 0.644 self 1.1857 3.2728 0.1405 0.1443 8.216 2.22e-16 *** smocc_zyban:self

Hey I'm not really sure what I should put on here, but I am having trouble with my R code. I am trying to get the p-values, R^2s etc for a number of different groups of variables that are all in one dataset. This is the code: #Stand counter st<-1 #Collections stands<-numeric(67) slopes<-numeric(67) intercepts<-numeric(67) mses<-numeric(67) rsquares<-numeric(67)

Amigos de la erre. He creado mi primer bucle con for para entrenar unos modelos con GAM. La respuesta es quasipoisson porque estoy trabajando con densidades de peces. Sin embargo, tengo un problema, no se muy bien como añadir el offset a la formula siguiente cuando creo el bucle. GAM.A1 <-gam ((DYO)~s(DMA,k=4)+ s(WOD,k=4)+s(CIN,k=4)+s(DRA,k=4)+s(DBR,k=4)

Hi, I am trying to run the command "forward.sel.par," however I receive the error message: "Error: could not find function 'simpleRDA2'." I have the vegan library loaded. The documentation on "varpart" has not helped me to understand why I cannot call this function. Maybe I am missing something obvious because I am still an 'R' novice. Below is a

I think you've badly misinterpreted the purpose of the R listserv with this request: https://www.stat.math.ethz.ch/mailman/listinfo/r-help says "The `main' R mailing list, for announcements about the development of R and the availability of new code, questions and answers about problems and solutions using R, enhancements and patches to the source code and documentation of R,

In R-3.6.0 autoprinting was changed so that print methods for the storage modes are not called when there is no explicit class attribute. E.g., % R-3.6.0 --vanilla --quiet > print.function <- function(x, ...) { cat("Function with argument list "); cat(sep="\n ", head(deparse(args(x)), -1)); invisible(x) } > f <- function(x, ...) { sum( x * seq_along(x) ) }

Letting a user supply the autoprint function would be nice also. In a way you can already do that, using addTaskCallback(), but that doesn't let you suppress the standard autoprinting. Having the default autoprinting do its own style of method dispatch doesn't seem right. Bill Dunlap TIBCO Software wdunlap tibco.com On Tue, May 21, 2019 at 10:50 AM Lionel Henry <lionel at

It also is a problem with storage.modes "integer" and "complex": 3.6.0> print.integer <- function(x,...) "integer vector" 3.6.0> 1:10 [1] 1 2 3 4 5 6 7 8 9 10 3.6.0> print(1:10) [1] "integer vector" 3.6.0> 3.6.0> print.complex <- function(x, ...) "complex vector" 3.6.0> 1+2i [1] 1+2i 3.6.0>

Dear all, I have a question concerning the p-value. When running coxph I get a p-value = 0. :confused: Can this be true? Why aren?t there decimal points? Is there a way to find out the exact p-value? Here is the output: ---------------------------------------------------------------------------------------------------- Call: coxph(formula = Surv(start, stop, status) ~ Albumin_gproL, data = dial,

Hi, All: I'm having trouble getting deparse(substitute(x)) inside print.generic to consistently I'm having trouble getting a print.something to work consistently. Consider the following toy example: # Define an object of class 'dum' k <- 1 class(k) <- 'dum' str(k) # as expected # Define print.dum print.dum <- function(x, ...)

Hi all, I am using a 'tricked' Cox Hazard regression model for discrete choice habitat modelling. However, I'm having a hard time understanding the meaning of the first line the following part of the summary() output: Rsquare= 0.307 (max possible= 0.475 ) Likelihood ratio test= 91.8 on 12 df, p=2.23e-14 Wald test = 26.3 on 12 df, p=0.00977 Score (logrank) test = 58.6 on 12 df,

Hi, I ran the example on pp. 799-800 from Machael Crawley's "The R Book" using package survival v. 2.36-5, R 2.13.0 and RStudio 0.94.83. The model is a Cox's Proportional Hazards model. The result was quite different compared to the R Book. I have compared my code to the code in the book but can not find any differences in the function call. My results are attached as well as a

R 2.8.1 Windows XP I am trying to plot the results of a coxph using plot(survfit()). The plot should, I believe, show two lines one for survival in each of two treatment (Drug) groups, however my plot shows only one line. What am I doing wrong? My code is reproduced below, my figure is attached to this EMail message. John > #Create simple survival object >

Hi. I am currently trying to run some Spearman correlations, and have encountered two issues. 1) When using cor.test() with a variable that includes ties, I get the "Cannot compute exact p-values with ties" error. I have read that this function now uses an asymptotic formula that allows for ties, so do not understand why I am getting this error. (I am running version 2.4.0.) I

Hi, Peter et al.: On 2017-01-09 4:24 AM, peter dalgaard wrote: > On 09 Jan 2017, at 10:53 , Spencer Graves <spencer.graves at prodsyse.com> wrote: > >> # Define an object of class 'dum' >> k <- 1 >> class(k) <- 'dum' >> str(k) # as expected >> >> # Define print.dum >> print.dum <- function(x, ...) >>

Hello, I am working on a variable selection problem and I wonder whether there is some function or package in R works similar to the 'PROC REG' in SAS? Thank you. Some facts about 'PROC REG': PROC REG in SAS first composes a crossproducts matrix. The matrix can be calculated from input data, reformed from an input correlation matrix, or read in from an SSCP data set. For each

This is a simple R program that I have been trying to run. I keep running into the "singular matrix" error. I end up with no sensible results. Can anyone suggest any changes or a way around this? I am a total rookie when working with R. Thanks, Rasika > library(survival) Loading required package: splines > args(coxph) function (formula, data, weights, subset, na.action, init,

I am trying to use the Windows version of R to do the following: - Have a Delphi front-end which gets the user input - Delphi will produce a R script, do a system call to R to calculate the statistics. Is there a better way? Maybe a library one can use in Delphi or Visual C? - How do one see a graph? A graph is plotted in a window which IMMEDIATELY closes afterwords so that one can't