similar to: why doesn't t.test value include standard error

I wish the t.test function in stats would return the standard error. It would be nicer for students if R simply reported the standard error used to calculate the t value. I trolled for this in r-help and got no answers, which I interpreted to mean that this is boring but possibly not wrong. Hopefully. I believe only simple changes are needed. In the source code src/library/stats/t.test.R file:

I'd like to learn the process of revising R functions & packages and then submitting proposed patches to the R Core team. Would someone be willing to mentor me through one example? For starters, consider an example. I'd like to revise the t.test function to return the stderr value to the user. We only need to change the "rval" in the third-from-the end line of

I have a dataframe ("mydf") that contains "differences of means". I wish to test whether these differences are significantly different from zero. Below, I calculate the t-statistic for each column. What is a "good" method to calculate/look-up the p-value for each column? mydf=data.frame(a=c(1,-22,3,-4),b=c(5,-6,-7,9)) mymean=mean(mydf) mysd=sd(mydf)

I'm writing to ask if R Core would make file.exists more Windows tolerant when the argument has a trailing slash. This has been discussed by users a few times here, I know it is not a new topic. But it is not a solved problem, yet. I acknowledge that CRAN packages exist which fix this by replacing file.exists(), but it seems more elegant to me to fix the problem in R itself. R Core goes to

I'm aware R-devel no longer includes zlib. This works find on up-to-date Linux systems. On the older Centos 6 cluster at KU, they have zlib tool old for R-devel. The R-devel configure fails thus: checking if zlib version >= 1.2.5... no checking whether zlib support suffices... configure: error: zlib library and headers are required In other software, I've seen this kind of thing, so

Yesterday I uploaded new rockchalk_1.8.97. Then I received email saying that I needed to confirm the submission. Here's the message. Dear Paul E. Johnson Someone has submitted the package rockchalk to CRAN. You are receiving this email to confirm the submission as the maintainer of this package. To confirm the submission to CRAN, follow or copy & paste the following link into your

This is a basic Windows system administrator problem, asked by a Linux guy who is helping out in a Windows lab. I want to keep R packages up to date on MS Windows 7 with a job in the "Task Scheduler". I have an R program that I can run (as administrator) that updates the existing packages and then installs all the new ones. I do not understand how to run that in a dependable way in

1. One general question for general discussion: Is HLM6 faster than lmer? If so, why? What should I watch for to spot the differences? I'm always advocating R to students, but some faculty members are skeptical. A colleague compared the commercial HLM6 software to lmer. HLM6 seems to fit the model in 1 second, but lmer takes 60 seconds. My first thought was that LM6 uses PQL by default,

I'm attaching a file foo.Rnw and I'm hoping some of you might run it through your R & latex systems to find out if the double-quotes in typewriter font turn out as black boxes (as they do for me). If you don't use Sweave, but you have a system with a working version of R and LaTeX, the file gives the instructions you need to use to process the file. The The file itself explains

Quite possibly, I should keep my big mouth shut and let the CRAN'ers answer this, but at face value, you submit a new version and are asked to confirm that problems with the previous version are fixed; isn't that as it should be? -pd On 08 Feb 2016, at 16:10 , Paul Johnson <pauljohn32 at gmail.com> wrote: > Yesterday I uploaded new rockchalk_1.8.97. Then I received email >

Hi Paul, This might be what you are looking for. https://stat.ethz.ch/pipermail/r-devel/2015-April/070951.html Jesper On Fri, Feb 12, 2016 at 6:13 PM, Paul Johnson <pauljohn32 at gmail.com> wrote: > I'm aware R-devel no longer includes zlib. This works find on > up-to-date Linux systems. > > On the older Centos 6 cluster at KU, they have zlib tool old for > R-devel.

> On Mar 31, 2016, at 1:00 PM, Paul Johnson <pauljohn32 at gmail.com> wrote: > > In the rockchalk package, I want to provide functions for regression > objects that are "well behaved." If an object responds to the methods > that lm or glm objects can handle, like coef(), nobs(), and summary(), > I want to be able to handle the same thing. > > It is more

Hello, everybody: I'm experimenting more with Sweave, R, and LyX. There's now an entry in the LyX wiki on using R, so anybody can do it! http://wiki.lyx.org/LyX/LyxWithRThroughSweave Now I notice this frustrating thing. I think I've done everything possible to make the Design library start up without any fanfare: <<echo=F,quiet=T,print=F>>= library(Design, verbose=F)

I set up a local repo for testing packages. My packages are not showing up from the repository when viewed by Linux clients. I suspect this is a web administrator/firewall issue, but it could be I created the repo wrongly. I am supposed to run write_PACKAGES separately in each R-version folder. Right? Maybe other novices can use these scripts, if they are not wrong :) Here's the file

Here's my little discussion example for a quadratic regression: http://pj.freefaculty.org/R/WorkingExamples/regression-quadratic-1.R Students press me to know the benefits of poly() over the more obvious regression formulas. I think I understand the theory on why poly() should be more numerically stable, but I'm having trouble writing down an example that proves the benefit of this. I

I'm writing to ask for a policy exception, or advice on how to make this package CRAN allowable. http://rweb.quant.ku.edu/kran/src/contrib/portableParallelSeeds_0.9.tar.gz Yesterday I tried to submit a package on CRAN and Dr Ripley pointed out that I had not understood the instructions about packages. Here's the part where the R check gives a Note * checking R code for possible

I'm looking for an R mentor. I want to propose a change in management of plot options xlim and ylim. Did you ever want to change one coordinate in xlim or ylim? It happens to me all the time. x <- rnorm(100, m=5, s=1) y <- rnorm(100, m=6, s=1) plot(x,y) ## Oh, I want the "y axis" to show above x=0. plot(x,y, xlim=c(0, )) ##Output: Error in c(0, ) : argument 2 is empty

I'm making some illustrations and it would be convenient to automatically shade the overlapping portions of circles. These illustrations are for Social Choice theory, a field in political science and economics. I've wrestled together some examples so you can see what I mean, but have not mastered the "color overlapping sections" problem (as you will see):

Suppose plotx <- "someName" modx <- "otherName" plotxRange <- c(10,20) modxVals <- c(1,2,3) It often happens I want to create a dataframe or object with plotx or modx as the variable names. But can't understand syntax to do that. I can get this done in 2 steps, creating the data frame and then assigning names, as in newdf <- data.frame( c(1, 2, 3, 4),

I hope you are doing well. For me, this was an unexpected problem. I've hoped for quite a few wrong things today, but I'm only asking you about this one. Why does ifelse(1, list(a, b, c), list(x, y, z)) return a list with only a, not list(a, b, c) as I hoped. I wish it would either cause an error or return the whole list, not just the first thing. Working example: > x <- 1 >