similar to: formula(model.frame(y~.^2, data=d)) does not return formula from terms attribute of the model.frame

Dear Martin, Since no one else has picked up on this, I?ll take a crack at it: The proposal is to define the S3 class of model-frame objects as c(?model.frame?, ?data.frame?) (not the formal class of these objects, even though this feature was coincidentally introduced in S4). That?s unlikely to do harm, since model frames would still ?inherit? data.frame methods. It's possible that some

I don't have a copy of SV4 (or SV3, where model.frame was introduced), but S+ 8.3 (based on SV4) puts the class "model.frame" on model.frame()'s return value but has no methods (in the default packages) for class "model.frame". Perhaps that is why R omitted the class. However, S+ 8.3's (and problably S's) formula.data.frame did look for a "terms"

The bug exists on R-1.9.0-alpha compiled the 10/3. Termplot has a problem if either the model only contains a single term or if asked to plot a single term. In addition there are problems with the option se = TRUE. Analysis: termplot starts with terms <- if (is.null(terms)) predict(model, type = "terms", se = se) else predict(model, type = "terms", se = se,

When formula() is applied to the output of model.frame() it ignores the formula in the model.frame's 'terms' attribute: > d <- data.frame(A=log(1:6), B=LETTERS[rep(1:2,c(2,4))], C=1/(1:6), D=rep(letters[25:26],c(4,2)), Y=1:6) > m0 <- model.frame(data=d, Y ~ A:B) > formula(m0) Y ~ A + B > `attributes<-`(terms(m0), value=NULL) Y ~ A:B This is in part

What is wrong in the following nested if-else statements: if (Condition_1) { # begin IF_1 statement_1 statement_2 statement_3 if (Condition_2) { # begin IF_2 a<- a +1 } # end IF_2 statement_4 statement_5 statement_6 statement_7 if (Condition_3) {

Hello list, Out of interest, I've been trying to compile R on WinXP with MinGW/MSYS. I've been following the instructions for compiling R on windows as closely as possible. Below is my configure statement, which seems to work okay except that no browser or pdf viewer is found, and html docs will not be compiled... I've tried enabling/disabling many different options, and I also add

Hi I try parallelising some code using the snow package and the following lines: cl <- makeSOCKcluster(8) pfunc <- function (x) (if(x <= (-th)) 1 else 0) ###correlation coefficient clusterExport(cl,c("pfunc","th")) cor.c.f <- parApply(cl,tms,c(1,2),FUN=pfunc) The parApply results in the error message: > cor.c.f <- parApply(cl,tms,c(1,2),FUN=pfunc) Error

note: whole function is below- I am sure I am doing something silly. when I use it like USGS(input="precipitation") it is choking on the precip.1 <- subset(DF, precipitation!="NA") b <- ddply(precip.1$precipitation, .(precip.1$gauge_name), cumsum) DF.precip <- precip.1 DF.precip$precipitation <- b$.data part, but runs fine outside of the function: days=7

Hi,all I only got the time data such as: tms<-c("19:30:23","18:39:10".....) I want to manipulate this time series data. For example, plus one second(or minute) or minus one second This data only has the time(h:m:s), without the date. I know that there are chron package, ISOPix class and the timeDate class, but all these class need the input of date. How can we

Look at the bottom of the message for my question #here is a little function that I wrote USGS <- function(input="discharge", days=7){ library(chron) library(gsubfn) #021973269 is the Waynesboro Gauge on the Savannah River Proper (SRS) #02102908 is the Flat Creek Gauge (ftbrfcms) #02133500 is the Drowning Creek (ftbrbmcm) #02341800 is the Upatoi Creek Near Columbus (ftbn) #02342500 is

Having some issues with: samba_dnsupdate Specifically: /usr/bin/nsupdate: cannot specify -g or -o, program not linked with GSS API Library I've looked through the script, and cannot find these options called. If anyone can point me to where they're called I'd appreciate it. Cheers, TMS III

Hi list, Sorry for the vague title, but here is the scenario. I?m writing an R package, let?s say, ?pkg1?, which contains 3 functions: f1, f2, f3. f2 and f3 are helper functions for f1, i.e. f1 calls f2 which in turn calls f3. f1 <- function(?) { ?. f2() ? } f2 <- function(?){ ? f3(?) ? } f3 <- function(...){ .... } Then, I wrote a new version of f3 and I want to

On Thu, 18 Mar 2004 k.hansen@biostat.ku.dk wrote: > The bug exists on R-1.9.0-alpha compiled the 10/3. > > Termplot has a problem if either the model only contains a single term > or if asked to plot a single term. In addition there are problems with > the option se = TRUE. I can't reproduce this in either R-devel or 1.8.1, and termplot hasn't changed since January. I do

Hello everyone, This IFELSE function call is not working properly. I do not receive an error message, but the actions are not executed conditional as I was hoping. Any assistance is appreciated. set.seed(12345) res1<-rbinom(10000,1,.1) rdata3<-transform(data.frame(res1),input1=rnorm(10000,50,10)) data3 #inducing correlation between res1 & input1

I will be out of town Tues 2/21 - Fri 2/24. If this is an emergency please call my cell at 818-400-6195, as I will only infrequently check my email during this period. Thanks, TMS III

Hi, I would like to use a TI dll used to communicate with TMS 320C6211 DSK with my own program written under Linux. Does anybody knows if it is possible (and how) to use (link) a Windows dll in a Linux program ?? Ciao, Riccardo Cerulli-Irelli

I'm running R 2.13 on Ubuntu 10.10 I have a data set which is comprised of character strings. site = readLines('http://www.census.gov/tiger/tms/gazetteer/zips.txt') dat <- c("01, 35004, AL, ACMAR, 86.51557, 33.584132, 6055, 0.001499") dat I want to loop through the data and construct a data frame with the zip code, state abbreviation, and city name in seperate columns.

I'm interested in plotting a y with an x factor as the combination of 2 factors and colour with respect to a third, which the code below does with interaction.plot(). However, this is because I redefine the x to be 1 factor. Is there a way of getting it to plot without redefining it, and ideally to not join up the lines BETWEEN levels a and b, but just join those between after and before for

Hello, In doing a series of non-linear estimations of a function which is a sum of a varying number of sinusoids, I would like to "autogenerate" the arguments needed by nls() depending on that number. For example, when there are two sinusoids: > nls( y ~ mu + A1 * cos(2*pi*f1*x - P1) + A2 * cos(2*pi*f2*x - P2), data = some.xy.data, start = list( mu=some.value0,

Having some difficulties determining which way to go here. Have an Ubuntu 10.2 server running S4Alpha11. The samba stuff is working well enough (though with a bit too much idle CPU consumption--haven't go there yet). However, I can't get proper responses for file system and user/group lookup. With AD in Samba 3 nsswitch.conf was passwd: files winbind group: files