similar to: rmaxwell in C

Dear all; A simple? question. I'm having a problem with a math expression in the legend of a plot and I haven't found the way to get this to work, so any help will be appreciate. Basically I want to include in the plot is the R-squared and its numerical value, so I tried this: R2c<-0.82879 # R-squared of calibration model plot(1:10,1:10) legend("topleft",

I have a column in a data frame that has a class of "Date" and a mode of "numeric". When I: max(df$Date) My output stays in Date format, i.e. "2006-01-03". However, when I run the following statment: tapply(df$Date, df$SomeFactor, max) my output looks like this: 9129 9493 9861 10226 10591 10956 11320 11687 12052 12417 The returned object is of

Hi ALL, Recently I confronted with a error in using package ProDenICA, when let k<p (exactly the example code with k=1, http://finzi.psych.upenn.edu/R/library/ProDenICA/html/ProDenICA.html), Error in solve.default(V, W) : 'b' must be compatible with 'a' I have examined the code of function ProDenICA(), but could not find the reason. And how large is the sample size N

Hi All! I'm new to R and I need to know is it possible for R to generate C/C++ source code, Java byte code or native Win32 DLL like MatLab? --- WBR, Vyacheslav.

Dear all, I have two questions. First, I am wondering whether we are able to use R to generate an exe file, or sth that can be executable outside R? Second, I am wondering whether read.csv can read a csv file with size of 300-400 gigabytes? Thank you very much! Min [[alternative HTML version deleted]]

The symptom triggering this email is that an S4 summary method sometimes refuses to be invoked, even when a package is explicitly loaded, if the first load of the package is implicit. It may or may not be specific to 'summary' methods and/or the 'VGAM' package. I've sent to R-devel because (i) it looks like some kind of bug to me, but I'm not sure; (ii) it's not

I just discovered that when the VGAM package (not on CRAN) is loaded, glm() doesn't work. This is because VGAM defines a family function() which gets found by glm() in place of the family function from stats. Then VGAM:::family returns an object which doesn't have a $family component, (it has a component $vfamily). I thought namespaces should protect us from this happening? Kjetil --

I would like to turn off all the messages during library(aPackage) or require(aPackage) I tried different commands: invisible, capture.output, sink but none of them is working. For example, loading VGAM, gives a lot of unnecessary messages: > library(VGAM) Attaching package: 'VGAM' The following object(s) are masked from package:splines : bs The

Dear R-help gurus (and T.Yee, the VGAM maintainer) - I've been banging my head against the keyboard for too long now, hopefully someone can pick up on the errors of my ways... I am trying to use optim to fit a zero-inflated negative binomial distribution. No matter what I try I can't get optim to recognize my initial parameters. I think the problem is that dnbinom allows either

Dear all , I am trying to install a package by name 'netmodels' keep on getting the following error : > install.packages("netmodels") Warning in install.packages("netmodels") : argument 'lib' is missing: using '/usr/lib64/R/library' also installing the dependency ‘VGAM’ trying URL

Hi, I was going through the examples in smartpred. It seems there's an unstated dependency on the fitted() function in package flexmix. n = 20 set.seed(86) x = sort(runif(n)) y = sort(runif(n)) library(splines) fit = lm(y ~ ns(x, df=5)) plot(x, y) lines(x, fitted(fit)) # won't work w/o prior loading of the flexmix package. newx = seq(0, 1, len=n) points(newx, predict(fit,

Hi Folks, I wonderif someone who is familiar with the details of vglm in the VGAM package can assist me. I'm new to using it, and there doesn;t seem much in the documentation that's relevant to the question below. Say I have a vector x of 0/1 responses and another vector y of 0/1 responses, these in fact being a bivariate set of 0/1 responses equivalent to cbind(x,y). E.g.

Hi! R users I would like to ask you where could we find the VGAM Package. I don't find it in the list of packages. Thak you for your help Lassana KOITA Etudes de S?curit? et d'Exploitation a?roportuaires / Safety Study & Statistical analysis Service Technique de l'Aviation Civile (STAC) / Civil Aviation Technical Department Direction G?n?rale de l'Aviation Civile (DGAC) /

Dear R developers: I am trying to fit a PPO model using vglm from the library VGAM, and get an error while executing the code. Here is the data, code, and error: Data: first row is the column names. a = age, and 1,2,3, 4 and 5 are condition grades. a 1 2 3 4 5 1 1 0 0 0 0 2 84 2 7 10 2 3 16 0 6 6 2 4 13 0 3 4 0 5 0 0 0 1 0 Library(VGAM)

Dear, When i multiply the y of a regression by 10, I would expect that the coefficient would be multiply by 10 and the z value to stay constant. Here some reproducible code to support the case. *Ex 1* library(mvtnorm) library(VGAM) set.seed(1) x=rmvnorm(1000,sigma=matrix(c(1,0.75,0.75,1),2,2))

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/20110905/43ff838a/attachment.pl>

Anyone know how to get p-values for the t-values from the coefficients produced in vglm? Attached is the code and output ? see comment added to output to show where I need p-values + print(paste("********** Using VGAM function gamma2 **********")) + modl2<- vglm(MidPoint~Count,gamma2,data=modl.subset,trace=TRUE,crit="c") + print(coef(modl2,matrix=TRUE))

Hello R users, I have a puzzle with the VGAM package, on my first excursion into generalized additive models, in that this very nice package seems to want to do either more or less than what I want. Precisely, I have a 4-component outcome, y, and am fitting multinomial logistic regression with one predictor x. What I would like to find out is, is there a single nonlinear function f(x) which acts

pkgA's NAMESPACE has importFrom(graphics, plot) exportClasses("A") exportMethods("plot") R/foo.R has setClass("A") setMethod("plot", "A", function(x, y, ...) {}) During R CMD INSTALL pkgA_1.0.tar.gz we are told ** preparing package for lazy loading Creating a generic function for 'plot' from package

Hello, I have been using mgcv to run GAM hurdle models, analyzing presence/absence data with GAM logistic regressions, and then analyzing the data conditional on presence (e.g. without samples with no zeros) with GAMs with a negative binomial distribution. It occurs to me that using the negative binomial distribution on data with no zeros is not right, as the negative binomial allows zeros.