search for: tms

...le 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, terms = terms) n.tms <- ncol(tms <- as.matrix(if (se) terms$fit else terms)) In this case terms is now a matrix of fitted values, with a single column for each term in the model or in the terms argument. Because the predict function returns a vector if called with a single ter...

...ence of using Least Square Fit with R I would appreciate your supervision and suggestion. I guess the returned coefficients of the oolynomial are: a = -1.3191398 b = 0.1233055 c = 0.9297401 Thank you very much in advance, Maura ###################################################### ## Main tms <- t(read.table("signal877cycle1.txt")) J <- ilogb(length(tms), base=2) + 1 y <- c(tms,rep(0,2^J - length(tms))) y.win <- tms.ext[1:17] ls.mat <- matrix(nrow=length(y.win),ncol=3,byrow=TRUE) dt <- 0.033 ls.mat[,1] <- 1 ls.mat[,2] <- seq(0,dt*(length(y.win)-1),dt) l...

Hello, Can anyone help me out with this problem? After performing logistic regressions and testing the significance with likelihood ratios, I have plotted the results using "termplot". I am wondering, how to get the names of my variables to appear on the x-axis rather than ascending numbers? I have used: termplot(fm, se=T, axes=FALSE, col.se="black") axis(1, 1:4,

..._NO_MATH_INLINES -g -O2 -c system.c -o system.lo gcc -I. -I../../src/include -I../../src/include -I/usr/local/include -DHAVE_CONFIG_H -D__NO_MATH_INLINES -g -O2 -c sys-unix.c -o sys-unix.lo sys-unix.c: In function `R_getProcTime':sys-unix.c:176: error: invalid use of undefined type `struct tms' sys-unix.c:177: error: invalid use of undefined type `struct tms' sys-unix.c:179: error: invalid use of undefined type `struct tms' sys-unix.c:180: error: invalid use of undefined type `struct tms' sys-unix.c: In function `Rf_InitTempDir': sys-unix.c:285: error: too many argume...

...xp prof workstation and maps their respective home folders as well as their respective shared folders on the gentoo server. For example, I have the following groups with their designated shares: Group Shared Dir ------------ ------------------ water /home/water tms /home/tms accounting /home/accounting and I have the users and groups User Group ------------ ------------------- bill tms , water, accounting jill water I am unable to get a client to join a domain unless I login and join...

...th, though it did until January. Could you provide an example? -thomas > Analysis: termplot starts with > > terms <- if (is.null(terms)) > predict(model, type = "terms", se = se) > else predict(model, type = "terms", se = se, terms = terms) > n.tms <- ncol(tms <- as.matrix(if (se) > terms$fit > else terms)) > > In this case terms is now a matrix of fitted values, with a single > column for each term in the model or in the terms argument. Because > the predict function returns a vector...

...end IF_2 } else { # ELSE_1 statement_11 } # end IF_1 It looks like R interpret does not like the above syntax. In fact in my script I have the following instructions: if (DonohoAplhaON){ tms <<- xx[,"sampamp"] #EXTRACT SIGNAL AMPLITUDE tmsLen <<- length(tms) J <<- ilogb(tmsLen, base=2) if (logb(tmsLen, base=2)%%2 > 0) { J <<- J + 1 } rm(xx) rawtms <- tms X <- Pre...

...> toothVC.lm <- lm(len ~ logdose, data=ToothGrowth, + subset=ToothGrowth$supp=="VC") > termplot(toothVC.lm) Error in xy.coords(x, y, xlabel, ylabel, log) : subscript out of bounds > traceback() 4: xy.coords(x, y, xlabel, ylabel, log) 3: plot.default(xx[oo], tms[oo, i], type = "l", xlab = xlabs[i], ylab = ylabs[i], xlim = xlims, ylim = ylims, main = main[i], col = col.term, lwd = lwd.term, ...) 2: plot(xx[oo], tms[oo, i], type = "l", xlab = xlabs[i], ylab = ylabs[i], xlim = xlims, ylim = ylims, main = main[i],...

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 the...

Hi; I have a question about speex on TMS320C54X: I have build and run the related project on Code Composer Studio IDE (ver. 3.1) and the male.snd (93.7KB) was its input file and malebits5x.dat (5.85KB) and maleout5x.snd (93.5KB) are output files. What is malebits5x.dat file? if it is encoded file, why I could'nt decode it? Al...

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 in do.call("fun", lapply(args, enquote)) : could not find function "fun". Any ideas? Best wishes Kristian ________________________________ Helmholtz Ze...

...uot;) stop(message = "Execution of move terminated.") } #end of if (sum(tmp) > 1) tmpos _ (1:length(search()))[tmp] } else tmpos _ pos #End of if (is.na(pos)) # Assign the correct directory name to dirN if it is not already assigned. tms _ search()[tmpos] dirN _ substring(tms, 6, nchar(tms)) #-------------------------------------------------- # The next if carries out the move, if this is possible. e.test _ NULL; for (i in 1:length(name)) e.test _ c(e.test, !exists(name[i], env=pos.to.env(tmpos)) ) if ( all(e.t...

