similar to: Stat textbook recommendation

I am writing a resampling program for multiple regression using lm(). I resample the data 10,000 times, each time extracting the regression coefficients. At present I extract the individual regression coefficients using brg = lm(Newdv~Teach + Exam + Knowledge + Grade + Enroll) bcoef[i,] = brg$coef This works fine. But now I want to extract the t tests on these coefficients. I cannot

I have been having problems with using scan(). I searched the archives and found someone with the same problem several years ago, but did not find a solution. I want to prompt the user for input, scan in that input, and then go on to use it as a variable. The simplest version of my program is cat("enter m","\n") m <- scan(n = 1, quiet = T) m If I enter these lines one

Sampling with and without replacement I seem unable to use "replace = F" when I want to sample without replacement. I would think that it comes down to "F is not a legitimate abbreviation for FALSE." except that Dalgaard (p. 118) uses F for FALSE and it works "pairwise.t.test(folate, ventilation, pool.sd = F)" I am having trouble when I try to sample a vector

Hi, I am getting some strange results using round - it seems that it depends if the number before the decimal point is odd or even .... For example: > round(1.5)[1] 2> round(2.5)[1] 2 While i would expect that round(2.5) be 3 and not 2. Do you have any explanation for that? I really appreciate your input, Monica

Dear R user I have 2 problems with lmer. The statistical consultance service of my university has recomended to me to expose those problems here. Sorry for this quite long message. Your help will be greatly appreciated... Gilles San Martin 1) anova() I fit a first model : model1 <- lmer(eclw~1 + density + landsc + temp + landsc:temp + (1|region) + (1|region:pop) + (1|region:pop:family),

Hello List, I'm trying to prepare some lecture notes on non parametric methods, and I can't manually reproduce the results of the wilcox.test function for ordinal data. The data I'm using are from David Howell's website, available here http://www.uvm.edu/~dhowell/StatPages/More_Stuff/OrdinalChisq/OrdinalChiSq.html If I run the wilcox.test function on the data I get a p-value of

Hello list, I wish to perform a randomization test on the F-statistics of a 2 way ANOVA but have not been able to find out how to do so - is there a package / function that can perform this that I am unaware of? FactorA has 6 levels (0,1,2,3,4,5) whereas FactorB has 3 (1,2,3). A sample: Resp. FactorA FactorB 2 0 2 3 1 2 1 2 2 0 3 2 0 4 2 0 5 2 4 0 1 6 1 1 1 2 1 0 3 1 1 4 1 0 5 1 2 0 2 3 1 2 1

Have Asterisk connected to a H323 compatible legacy PBX using QSIG protocol and IP trunks. I can call to Asterisk, and from Asterisk using X-Lite softphone but whenever either end picks up, the calls disconnects. No gatekeeper is installed. I have attached a copy of my h323 logfile for debugging. What do you suggest what change needs to take place to keep calls connected? 11:33:19:864

#I am trying to understand how R fits models for contrasts in a #simple one-way anova. This is an example, I am not stupid enough to want #to simultaneously apply all of these contrasts to real data. With a few #exceptions, the tests that I would compute by hand (or by other software) #will give the same t or F statistics. It is the contrast estimates that R produces #that I can't seem to

I've installed R 1.8.1 on OSX 10.3 (Panther) using the RAqua.pkg and all indications were that the installation was successful. However, after double clicking the R icon in the Applications folder, nothing appears to happen. The following message appears on the console: Mac OS X Version 10.3 (Build 7B85) 2003-12-11 08:45:31 -0500 prefs written dyld:

I've installed R 1.8.1 on OSX 10.3 (Panther) using the RAqua.pkg and all indications were that the installation was successful. However, after double clicking the R icon in the Applications folder, nothing appears to happen. The following message appears on the console: Mac OS X Version 10.3 (Build 7B85) 2003-12-11 08:45:31 -0500 prefs written dyld:

Hi All, I am struggling for 2 weeks solving authentication problem in dovecot. logs from /etc/mail/maillog Oct 15 18:00:35 localhost dovecot: auth: Debug: Loading modules from directory: /usr/lib64/dovecot/auth Oct 15 18:00:35 localhost dovecot: auth: Debug: Module loaded: /usr/lib64/dovecot/auth/libauthdb_ldap.so Oct 15 18:00:35 localhost dovecot: auth: Debug: Module loaded:

I've downloaded and installed R 1.9.0 under OS X 10.3.3. I then removed my previous R 1.8.0 installation by following the instructions in the readme file, e.g., 1. drag the "StartR" icon into the Trash; 2) go inside your Home directory and search for the Library folder. Inside the Library folder in your Home directory you'll find a RAqua directory: drag the RAqua

Hello everybody, I use the sweetpotato database included in R package: data(sweetpotato) This dataset contains two variables: yield(continous variable) and virus(factor variable). Due to Levene test is significant I cannot assume homogeneity of variances and I apply Welch test in R instead of one-way ANOVA followed by Tukey posthoc. Nevertheless, the problems come from when I apply posthoc

Dear R News, I am inquiring whether the articles that are published in the 'R News' have undergone a peer-review/refereeing process, for the purpose of reporting the academics publications here at the University of Sydney to the Australian Government. Thank you for your assistance. Regards Glenn Glenn Howell| Faculty Research Assistant FACULTY OF ECONOMICS AND BUSINESS

Hello, I'm trying to install a package that requires a Fortran compiler (Hmisc) using R CMD INSTALL. I downloaded the package source onto my Desktop, unzipped it, and then typed: R CMD INSTALL /Users/brianbeckage/Desktop/Hmisc * Installing *source* package 'Hmisc' ... ** libs g77 -fno-common -g -O2 -c cidxcn.f -o cidxcn.o g77 -fno-common -g -O2 -c cidxcp.f -o cidxcp.o g77

1) When adding additional data sets to a plot using "plot" followed by "lines", is there a way to automate the scaling of the axes to allow for all data sets to fit within the plot area? 2) I attempted to solve this by setting xlim=c(min(c(data1,data2,data3)),max(c(data1,data2,data3))) however, there are some NAs and Infs in these data sets, and min(data1) and max(data1) both

Is it possible to give a temporary object a name that varies with each run of a foreloop? For example, I want to fill a matrix every time I run a loop, and I want a new matrix with each run, with an appropriate new name. i.e.: for(i in 1:5){... matrix.i<-some values ...} so that in the end I would have: matrix.1 matrix.2 matrix.3 matrix.4 matrix.5 Thanks, Ben Osborne -- Botany Department

Hi Timo & everyone, I'm trying out a 2-node director setup, but I keep getting the following error: Oct 3 16:11:29 imapdir1 dovecot: master: Dovecot v2.0.15 starting up (core dumps disabled) Oct 3 16:11:34 imapdir1 dovecot: director: Error: connect(132.198.100.150:9090) failed: Invalid argument Oct 3 16:11:41 imapdir1 last message repeated 3 times Both nodes report this error. The

Hello all, I know I need to read the manual but I don't even know where I need to look in the manual. I'm haveing a little bit of a problem getting XP to print to my linux printer. I've got the smb.conf file set up to allow xp to access the corect share and it even sees the printer but it apears it doesn't have permision to print to it. What might I need to do to