search for: juding

Hi, I am using function mars() in package mda to find knots in a whole bunch of predictor variables. I hope to be able to replicate all or some of the basis functions that the MARS software from Salford Systems creates. When I ran mars() on a small dataset, I was able to get the knots. However, when I tried running mars() on a larger dataset (145 predictor variables), for a different

R-Devel-ers: I have an S4 method that simply formats an object: setGeneric("formatMe", function(x) standardGeneric("formatMe")) setMethod("formatMe", "ANY", function(x) format(x)) If I issue the above in an R session, then define an S4 class with its own format method, I get the desired result: > setClass("A",contains="numeric")

Here is the website link: https://jobs.stjude.org/css_external/CSSPage_JobDetail.ASP?T=20141205131210& The position is particularly responsible for GPU (CUDA) programming to translate statistical algorithms into R CUDA package, as well as the maintenance of the R CUDA packages. The position calls for willingness to learn and adapt new methods and technologies for parallel and distributed

Hi All, I am running quantile regression in a "for loop" starting with 1 variable and adding a variable at a time reaching a maximum of 20 variables. I get the following warning messages after my "for" loop runs. Should I be concerned about these messages? I am building predictive models and am not interested in inference. Warning messages: 1: In

Hi, I am trying to read a SAS version 9.1.3 SAS dataset into R (to preserve the SAS labels), but am unable to do so (I have read in a CSV version). I first created a transport file using the SAS code: libname ces2 'D:\CES Analysis\Data'; filename transp 'D:\CES Analysis\Data\fadata.xpt'; /* create a transport file - R cannot read file created by proc cport */ proc

I think I'm missing something tonight in the usage of 'legend': plot(0, type="n") legend(locator(1), month.abb[1:5], pch=15, col=1:5) gives me something similar to what I want. But legend(locator(1), month.abb[1:5], fill=T, col=1:5) gives me 5 black boxes. What am I doing wrong? Thank you, > version _ platform i686-pc-linux-gnu arch i686 os

Hi, Is there a way to specify a vector of colors for the tick annotation in a call to axis(), to achieve the x-axis here? <<Rplot003.png>> Thanks, Mark Wilkinson Informatics Analyst St. Jude Children's Research Hospital Department of Pharmaceutical Sciences The opinions expressed here are my own and do not necessarily represent those of St. Jude Children's Research

Hi, I posted a while ago about 'overlap' problems using png/jpeg. If what Patrick Connolly suggests is truly happening, I think the following may be related. My new problem is with the following code (the overlap is still there if I use png() instead of postscript(), compounded by the difficulty described below): tmp <- matrix(runif(16000), nrow=16)

When I execute the following code, it works just like I want it to: three pages of nine (or fewer) plots. However, when I execute the code with the first and last lines uncommented, I get three pages (files), but the 2nd & 3rd pages have overlapping plots. It's like a new page wasn't created. I'm pretty sure I've either misplaced or left out a crucial call to some

Sorry for these multiple postings. I solved the problem using na.omit() to drop records with missing values for the time being. I will worry about imputation, etc. later. I calculated the sum of squared errors for 3 models, linear regression, neural networks, and support vector machines. This is the first run. Without doing any parameter tuning on the SVM or playing around with the number of

Hi mates, I have a problem. I am new to R and want to conduct the Fama/French asset pricing test. As I am from Germany, I cannot use the already computed factors from French's website, but need to compute them myself. So I have to sort a number of stocks into different portfolios using one factor, then subdivide these portfolios into several more using another factor, then compute portfolio

Since doing "yum update" this morning, an update that installed several packages, the updater shows me 3 packages to install: hpijs-1:3.12.4-4.el6_4.1 (x86_64) hplip-common-3.12.4-4.el6_4.1 (x86_64) hplip-libs-3.12.4-4.el6_4.1 (x86_64) but when I attempt to actually run the update I get: hplip-libs-3.12.4-4.el6.i686 requires libsane.so.1 hplip-libs-3.12.4-4.el6.i686 requires

As a follow-up to my email below: The input data frame to nnet() has dimensions: > dim(coreaff.trn.nn) [1] 5088 8 And the predictions from the neural network (35 records are dropped - see email below for more details) has dimensions: > pred <- predict(coreaff.nn1) > dim(pred) [1] 5053 1 So, the following line of R code does not work as the dimensions are

Hi, I am using the rpart package to build a classification tree. I did manage to build a tree with data on a previous project. However, when attampting to build a tree on a project I am working on, I seem to be getting the error shown below: > nhg3.rp <- rpart(profitresp ~., nhg3, method="class") > plot(nhg3.rp, branch=0.4, uniform=T); text(nhg3.rp, digits=3) Error in

I need to use the Bhattacharyya distance metric to determine population separation. Has anyone written a Bhattacharyya distance metric function in R? -- View this message in context: http://old.nabble.com/Bhattacharyya-distance-metric-tp26221259p26221259.html Sent from the R help mailing list archive at Nabble.com.

Hello, it wanted to know how I can extract of a dates frame the values peaks according to an interval that I establish. For example if dates are: 1 23 2 4 3 56 4 7 5 99 6 33 extract the date i wanted to divide into intervals of 2 an d to take alone the numbers 23, 56 and 99 of those 3 intervals. Thanks Ruben

Good morning, I am merging two datasets and I would like to save the non-matching rows in a separate file. The problem is how to retrieve the non-matching rows in R. Example: DATASET A code nomi A1 Franco A2 Mario A3 Andrea A4 Sandro A5 Luca DATASET B code book A1 Guerra e Pace A1 Storia di Roma A2 La coscienza di Zeno A4

Hi All, I am doing econometric modeling of panel data (fixed effects). We currently use Eviews to do this, but I have discovered a bug in Eviews 7 and am exploring the use of R to build panel data models / longitudinal data models. I looked at the plm package but do not see how I can incorporate AR terms in the model using the plm package. I have an Eviews model with two AR terms, AR(1) and

Hi, I use nnet for my classification problem and have a problem concerning the calculation of the final value for my validation data.(nnet only calculates the final value for the training data). I made my own final value formula (for the training data I get the same value as nnet): # prob-matrix pmatrix <- cat*fittedValues tmp <- rowSums(pmatrix) # -log likelihood

Hi, I apologize if this has previously been posted. I've just subscribed to the R-help digest. I'm writing a plotting function that uses layout() to plot several different plots on the same device. This function uses plot(), image(), and a custom function that uses text(). Each cell of the layout needs different par() parameters, so what I'd like to do is pass them as lists: