similar to: calculation of the p value in ks.test()

Hi, I'm using the ks.test() on two vectors. I looked up the reference and also coded up a version of the two sample Smirnov test. My question is that how can I decide from the output of R that the two vectors x & y come from the same distribution? Am I correct in assuming that smaller D values indicate that they come from the same distribution? In addition how can I use the p value that

Hi, is there a way to plot distribution curves (say normal or chi sq etc) from within R? For example I looked up the *chisq family of functions but I'm not sure as to how I would use them to generate a plot of the chi sq distribution (for arbitrary d.o.f). Thanks, ------------------------------------------------------------------- Rajarshi Guha <rajarshi at presidency.com>

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hi, I'm a new user of R and have managed to make a plot of a histogram. Is there any way I set the title and axes labels and then save the plot as an image (png/gif)? Thanks - -- - ------------------------------------------------------------------- Rajarshi Guha <rajarshi at presidency.com> <http://jijo.cjb.net> GPG

Hi, I'm plotting 4 graphs on one page (2x2 matrix) but I cant seem to get the title for the whole page right. I'm doing: op <- par(mfrow = c(2,2), pty="s") hist(var$V2, breaks="FD",main="Euclidean Metric", xlab="Sum of 3NN ... hist(var$V2, breaks="FD",main="Manhattan Metric", xlab="Sum of 3NN ... hist(var$V2,

I'm running Debian/Jessie on a server which is used as an AD DC. I'm getting errors like this: Hit http://ftp.ca.debian.org jessie InRelease Hit http://ftp.ca.debian.org jessie/main Sources Hit http://ftp.ca.debian.org jessie/main amd64 Packages Hit http://ftp.ca.debian.org jessie/main Translation-en Err http://ftp3.nrc.ca jessie InRelease Err http://ftp3.nrc.ca jessie-updates InRelease

Hi All, I am trying to use A Gaussian quadrature over the interval (-infty,infty) with weighting function W(x)=exp(-(x-mu)^2/sigma) to estimate an integral. Is there a way to do it in R? Is there a function already implemented which uses such weighting function. I have been searching in the statmode package and I found the function "gauss.quad(100, kind="hermite")" which uses

and 'host google.com 8.8.8.8' also gives the expected output? Like this: > root at pf:~# host google.com 8.8.8.8 > Using domain server: > Name: 8.8.8.8 > Address: 8.8.8.8#53 > Aliases: > > google.com has address 74.125.136.100 > google.com has address 74.125.136.113 > google.com has address 74.125.136.101 > google.com has address 74.125.136.138 >

Hello to all, I am having trouble with intregrating a complicated uni-dimensional function of the following form Phi(x-a_1)*Phi(x-a_2)*...*Phi(x-a_{n-1})*phi(x-a_n). Here n is about 5000, Phi is the cumulative distribution function of standard normal, phi is the density function of standard normal, and x ranges over (-infty,infty). My idea is to to use quadrature to handle this integral. But

I have a function written for Splus, when I run it in R I obtain get an error because the function has the elements "0.d0" and "2.d0". How can I change it to run in R? The function can be found in page 230 from http://www.stat.wisc.edu/~mchung/teaching/stat471/stat_computing.pdf Function is as follows: gauher <- function(n) {# Gauss-Hermite: returns x,w so that

Full_Name: Robin Hankin Version: 2.4.0 RC OS: MacOSX 10.4.7 Submission from: (NULL) (139.166.242.29) The first cut line described in Trig.Rd for asin() is incorrect in the ascii version of the manpage. The Rd file reads: For \code{asin()} and \code{acos()}, there are two cuts, both along the real axis: \eqn{\left(-\infty, -1\right]}{\(-Inf, 1\]} and Note the inconsistency between the

-----BEGIN PGP SIGNED MESSAGE----- Hash: SHA1 Hi, First, as far as I can tell, this topic hasn't been mentioned in the mailing lists in awhile (atleast not as far as google can tell). Is there any client/program that can send "net start" requests to windows boxes. The ability to start a program as a service on a windows box remotely would obviously be very handy. I am aware of

Does rsync -aAXv /dir1 /dir2, where dir2 is a mounted FAT formatted SD card, and then: rsync -aAXv /dir2 /dir1, to restore the files in dir1 from the FAT formatted SD card, restore the ACL and owner/permissions of the files and directories in dir1? Thanks, John -- John Conover, conover at rahul.net, http://www.johncon.com/

Greetings, Sorry, the last message was sent by mistake! Here it is again: I encounter a strange problem computing some numerical integrals on [0,oo). Define $$ M_j(x)=exp(-jax) $$ where $a=0.08$. We want to compute the $L^2([0,\infty))$-inner products $$ A_{ij}:=(M_i,M_j)=\int_0^\infty M_i(x)M_j(x)dx $$ Analytically we have $$ A_{ij}=1/(a(i+j)). $$ In the code below we compute the matrix

Hi, I'm running a calculation in two ways. The first way is to employ vectors and evaluate a function in one go. The second way is to break down the function into pieces and combine the pieces to the final answer. Algebraically, they should give me the same result. But the final vector differs significantly. I must be missing something very obvious, but I just cannot see it xx <-

Hi there. Thanks for your time in advance. I am using R 2.2.0 and OS: Windows XP. My final goal is to calculate 1/2*integral of (f1(x)^1/2-f2(x)^(1/2))^2dx (Latex codes: $\frac{1}{2}\int^{{\infty}}_{\infty} (\sqrt{f_1(x)}-\sqrt{f_2(x)})^2dx $.) where f1(x) and f2(x) are two marginal densities. My problem: I have the following R codes using "adapt" package. Although "adapt"

We have a Samba-4.3 AD-DC hosted on a FreeBSD-10.3 system. We have just discovered that for some reason our domain workstations are not synchronising their time with the AD-DC. We have not altered the default GPO nor have we created a special GPO to handle time services. I have been working on several of the domain clients and I cannot seem to get them to source from the Samba server. I have

Dear list, I'm seeking some advice regarding a particular numerical integration I wish to perform. The integrand f takes two real arguments x and y and returns a vector of constant length N. The range of integration is [0, infty) for x and [a,b] (finite) for y. Since the integrand has values in R^N I did not find a built-in function to perform numerical quadrature, so I wrote my own after

Hi, I had some Java code that worked with SJava 0.65 under R 1.8.*. I'm trying to get it to work with R 2.0.0. My JVM is Sun 1.5.0 (running on Fedora Core 2) I downloaded and installed SJava 0.65 in my personal directory using R CMD INSTALL SJava-0.65.tar.gz -l ~/src/Rlibrary and before starting R I sourced RJava.bsh. However after loading the SJava library if do javaConfig() I get

Hi, I have 3 lists, x, y, z and I'd like to perform a calculation over all the lists simultaneously. If it were a single list I could use lapply, but for more than one list I'm using a for loop. Is there an idiom that would let me use something like lapply, but the function specified to lappy would have access to an element from each list? (In Python, I would have used for a,b,c in

>>>>> Jan van der Laan <rhelp at eoos.dds.nl> >>>>> on Fri, 6 Oct 2017 12:13:39 +0200 writes: > It is actually model.matrix that crashes, not glm. Same > crash occurs with e.g. lm. > model.matrix(dob_mon ~ dob_day*dob_mon, data = tab) > also crashes R. Yes, segmentation fault. It only happens when these are *logical*