...(.*)" L.DD <- grep(pat, L, value = TRUE) library(gsubfn) DD <- strapply(L.DD, pat, c, simplify = rbind) DDdf <- data.frame(gauge = as.numeric(DD[,1]), gauge_name = DD[,2]) both <- merge(DF, DDdf, by = "gauge", all.x = TRUE) dts <- as.character(both[,"date"]) tms <- as.character(both[,"time"]) date_time <- as.chron(paste(dts, tms), "%Y-%m-%d %H:%M") DF <- data.frame(Date=as.POSIXct(date_time), both) #change precip to numeric DF[,"precipitation"] <- as.numeric(as.character(DF[,"precipitation"])) precip.1 &...

Dear R-helpers, In my data set, I have a time variable 'RecordTime' whose class property is 'times'. When I list my data set, I see the values of RecordTime is like 10:20:30 in a 'h:m:s' format. Suppose I want to choose all the data after 10 o'clock, then use subset(data,RecordTime>10:20:30) just doesn't work. I noticed that all the values of RecordTime seem

...t; bug.data_ data.frame( x=c(1:3, NA), y=runif( 4)) test> lm.tpbug_ lm( y~x, data=bug.data) test> termplot( lm.tpbug) Error in xy.coords(x, y, xlabel, ylabel, log) : subscript out of bounds A simple fix seems to be add the commented line below, near the top of `termplot': n.tms <- ncol(tms <- as.matrix(if (se) terms$fit else terms)) mf <- model.frame(model) if (is.null(data)) data <- eval(model$call$data, envir) if (is.null(data)) data <- mf data <- data[ dimnames( tms)[[1]], ] # MVB nmt <- colnames(t...

...d) > formula(m) y ~ x + f3 In S+ formula(m) gives formula given to model.frame(), but in R you have to do the following get that formula: > formula(attr(m, "terms")) y ~ (x + f3)^2 Would it break anything to add to the top of formula.data.frame something like if (!is.null(tms <- attr(x, "terms"))) { return(formula(tms)) } so that formula() would retrieve the formula buried in model.frame's output? Bill Dunlap Spotfire, TIBCO Software wdunlap tibco.com

...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

...(.*)" L.DD <- grep(pat, L, value = TRUE) library(gsubfn) DD <- strapply(L.DD, pat, c, simplify = rbind) DDdf <- data.frame(gauge = as.numeric(DD[,1]), gauge_name = DD[,2]) both <- merge(DF, DDdf, by = "gauge", all.x = TRUE) dts <- as.character(both[,"date"]) tms <- as.character(both[,"time"]) date_time <- as.chron(paste(dts, tms), "%Y-%m-%d %H:%M") DF <- data.frame(date_time, both) library(ggplot2) #discharge if(input=="discharge"){ qplot(as.POSIXct(date_time), discharge, data=DF, geom="line")+facet_wrap(~ga...

I have a large dataset similar to this: ID time result A 1 5 A 2 2 A 3 1 A 4 1 A 5 1 A 6 2 A 7 3 A 8 4 B 1 3 B 2 2 B 3 4 B 4 6 B 5 8 I need to extract a number of features for each individual in it (identified by "ID"). These are: * The lowest result (the nadir) * The time of the nadir - but if the nadir level is present at >1 time point, I need the minimum and maximum time of nadir

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. Given the size of this data set, I was wondering if...