similar to: Extracting variable names of final model in stepAIC

Hi. I uses rpart to build a regression tree. Y is continuous. Now, I try to predict on a new set of data. In the new set of data, one of my x (call Incoterm, a factor) has a new level. I wonder why the error below appears as the guide says "For factor predictors, if an observation contains a level not used to grow the tree, it is left at the deepest possible node and

Can I say that RPART is a modified algo of CART and PARTY a modified of CHAID? Thanks. ---- Chua Siang Li Consultant - Operations Research Acceval Pte Ltd Tel: 6297 8740 Email: siang.li.chua at acceval-intl.com Website: www.acceval-intl.com This message and any attachments (the "message"...{{dropped:12}}

Hello there. I have fitted a rpart model. > rpartModel <- rpart(y~., data=data.frame(y=y,x=x),method="class", ....) and can use rpart$where to find out the terminal nodes that each observations belongs. Now, I have a set of new data and used predict.rpart which seems to give only the predicted value with no information similar to rpart$where. May I know how

Hi. I need some advises on how to use R to find pi (i is the index) with the following objective function and constraint: max (sum i)[ f(ai, bi, pi) * g(ci, di, pi) * Di ] s.t. (sum i)[ f(ai, bi, pi) * Di * pi] / (sum i)[ f(ai, bi, pi) * Di ] <= constant f and g are diffentiable. So, I am thinking of optim with method = "BFGS"? But wonder how to include the

Hello there. Is there any function in R that can do cluster on a set of data that has both categorical and numerical variables? thanks. siangli

Dear community, I'm having this silly problem. I've a linear model. After fixing it, I wanted to know which data had studentized residuals larger than 3, so i tried this: d1 <- cooks.distance(lmmodel) r <- sqrt(abs(rstandard(lmmodel))) rstu <- abs(rstudent(lmmodel)) a <- cbind( mydata, d1, r,rstu) alargerthan3 <- a[rstu >3, ] And suddenly a[rstu >3, ] has

Dear R-helpers, I would appreicate if someone can help me on the transfer parameter in ARIMAX and also see what I am doing is correct. I am using ARIMAX with 2 Exogeneous Variables and 10 years data are as follows: DepVar Period, depVar, IndepVar1 Period, indepVar1, IndepVar2 Period, indepVar2 Jan 1998,708,Jan 1998,495,Jan 1998,245.490 Feb 1998,670,Feb 1998,421.25,Feb 1998,288.170 Mar

Hi, A following function was kindly provided by GGally?s maintainer, Barret Schloerke. function(data, mapping, ...) { p <- ggplot(data = data, mapping = mapping) + geom_point(color = I("blue")) + geom_smooth(method = "lm", color = I("black"), ...) + theme_blank() + theme(panel.border=element_rect(fill=NA, linetype =

Dear R-helpers. Would appreciate if someone can explain the warning messages below, after auto.arima. I couldn't find any clue in the archived help. Also, how do I retrieve the AICs of each tried model in auto.arima? The purposes are (1) to output to a text file, and (2) to find the 2nd best model by finding 2nd lowest AIC instead of eyeballing thru the value at the console

Hi, I am trying to build a binary of Rquantlib package 0.3.5 under Window 7 64bits. I am using mingw gcc version 4.5.1 (tdm64-1) which is a experimental version of 64 bits ming gcc. (http://tdm-gcc.tdragon.net/). I have successfully compiled boost 1.45 and Quantlib 1.0.1 with this gcc. And then when I tried to compile Rquantlib 0.3.5, it reports some error. The first error is like below.

Dear Help: I am trying to follow Professor Bates' recommendation, quoted by Professor Crawley in The R Book, p629, to determine whether I should model data using the 'plain old' lm function or the mixed model function lmer by using the syntax anova(lmModel,lmerModel). Apparently I've not understood the recommendation or the proper likelihood ratio test in question (or both) for I

Dear Community, I want to identify outliers in my data. I don't know how to use identify command in the plots obtained. I've gone through help files and use mahalanobis example for my purpose: NormalMultivarianteComparefunc <- function(x) { Sx <- cov(x) D2 <- mahalanobis(x, colMeans(x), Sx) plot(density(D2, bw=.5), main="Squared Mahalanobis distances, n=nrow(x),

hi i have switch to using oddcast under window platform and now using icecastwin32 as server in this case is there any possible way to let icecastwin32 to stream Wma? pls advise thanks you -----Original Message----- From: Geoff Shang [mailto:geoff@hitsandpieces.net] Sent: Wednesday, February 08, 2006 1:06 PM To: Neo Jiun Siang - R&D Cc: icecast@xiph.org; njsalan@hotmail.com Subject: RE:

thanks alot Geoff by the way is there any possible ways to add plug-in to icecast2 server to stream Wma? regardless of the license problem -----Original Message----- From: Geoff Shang [mailto:geoff@hitsandpieces.net] Sent: Tuesday, February 07, 2006 8:08 PM To: Neo Jiun Siang - R&D Cc: icecast@xiph.org; njsalan@hotmail.com Subject: RE: [Icecast] guideline for plug-in to icecast2 server Neo

Hello I want to compare the slope (let's say 'b') of a linear model obtained with lm() to a theoretical value (let's say 'th'). To do so, I think I should compute a 't value' using something like : (b - 'th')/standard.deviation(b) and then look at the p-value of this computed t. I don't understand how to do this in a simple way, just using lm()

hi actually im trying to add feature to ices2 to able to stream mix file format such as ogg and mp3 in the same playlist.m3u without extra decoding and encoding that function like ezstream. the reason i dont want to use ezstream as source client becoz after going thru decoding and encoding, the audio quality drops. Therefore im trying to find out is there any easy way out to add plug-in to

Helo there. I am using kmeans of base package to cluster my customers. As the results of kmeans is dependent on the initial centroid, may I know: 1) how can we specify the centroid in the R function? (I don't want random starting pt) 2) how to determine the optimal (if not, a good) centroid to start with? (I am not after the fixed seed solution as it only ensure that the

hi all i would like to find out is that any ways to add plugin to icecast2 and ice2 to stream mp3 .Is there a list of gulidline for adding feature to present icecast2 and ices2. i will a appretiated for the help from icecast memeber and developer thank you .

Dear helpeRs- I'm using a for loop to create a series of models. I'm trying to assign a name to each model created, using the loop index. The loop gets stuck at the name of the model, giving the error "target of assignment expands to non-language object". The linear model runs without error; only the name is problematic. Here is the current loop syntax. The use of dat

I've just recently converted my Windows pc to Mandrake Linux & I've got everything I need working except for this software - Equis Metastock. I've managed to install it. However, on execution, it dies with the following msg: "A CreateFont error occurred while loading the Symbols font: C:\Program Files\Equis\MetaStock\equis.ttf". On the wine console, the following msg