similar to: Standard error of standard deviation: bootstrap or theoretical results?

I am wanting to plot a 95% confidence band using segplot, yet I am wanting to have groups. For example if I have males and females, and then I have them in different races, I want the racial groups in different panels. I have this minor code, completely made up but gets at what I am wanting, 4 random samples and 4 samples of confidence, I know how to get A & B into one panel and C&D in to

Hi, I am new to R. My problem is with the ordered logistic model. Here is my question: Generate an order discrete variable using the variable wrwage1 = wages in first full calendar quarter after benefit application in the following way: * wage*1*Ordered *= 1 *if*0 *· wrwage*1 *< *1000 2 *if*1000 *· wrwage*1 *< *2000 3 *if*2000 *· wrwage*1 *< *3000 4 *if*3000 *· wrwage*1 *<

I am using JRI (Java R Interface) library in order to call R from within my Java application. But since the "rmu1" and "rmu2" ,see the following code, are objects of type S4 once i run the application the value of Null will be returned for both of them. On this regard, i would appreciate it if anyone can tell me how i am going to display and/ or convert these objects to Java

Please help me in writing the R code for this problem. I've been solving this for 4 days. It was hard for me to solve it. It's a simulation problem in R. The problem is My true model is a normal mixture which is given as 0.5 N(-0.8,1) + 0.5 N(0.8,1). This model has two components. I will get a random sample of size 100 from this model. I will do this 300 times. That means, I will have

I continue to struggle with mixed models. The square zero version of the problem that I am trying to deal with is as follows: A number (240) of students are measured (tested; for reading comprehension) on 6 separate occasions. Initially (square zero) I want to treat the test time as a factor (with 6 levels). The students are of course ``random effects''. Later I want to look at

Dear all, I try to fit least absolute deviation (LAD) regression with R. I couldn't find any command to carry out LAD. What I am doing now is using nlm() and setting the function argument to minimize the sum of the absolute deviation. But for some data it keeps giving out the error message: Error in nlm(absdev, b.start) : non-finite value supplied by nlm. (there is no na in my dataset.) I

Dear R users, I want to retrieve "..." argument values within a function. Here is a small exmaple: myfunc <- function(x, ...) { if (hasArg(ylim)) a <- ylim plot(x, ...) } x <- rnorm(100) myfunc(x, ylim=c(-0.5, 0.5)) Error in myfunc(x, ylim = c(-0.5, 0.5)) : Object "ylim" not found > I need to retrieve values of "ylim" (if it is defined

Thank you for helping me! I try the "pretty" funtion to select the x-axis position value.Then I use the "format" funtion. xax.pos <- pretty(as.numeric(x$x.name)) format(xax.pos,'%d %b %y') > xax.pos [1] 1091600000 1091800000 1092000000 1092200000 1092400000 1092600000 1092800000 [8] 1093000000 There are something wrong. I found the xax.pos has been changed to

UPDATE: It turned out QEMU still dumped core with that patch applied. That said, I'm still interested to know how to report issues with QEMU from centos-virt repo. Thanks, Ray > -----Original Message----- > From: Xiong, Huan > Sent: 2018?6?22? 11:12 > To: 'centos-virt at centos.org' <centos-virt at centos.org> > Subject: Q: how can I request to add an upstream

I am trying to use Fisher's z' transformation of the Pearson's r but the standard error does not appear to be correct. I have simulated an example using the R code below. The z' data appears to have a reasonably normal distribution but the standard error given by the formula 1/sqrt(N-3) (from http://davidmlane.com/hyperstat/A98696.html) gives a different results than sd(z). Can

Hi everybody, I have trouble installing my own package under R 2.1.0 (it is fine under R 2.0.1). My OS is Windows NT. I installed my package 'mag' by using menu "Packages/Install package from local zip files....". It's fine (html package description was updated). But when I typed in library(mag) it gave me error: Error in library(mag) : there is no package called

I am using GLM to calculate logit models based on cross-sectional data. I am now down to the hard work of making the results intelligible to very average readers. Is there any way to calculate a psuedo analoque to the R^2 in standard linear regression for use as a purely descriptive statistic of goodness of fit? Most of the readers of my report will be vaguely familiar and more comfortable with

Hi everyone, I want to use the Fisher r-to-z transformation as part of a hypothesis test of r. I can't find an R function that can do that. Am I missing it? -Tim -- Tim Wilson | Visit Sibley online: | Check out: Henry Sibley HS | http://www.isd197.org | http://www.zope.com W. St. Paul, MN | | http://slashdot.org wilson at visi.com | <dtml-var

----- Original Message ----- > From: "Zhang Huan" <zhanghuan at open-fs.com> > To: "Raghavendra G" <raghavendra at gluster.com> > Cc: "GlusterFS Maintainers" <maintainers at gluster.org>, "Gluster Devel" <gluster-devel at gluster.org>, "Kaushal Madappa" > <kmadappa at redhat.com> > Sent: Tuesday, May 30,

Hello. I have an XYY series that I would like to graph with matplot() or some other single function that will do the trick. The X in question is a vector of POSIXt values obtained from strptime(). Is it possible to tell matplot() how to handle POSIXt x values? I have examined the examples at http://lark.cc.ukans.edu/~pauljohn/R/statsRus.html#5.22 , but would prefer not have to overlay the

Hi everybody, I have trouble in installing packages from source code by following Section 5.1 in manual R-admin.pdf . I am using R 2.1.0 and Win NT. Following the Windows toolset section in the manual, I download the tool set package from: http://www.murdoch-sutherland.com/Rtools/tools.zip and unzip under C:\tools I also downloaded Perl (Windows Port) and installed it.

Hi Folks, I wonderif someone who is familiar with the details of vglm in the VGAM package can assist me. I'm new to using it, and there doesn;t seem much in the documentation that's relevant to the question below. Say I have a vector x of 0/1 responses and another vector y of 0/1 responses, these in fact being a bivariate set of 0/1 responses equivalent to cbind(x,y). E.g.

Hi All, As the in-place help said, --in-place Do not create an output virtual machine in the target hypervisor. Instead, adjust the guest OS in the source VM to run in the input hypervisor. This mode is meant for integration with other toolsets, which take the responsibility of converting the VM configuration, providing for rollback in case of errors, transforming the storage, etc. I want

Full_Name: Huan Huang Version: 1.9.1 OS: Win NT Submission from: (NULL) (155.140.122.227) Traceback: Traceback (most recent call last): File "/scratch/local/app/mailman-sfs/2.1.5/scripts/driver", line 87, in run_main main() File "/scratch/local/app/mailman-sfs/2.1.5/Mailman/Cgi/options.py", line 91, in main language = cgidata.getvalue('language') File

Hi! I am calculating correlation between two variables: 1. X versus Y 2. X versus Y(with a 3 steps lag) I would like to test if the correlation increase/decrease from 1 to 2 is significant or not. Is there any function in R to do this? any hints? Thanks for help :) David Ria?o Center for Spatial Technologies and Remote Sensing (CSTARS) University of California 250-N, The Barn One Shields