similar to: Question about mars() -function

Can you help me interpret the output I get with "polymars"? I'd like to find a regression model to predict 2 dependent variables (then called Y1 and Y2) with 2 independent variables (then called X1 and X2). Here is the output: polymars(responses = data[, 13:14], predictors = data[, 2:3]) Model fitting 0/1 size RSS 1 RSS 2 GCV 1 1 1 0.2486744

Hi, I have a set of covariance matrices but not the original data. I want to carry out some exploratory factor analysis. So, I am trying to construct a covariance matrix list as the input for factanal. I can construct a list which includes the cov, the centers, and the n.obs. But it doesn't work. I get an error that says "Error in sqrt(diag(cv)) : Non-numeric argument to mathematical

Hi, I have the following problem (library Hmisc loaded, iris data loaded, R Version 1.7.0 (2003-04-16), packages updated, running on a linux Debian i386): > summary(Species~Sepal.Length,method="reverse")->a > a Descriptive Statistics by Species +------------+-----------------+-----------------+-----------------+ | |setosa |versicolor |virginica

Hi, I am having a problem figuring out how to extract a subset of rows. I have a list with 68 similar data.frames. Each data.frame is 500 rows by 5 columns. I want to take one row from each data.frame based upon the data in a particular column (i.e. it matches a symbol). For example: > str(database) List of 68 $ X2011.01.11:'data.frame': 500 obs. of 5 variables: ..$ Symbol

Dear all Suppose the object below: > require(Hmisc) > require(plyr) > x <- dlply(iris, .(Species), describe) How can I print the object without displaying the attributes? I inspected ?print and ?print.default with no luck. > x $setosa x[, "Sepal.Length"] n missing unique Mean .05 .10 .25 .50 .75 50 0 15 5.006 4.40 4.59

Hi all, I am using "lm" to fit some anova factor models with interactions. The default setting for my unordered factors is "treatment". I understand the resultant "lm" coefficients for one factors, but when it comes to the interaction term, I got confused. > options()$contrasts unordered ordered "contr.treatment"

Hello everybody, I'm running a panel regression with dummy variables controlling for years. While this works fine for most of my models, in some cases I get the error message "In sqrt(diag(object$vcov)) : NaNs wurden erzeugt" and missing results for some or all of my year dummies: Estimate Std. Error z-value Pr(>|z|) year2000

http://biomet.oxfordjournals.org/cgi/reprint/86/3/677 biometrika1999 http://biomet.oxfordjournals.org/cgi/reprint/94/4/1006 biometrika2000 Hi All: I just want to try some luck. I am currenly working on my project,one part of my project is to reanalysis the kenward cattle data by using the method in Mohsen's paper,but I found I really can get the same or close output as he did,so,any

Dear list, I'm trying to learn how to use nlme to be able to fit ad analyse mixed-model logistic regressions. In order to keep things simple, I started by the simplest possible model : a one (fixed-effect ...) continuous variable. This problem is, of course, solved by glm, but I wanted to look at a "hand-made" nls fit, in order to be able to "generalize" to nlme

Dear list members, I am trying to estimate parameters of the AR(1)-GARCH(1,1) model. I have one additional dummy variable for the AR(1) part. First I wanted to do it using garchFit function (everything would be then estimated in one step) however in the fGarch library I didn't find a way to include an additional variable. That would be the formula but, as said, I think it is impossible to add

Dear all: I converted the columns (Baci, Meti, Fungii, Protozoai) into integers (using excel) and then imported the data (.txt) into R. Interestingly, the other three variables were loaded as INT, but the 'Baci' one continued as Num. I imported the data using the following command line: X <- read.delim(file.choose(), header = TRUE, dec =

Hi everyone, I'm sorry if this turns out to be more a statistical question than one specifically about R - but would greatly appreciate your advice anyway. I've been using a logistic regression model to look at the relationship between a binary outcome (say, the odds of picking n white balls from a bag containing m balls in total) and a variety of other binary parameters:

Dear R people, This is probably a FAQ, and a silly question besides, but... a lot of R packages have datasets contained within them. cf. spatial has "topo". Once I have loded the appropriate library, how do I access the datasets contained within? I know this must be very simple, but I have spent some (fruitless, frustrated time) trying to do this, so am asking for help.

Ok. I did some quality tests. First off; never do quality tests with ints. I had serious problems interpreting my results until it dawned on me that the signal differences were just 0 or 1. So, after a lot of scratching my head, these are done comparing the result from the _float versions (which is how both resamplers work internally anyway). What I did was this: Load speex_wb.wav as one

>> Yes, I plan to use it in a VoIP environment if I can get latency reduced to >> an acceptable level :) >> The latency depends directly on the overlap parameter, which also controls >> the quality. Higher quality => higher latency. You could set the overlap to >> 0, but that would give you some nasty artifacts. >> You can also resample with smaller block

Can you explain why you need them as 'integer', A floating point representation can hold a value upto ~4.5e15 as an "integer" keeping the precision that you might need. Jim Holtman Data Munger Guru What is the problem that you are trying to solve? Tell me what you want to do, not how you want to do it. On Tue, Apr 26, 2016 at 1:11 PM, Andr? Luis Neves <andrluis at

Hallo, I have this dataset with repeated measures. There are two within-subject factors, "formant" (2 levels: 1 and 2) and "f2 Ref" (25 levels: 670, 729, 788, 846, 905, 1080, 1100, 1120, 1140, 1170, 1480, 1470, 1450, 1440, 1430, 1890, 1840, 1790, 1740, 1690, 2290, 2210, 2120, 2040, 1950), and one between-subject factor, lang (2 levels:1 and 2). The response variable

Dear R-help, I am having trouble with your scatterplot3d program. For help with this problem I was directed to your address by Martin Maechler at " r-core-bounces at r-project.org." I'm also sending a CC to " r-core-owner at r-project.org" as I'm not yet certain of the proper address to use for this. I have R version 2.8.1 and have downloaded 'scatterplot3d.'

R Help - I'm trying to run a correlation matrix with a covariate of "age" and will at some point will also want to covary other variables concurrently. I'm using the "psych" package and have tried other methods such as writing a loop to extract semi-partial correlations, but it does not seem to be working. How can I accomplish this? library(psych) > set.cor(y =

Ok. I`m trying to run a Poisson glmm with an observation-level random intercept. But I`m getting the following error for the 'Baci' variable: 'Error: (maxstephalfit) PIRLS step-halvings failed to reduce deviance in pwrssUpdate'. I guess this message is because the baci variable is not a an integer, and cannot be transformed into an integer as R has a threshold of 2x10^9 even in