similar to: Change class factor to numeric

Hello I am able to extract partial regression coefficients from a fitted model object "model", i.e. model <- lm(var.sel.gkm, weights = count.gkm, data = DATA) summary(model) write.table(model$coef, file = "C:/coef_CO_gkm.txt", row.names = TRUE, col.names = TRUE) I was wondering if anyone could advise me how to extract other object items such as std. error, t-values

Dear List members I am using the regsubsets function to select a few predictor variables using Mallow's Cp: > sel.proc.regsub.full <- regsubsets(CO2 ~ v + log(v) + v.max + sd.v + tad + no.stops.km + av.stop.T + a + sd.a + a.max + d + sd.d + d.max + RPA + P + perc.stop.T + perc.a.T + perc.d.T + RPS + RPSS + sd.P.acc + P.dec + da.acc.1 + RMSACC + RDI + RPSI + P.acc + cov.v + cov.a +

I am hoping to get some advise on the following: I am looking for an automatic variable selection procedure to reduce the number of potential predictor variables (~ 50) in a multiple regression model. I would be interested to use the forward stepwise regression using the partial F test. I have looked into possible R-functions but could not find this particular approach. There is a function

Thank you all for your (fast) comments. Unfortunately I could not make the advise work: > mass [1] 800 800 800 800 800 800 800 800 800 800 800 800 800 800 800 800 800 800 910 910 910 910 910 910 910 [26] 910 910 910 910 910 910 910 910 910 910 910 910 910 910 1,020 1,020 1,020 1,020 1,020 1,020 1,020 1,020 1,020

Hello, The estimate of glm dispersion can be based on the deviance or on the Pearson statistic. I have compared output from R glm() to another statastical package and it appears that R uses the Pearson statistic. I was wondering if it is possible to make use R the deviance instead by modifying the glm(...) function? Thanks for your attention. Kind regards, Robin Smit This e-mail and its

Hello I am trying to use all subsets regression on a test dataset consisting of 11 trails and 46 potential predictor variables. I would like to use Mallow's Cp as a selection criterion. The leaps function would provide the required output but does not work with this many variables (see below). The alternative function regsubsets should be used, but I am not able to define the function in

Dear list members, I apologize for putting this (probably) very basic question on the mailing list. I have scanned through the R website (using search) but did not found an answer. (code included below) A factor matrix is simply extracted (which can then subsequently be exported using write.table) by FACT$loadings[1:6,]. I would also like to specifically extract and export

Hello, I have got a simple "cosmetic" question. I have created a bar plot with confidence intervals using: barplot(mean, ylim = c(0,0.2), las = 3, space = 0) arrows(1:17 ,X95p_low, 1:17, X95p_high, length = 0.07, angle = 90, code = 3, lty = 1) Unfortunately, the confidence bars align with the right side of the bars and I would like to shift these to the middle of the bars. Is

Dear mailing list member I have got the following question: The function "cor.test(x1, x2)$p.value" computes the statistical significance of the correlation between vectors x1 and x2. However, I have got a data matrix with many variables. Instead of typing in all possible combinations, I was wondering if there is perhaps an R-function that produces a matrix showing the

Hi to all, very strange! I tried to debug why imap process takes 100% CPu and it is very slow when I move an email from an imap folder to another imap folder and then I typed: strace -tt -o /log-imap.txt -p <process number> well, during the trace the imap process doesn't crash and it is very fast. If I stop the tracing, I have again the problem and sometimes the crash of the

Hi, I'm just a beginner, who has just encountered a problem! 1. - I wanted to load a csv file with names in the rows (1st column) and and numbers in the 2nd til 10th column. The file contains names in the headers. - I used; a <- as.matrix(read.table("filename", sep=',", row.names=1, header=TRUE) Question; 1 - I would like to select the first four columns 2 -

Does anyone have nice quality controlls for affy arrays,.... Can't find any tools as are being used for 2 dye arrays..... cheers, marinus This e-mail and its contents are subject to the DISCLAIMER at http://www.tno.nl/disclaimer/email.html [[alternative HTML version deleted]]

Thank you for your reply. I am aware of the good reasons not to use the deviance estimate in binomial, Poisson, and gamma families. However, for the inverse Gaussian, the choice seems to me less clear cut. So I just wanted to compare two different options. I have used the dispersion parameter to compute the standardized deviance residuals: summary(model.gamma)$deviance.resid

On Wed, Feb 11, 2015 at 06:41:20PM +0800, Chen Hanxiao wrote: > + /* Output of `btrfs scrub -R status' is like: > + * > + * scrub status for 346121d1-1847-40f8-9b7b-2bf3d539c68f > + * scrub started at Mon Feb 2 17:39:38 2015, running for 93 seconds > + * data_extents_scrubbed: 136670 > + * tree_extents_scrubbed: 30023 > + *

While setting my first couple IP phones, I set their voicemail buttons to an extension that runs VoicemailMain. exten => 8500,1,Wait(1) ; voicemail exten => 8500,2,VoicemailMain ; exten => 8500,3,Hangup ; I would like to be able to pass the mailbox number allowing each phone to go in directly but I'd rather tno have

Peoplefone AG jest firm? oferuj?c? us?ugi telekomunikacyjne oparte na technologii Voice over IP (VoIP) w wyj?tkowo niskich stawkach. Peoplefone jest certyfikowanym partnerem firm Siemens<http://www.siemens.ch/index.jsp?sdc_p=c175fi1012637lmno1012637psuz1&sdc_sid=1113876080&>i AVM/FRITZ!Box <http://www.fritz-shop.ch/>. W zwi?zku z naszym dynamicznym rozwojem, do nowo tworzonej

oops, Forgot to cc to the list. Regards, Mike -----Original Message----- From: dr mike [mailto:dr.mike at ntlworld.com] Sent: 24 February 2005 19:21 To: 'Spencer Graves' Subject: RE: [R] Bayesian stepwise (was: Forward Stepwise regression based onpartial F test) Spencer, Obviously the problem is one of supersaturation. In view of that, are you aware of the following? A Two-Stage

I have defined a function to compute the value of a beta distribution of the second kind (the existing beta distribution of th stats package is the beta distribution of the first kind). It works perfectly for a single value, but I want to apply it to a vector of 22 000 values. I can use a loop for the calculation of each value but it runs very very slowly. So, what can I change ? Hers's the

I've enconuntered this problem with the last cran version of geoR: > library(geoR) > day <- rep(1:2, each=5) > coords <- matrix(rep(runif(10),2), 10, 2) > data <- rnorm(10) > data[1] <- NA > as.geodata(cbind(coords, data, day), realisations=4) as.geodata: 1 points removed due to NA in the data Errore in as.geodata(cbind(coords, data, day), realisations = 4) :

Thanks Bart. Unfortunately "The VIB module does not comply with the security recommendations imposed by VMWare. You lose VMWare support by installing this package." This pretty much rules out this option. :( We need vCenter support, as well. nomad On Thu, May 14, 2020 at 10:06 AM Bart J. Smit <bart at smits.co.uk> wrote: > This shuts down the hosts directly (by-passing