similar to: using cvlm to do cross-validation

Hi All, when doing CVlm cross validation, I would just like to view the overall MS but i am having trouble getting just this value. I'm kina new to R. Any help really appreciates. Cheers, Ben -- View this message in context: http://r.789695.n4.nabble.com/Returning-overall-ms-for-CVlm-tp4664690.html Sent from the R help mailing list archive at Nabble.com.

I would like to run a loop in R. I have never done this before, so I would be very grateful for your help ! 1. I have a sample set: 25 objects. I would like to draw 1 object from it and use it as a test set for my future external validation. The remaining 24 objects I would like to use as a training set (to select a model). I would like to repeat this process until all 25 objects are used as a

Can anyone help explain to me why the two codes below have different result? I thought I can use log(time)~. to replace log(time)~dist+climb+timef.I am using CVlm from DAAG package. I think nihills is preloaded with the package. Thanks in advance. > CVlm(df=nihills, form.lm=formula(log(time)~.),plotit="Observed",m=2)Analysis of Variance Table Response: log(time) Df Sum Sq

Hello, I used the cv.lm function to validate a linear regression model fit<-lm(y ~ x1+x2+x3+x4+0, data=mydata) without intercept I tried to validate the model by performing a leave one out cross validation procedure usinfg the cvlm function: CVlm(df=mydata, fit, m=196) But how can I get the adjusted R? from the output of this function. Or is there any other function to

Dear all, I've tried a multiple regression, and now I want to try a cross-validation. I obtain this error (it must be sth related to df) that I don't understand, any help would be appreciated. cv.lm(df= dat, lm2.52f, m=3) Error en `[.data.frame`(df, , ynam) : undefined columns selected lm2.52f is my lm object, dat is a dataframe where the variables involved in .lm are I tried CVlm

Guys, I have setup my hard disc with 3 partitions: 1- 256MB on /boot; 2- 2GB on / for my dom0 (Debian 5.0) (eth0 default bridge for guests LAN); 3- 498GB exported with vblade-persist to my network (eth1 for the AoE protocol). On dom0 hypervisor01: vblade-persist setup 0 0 eth1 /dev/sda3 vblade-persist start all How to create a CVLM VG with /dev/etherd/e0.0 on each of my dom0s? Including the

Dear community, I have fitted a model using comands above, (rlm, lmrob or lmRob). I don't have new data to validate de models obtained. I was wondering if exists something similar to CVlm in robust regression. In case there isn't, any suggestion for validation would be appreciated. Thanks, user at host.com -- View this message in context:

Dear R users, we are trying to do some parallel computing using library(snow). In particular we have a cluster with 3 nodes >cl <- makeCluster(3, type = "MPI") 3 slaves are spawned successfully. 0 failed. and we want to compute the function op_mat (see below) first with the master and then with the cluster using system.time for checking the computational performance.

Hi, I am trying to find m breakpoints of a linear regression model. I used the segmented package. It works fine for small number of predicators and breakpoints.(3 r.v. 3 points). However, my model has 14 variables it even would not work even for just one breakpoints!. The error message is always estimated breakpoints are out of range. Since my problem is time related problem. So I

Having spent the last few weeks trying to decipher R, I feel I may finally be getting somewhere, but i'M still in need of some advice and all my tutors seem to be on holiday! Basically a bit of background, I have data collected on a population of Lizards which includes age,sex, and body condition. I collected data myself this year and I have data previously collected from 1999, 2002 and

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 =

Hi, I have a line such as: myline <- " 0.100000 1.5000 0.6000 538 0.369404" and I would like to put the numbers into a vector. Some combination of tabs and spaces occur between the numbers. I tried: try1 <- strsplit(myline,"[[:blank:]]+") > try1 [[1]] [1] "" "0.100000" "1.5000"

Hi, I am trying to make a multi-class classification tree by using rpart. I used MASS package'd data: fgl to test and it works well. However, when I used my small-sampled data as below, the program seems to take forever. I am not sure if it is due to slowness or there is something wrong with my codes or data manipulation. Please be advised ! The data is described as the output from str()

All, When using grid.text it seems my supplied (x,y) coordinates are being plotted only in npc (normalized parent coordinates) where (.5,.5) is the center of the graph. How do I allow (x,y) to be coordinates corresponding to the (x,y) values in the graph? The examples in ?grid.text seem to do this. (I think the problem lies in how x and y are defined with respect to unit ). Any help much

All, How does one replace plot symbols with say subject IDs when using xyplot? Or superimpose them next to plot symbols? I searched the archives under various key words but haven't had much. Any suggestions or links much appreciated. Sample code below. David junk.frm = data.frame(ID = rep(1:16, each = 2), x, y, z = rep(c("D", "P"), 16)) y = c( 0.4, 0.6, -0.1,

All, I'm having problems w/ a simple attempt to subset an xyplot. The first plot below is a plot of y versus x for certain values of a third categorical variable z. Now I'd like to further restrict this to certain values of variable y. Neither of the two attempts below work. Any suggestions much appreciated. (note: I don't want to merely use ylim since I have a loess plot and I

Hi there, as a student sociology, I'm starting to learn about SEM. The course I follow is based on LISREL, but I want to use the SEM-package on R parallel to it. Using LISREL, I found it to be very usable to be able to see the total direct and total indirect effects (standardized and unstandardized) in the output. Can I create these effects using R? I know how to calculate them

Hi all, Im a little bit confused concerning the effect() command, effects package. I have done several glm models with family=quasipoisson: model <-glm(Y~X+Q+Z,family=quasipoisson) and then used results.effects <-effect("X",model,se=TRUE) to get the "adjusted means". I am aware about the debate concerning adjusted means, but you guys just have to trust me - it

Hi, LLVM 2.1-pre1 test results: Linux (SUSE) on x86 (P4) Release mode, but with assertions enabled LLVM srcdir == objdir # of expected passes 2250 # of expected failures 5 I ran the llvm-test suite on my desktop while I was also working on that PC, so don't put too much trust in the timing info. Especially during the "spiff" test the machine was swapping

Hello R users, For the last few days I am struggling with the following task: my data.frame: A1 A2 A3 B1 B2 B3 58.81 53.292 54.501 13.013 17.39 19.407 56.02 56.251 54.033 20.099 13.15 10.411 55.376 53.099 57.625 13.396 21.031 13.22 58.584 53.194 54.218 13.038 16.854 19.289 55.7 55.921 53.847 19.942 13.153 9.828 55.093 52.934