similar to: questions about performing Robust multiple regression using bootstrap

Dear Faiz, Bootstrapping R^2 using Boot() is straightforward: Simply write a function that returns R^2, possibly in a vector with the regression coefficients, and use it as the f argument to Boot(). That will get you, e.g., bootstrapped confidence intervals for R^2. (Why you want that is another question.) See the example in ?Boot that shows how to bootstrap the estimated error variance (without

I am trying to perform goodness of fit test using R. I am using this website http://wiener.math.csi.cuny.edu/Statistics/R/simpleR/stat013.html for help. However, I am unable to carry out the test successfully. My code follows. It is taken from the website just mentioned. freq=c(22,21,22,27,22,36) # frequencies obtained after rolling the dice 150 times. prob=c(1,1,1,1,1,1)/6 # specify expected

Hi all, I want to make a contingency table in R. I want to tabulate two variables, one as the independent and second as the dependent variable. The IV has two categories, namely, birth complications, and no birth complications. The frequency of birth complication category is fifty, and the frequency of no birth complication category is 34. The categories and frequencies of DV follows.

Dear list members, I want to save residuals above or less than a certain value to an R object. I have performed a multiple linear regression, and now I want to find out which cases have a residual of above + 2.5 and ? 2.5. Below I provide the R commands I have used. Reg<-lm(a~b+c+d+e+f) # perform multiple regression with a as the dependent variable. Residuals<-residuals(reg) # store

Hi all, I would like to read the results of the commands type in the terminal window in Microsoft Word. As a blind user my options are somewhat limited and are time consuming if I want to see the results of the commands that I have type earlier. for example if my first two commands were x<-c(1,2,3,4,5) mean(x) and I have typed ten more commands after the first two commands it is not easy

Residuals are stored as a numeric vector. The R software comes with a document "Introduction to R" that discusses basic math functions and logical operators that can create logical vectors: abs( stdresiduals ) > 2.5 It also discusses indexing using logical vectors: stdresiduals[ abs( stdresiduals ) > 2.5 ] Note that in most cases it is worth going the extra step of making your

Dear list, I have the following how-to-do it in R, questions. Suppose I have ten independent variables, and one dependent variable. I want to find the Pearson correlation of all the IVs with the DV, but not the correlation between the IVs. What I know so far, about R, that I have to type the cor () function ten times, each time requesting for a correlation between one IV and the DV. I was

Hi For first question, maybe I am completely wrong but cor(swiss[,-1], swiss[,1]) should give you what you want in one step. Second question Without an example it is hard to say but maybe aggregate is the way forward. > aggregate(iris[,1:4], list(iris$Species), function (x) c(mean=mean(x), sd=sd(x))) Group.1 Sepal.Length.mean Sepal.Length.sd Sepal.Width.mean Sepal.Width.sd 1

Hi all, I am trying to use type message with sink, like this: sink("all.Rout", type="message") 1+3 sink() readLines(con = "all.Rout") So to get the following output: > 1+3 [1] 4 Obviously this doesn't work. I tried some variations (based on the explanations in the help) but am missing something on how to make it work. Any suggestions? (p.s: I need

After reading more, I understand I didn't formulate my last question correctly, so please allow me to rephrase: What I am looking for is a way to save the R console session output. That is, a command that would combine the results of using: ?sink # And ?savehistory My motivation for this is that doing it will allow someone who is a blind user of R to be able to easily export his results to

Hi all, I've got some questions to post in regard to running asterisk in a production-grade environment, specifically targeting high-density IVR applications. No VoIP involved, just straight PSTN -> * and perhaps the occasional outdials or agent-based predictive dialing. 1) Which user would you run * under? 2) What other security-related issues do you have to resolve? 3) How do you handle

Hello, I''m experimenting with a xhtml templating language and have some questions about ActionPack, particularly to do with compilable templates. I''m mostly looking for further information about compilable templates. The templating language currently works, though it isn''t a compilable template handler. I''ve yet to benchmark it, but it seems to be running

hi there, i've been able to successfully run asterisk with the Voicetronix OpenLine4 card, it can accept calls and function normally. The only problem I'm experiencing so far is getting the card to outdial to a third party. What I'm trying to achieve is basically call bridging, where the caller dials in to asterisk, some IVR plays and then attempts to perform a "transfer"

I am fitting quantile regression models using data collected from a sample of 124 patients. When modeling cross-sectional associations, I have noticed that nonparametric bootstrap estimates of the variances of parameter estimates are much greater in magnitude than the empirical Huber estimates derived using summary.rq's "nid" option. The outcome variable is severely skewed, and I am

Hi, I''m looking for past efforts to integrate Rails and Needle. I am thinking of doing it myself but would be interested in seeing what others have done before. If there are examples of other IoC frameworks being used with Ruby I''d be interested in hearing about them too. I find IoC to be a helpful concept in general and, having glanced at a few existing Rails app.s, I

I have an advanced question about bootstrapping. There are two datasets. In each bootstrap iteration, I would like to sample One observation per cluster from the first dataset. N observations with replacement from the second dataset. Right now I am using dplyr::sample_n() for first dataset, with this sampling embedded in the program that boot() from the boot package is running to sample the

Apologies if this messages appears more than once, somehow sending it via gmane's nntp-server doesn't work as expected. Been trying to send this message for 2 days now, without success so far... --- Hello, I'm exploring the capablities of Asterisk and must say I'm really impressed! However, I don't need most of the options, but can't figure out the things that appear

I am trying to get robust standard errors in a logistic regression. Is there any way to do it, either in car or in MASS? Thanks for the help, Celso [[alternative HTML version deleted]]

Hello, Hacker I install VoiceTronix OpenSwitch 12 port PCI Telephone Card, and setting vpb.conf, extensions.conf following My problem is: When i dial to fxo(channel 9-12), it is ok, but when i continue press exten 102, the channel crach with error messages following exception caught: VPBAPI_DIAL_INVALID_LEVEL, file: vpbdial.cpp line:872 Do i ignore some setting for VoiceTronix OpenSwitch12

Hello,   Anyone familiar with robust ANOVA in R? Please help.   I am conducting a between-between-within (2 X 2 X 2) ANOVA, with the focus on the significance of the three-way interaction. Because of very uneven sample sizes, and the violation of covariance homogeneity, I need to use some sort of robust ANOVA method, and I prefer bootstrap. I am aware of Rand Wilcox's robust package that can