similar to: Help - lm, glm, aov results inconsistent with other stati stical package

Hello, I am new to using R and I am having problems get boot() to work properly. Here is what I am trying to do: I have statistic called "cs". cs takes a data matrix (154 x 5) and calculates 12 different scores for me. cs outputs the data as a vector (12 x 1). cs doesn't really use weights, per se, however I have included this as one of the 2 arguments cs can take. I try

Hello, I 'm sure there must a be a simple explanation for what I'm doing wrong but I am stumped. I am a novice R user and this has shaken my confidence in what I'm doing! I am trying to run a simple ANCOVA using the model y~A*x, where y and x are continuous and A has two levels. Everything seems fine until I compare the output with what I get from another statistical package

Hello all, I'm trying to establish some confidence intervals on predictions I am making using >predict(nls(...)) and predict.nls (unfortunately) does not utilize the se.fit option. A little more background is that I am trying to match the output with older SAS routines to maintain consistency. Because predict.nls does not provide se's for individual predictions, I have been using a

Hi, My R version is old, and I would like to update it. How can I update R and keep all the libraries? Thanks Best, Yunfei Li -------------------------------------------------------------------------------------- Research Assistant Department of Statistics & School of Molecular Biosciences Biotechnology Life Sciences Building 427 Washington State University Pullman, WA 99164-7520 Phone:

All, I am working with NIR spectral data and it was great to find that the example in ?plsr also used spectral data. Unfortunately, I am having difficulty figuring out how the "yarn" dataset is structured to allow for the plsr model to read: library(pls) data(yard) yarn.oscorespls <- mvr(density ~ NIR, 6, data = yarn, validation = "CV", method = "oscorespls")

I am attempting to create a formula using as.formula for a PLS analysis. I have used the code below successfully, but in a previous R version and with many fewer predictors. Any help getting all of these predictors into one formula would be greatly appreciated. TC.fmla <- as.formula(paste("TC ~ ", paste(vars, collapse= "+"))) As I mentioned, this code worked fine in a

Hi, I am studying using SVM functions of e1071 package to do prediction, and I found during the training data are "factor" type, then svm.predict() can predict data directly by categories; but if response variables are "numerical", the predicted value from svm will be continuous quantitative numbers, then how can I connect these quantitative numbers to categories? (for

I am attempting to run principal components analysis on a dataset of spectral reflectance (6 decimal places). I imported the data using read.table and there are both column and row headers. When I run princomp I receive the following error: Error in cov.wt(z) : 'x' must contain finite values only Where am I going wrong? Ross

All, I have been using a data set to build pls models for three different soil properties. Two of the three models run fine; however I receive the following error for the final model. > libs.IC.cal <- mvr(libs.IC.fmla, data = libsdata.cond.cal, ncomp=20,validation = "LOO", method = "oscorespls") Error in colMeans(x, n, prod(dn), na.rm) : 'x' must

I am attempting to write a routine where I can run PAM (partition around mediods) on a dataset containing multiple soil cores and PCA spectral data from several depths per core. I want to run PAM on each individual core, so I need to group the data by core to run the analysis. Below is an example of my data structure: Lab.id PC1 PC2 PC3 MAT057.2.5 2.438454966 -1.011182986 -3.040881377

Hi, I tried to add a progress bar to my script, but it seems the "txtProgressBar" function will wipe out what I want to print before it(for example - the "Hello" in following script), and the parameters "title" and "label" do not work either. How can I put a title for the progress bar? or is there any other better function I can use? testit <-

Hello all, I'm sure I'm missing something simple here, but I can't figure out how to modify the glm.fit() function and then get R to use it (sort of). I'm doing something along the lines of: glm.fit<-edit(glm.fit) # add something trivial to the top of the glm.fit function like: print("Hello world!") #now have a modified glm.fit in position 1/.GlobalEnv

Hi all, I have a question regarding how to properly analyze a data set and then how to perform the analysis in R. First, I have data that I would like to analyze using a mixed GLM (I think this is the most appropriate method, but I am unsure). In a mixed model (y = X*beta+Z*gamma+epsilon), I would like to structure the variance matrices of gamma, G, and the error, R, to take advantage of all my

I am new to Linux (not to R) and recently installed Mandriva Spring 2007 on my partitioned hard drive. My next objective is to install R in the Linux environment, unfortunately Mandriva is not one of the Linux distributions available for download... Could someone please let me know which distribution I should use? Thanks. Jonathan --------------------------------- Boardwalk for $500?

snip from make.log: Compiling libsmb/libsmbclient.c with -KPIC "libsmb/libsmbclient.c", line 3249: warning: argument #3 is incompatible with prototype: prototype: pointer to char : "libsmb/libsmbclient.c", line 1167 argument : pointer to const char "libsmb/libsmbclient.c", line 3843: warning: argument #5 is incompatible with prototype:

Hi, i'm having trouble using Statis in ADE4 package. I want to study a matrix with 18x414...could any body help me?!?!' Urgent Thanks Miguel Ribeiro Portugal [[alternative HTML version deleted]]

Noticed that the sites has not been responding correctly for a few days, and was wondering if anyone knew what was up. traceroute and ping get to the ip, but website times out after a very long time. +----------------------------------------------------------+ Michael D. Setzer II - Computer Science Instructor Guam Community College Computer Center mailto:mikes

Users of R may find the "Message from the President" in the latest issue of Liaison, the newsletter of the Statistical Society of Canada, to be of interest. (http://www.ssc.ca/documents/pubs/Liaison17.2.pdf). In this message Jim Ramsey, President of the SSC, describes his purpose in inviting Robert Gentleman to give the President's Invited Address at their annual general meeting.

Also Consider Bonferroni Hochberg Holm type procedures or . Dunn OJ. Multiple contrasts using rank sum tests. Technometrics 1964;6:241#/52. [[alternative HTML version deleted]]

All, I am a new user to Linux but I am familiar with R. I have previously used and installed R on a Windows platform without problems. I recently set up a dual boot system (XP_64, Mandriva) to run R on a Linux platform in order to more efficiently handle large datasets. I have not done compiling before, but read the R instructions and followed to my best ability. I downloaded the most recent