similar to: Returning "overall ms" for CVlm

Hello, I did a cross-validation using cvlm from DAAG package but wasn't sure how to assess the result. Does this result means my model is a good model? I understand that the overall ms is the mean of sum of squares. But is 0.0987 a good number? The response (i.e. gailRel5yr) has min,1st Quantile, median, mean and 3rd Quantile, and max as follows: (0.462, 0.628, 0.806, 0.896, 1.000, 2.400) ?

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

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

Hiya, I am trying to look at the similarities between a number of sequences, for example i am trying to see how similar "ababbbassdaa" is to "addffggssbbsbbs" I was wondering is the some way for me to see how similar they are in terms of, for example, number of a's, number of b's, how often a and ab are consecutive, how often abab is together etc. Any advice

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

Hi everyone, I was trying to mode the following exercise using R. The question: Set up a one pool model using numericintegration. The model will run from time 1 to time 30 using a time step of 1.The pool (A) will be fed by flux "inA" at a rate of 5 units per hour anddrained by flux "outA" at a rate of 20% per hour. At time 0, A has 5units. At time 30, what is the pool size of

Hi R help, I wanted to simulate two pool model (A&B) using deSolve package for time 0 to 12 by 1.? Initial values of the state variables are A=5, B=3. The fluxes are as follows1) Flux into A= 5 units per unit time?2) Flux from A to B= 0.33) Flux out of A=0.1?4) Flux from B to A=0.35) Flux out of B=0.3 Here is the R code I compiled to estimate the size of A and B and graph the output

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

Add this: --delete-excluded See the manual. Rainer -----Urspr?ngliche Nachricht----- Von: rsync-bounces+rsync=diplan.de@lists.samba.org [mailto:rsync-bounces+rsync=diplan.de@lists.samba.org]Im Auftrag von Mandrake Pe Gesendet: Mittwoch, 29. Oktober 2003 13:33 An: rsync@lists.samba.org Betreff: trying to figure out how the --delete option works. I use this comand to back up my entirely drive

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:

Hi there, I would like to make a 'bumps chart' like the ones described e.g. here: http://junkcharts.typepad.com/junk_charts/bumps_chart/ Purpose: I'd like to plot the proportion of people in select countries living for less then one USD pr day in 1994 and 2004 respectively. I have already constructed a barplot - but I think a bumps chart would be better # The barplot and data

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

I have two Linux computers suse and RH enterprises. The first one (suse) works as PDC and everything is OK with PDC, all Shares work correctly. Now I have installed second samba Server on RH Enterprise and try set it as domain memember, so I set the security on domain. As I read, if this second samba stay a member of domain every share I try to rich would be checked with password on domain I mean

I use this comand to back up my entirely drive to a second drive. Because the second drive is smaller I have excluded my movies directory. I also excluded the /mnt/ diectory, becous I dont want to back up the backup. This part works just fine. It's the --delete option I cant get to work. I want rsync to delete files on the backup(destination) drive, that I have deleted or moved on the

Dear all, I have a 3x2 plot and in addition to the title of the individual plots I would like to have an overall title for each row. I managed to get an overall title for the whole plot matrix with mtext: par(mfrow=(c(3,2)), mar=c(6.4,4.5,4.2, 1.8), oma=c(0,0,3,0)) for (i in 1:6) barplot(sample(1:10,5), main=paste("Plot #",i)) mtext("Overall Title Row 1", outer=TRUE) but I

Dear useRs, I am currently fitting an advanced failure time model using G?ran Brostr?m's excellent "eha" library with the "aftreg" command. My question: How do I interpret the "Overall p-value", that is reported at the very bottom of the output? I already figured out it must be a chi-square test, but I am wondering what a p-value < 0.01 means: Does it mean

Hello fellow R users, I am having a problem finding the estimates for some overall treatment effects for my mixed models using 'lme' (package nlme). I hope someone can help. Firstly then, the model: The data: Plant biomass (log transformed) Fixed Factors: Treatment(x3 Dry, Wet, Control) Year(x8 2002-2009) Random Factors: 5 plots per treatment, 5 quadrats per plot (N=594 (3*5*5*8)

In a CentOS 5 dom0, how does one find out the overall system load? In the early days of Xen, I believe an uptime on dom0 provides an overall load for the server. This doesn't seem to be the case anymore. -------------- next part -------------- An HTML attachment was scrubbed... URL: <http://lists.centos.org/pipermail/centos/attachments/20070525/e4c30b4b/attachment.html>

Hi all, I know the overall display precision can be changed in R... but what about overall calculation precision? Thank you! [[alternative HTML version deleted]]

Hello, I have two classifications. How can I compare the overall-accuracy of these classifications to each other? Is there a possibility within R to test if the achieved overall-accuracy for the classifications are differing significantly? Additionaly, are the McNemar-Test and Broker-Test implented in a package of R? Thank in advance for your help, Markus