similar to: (PR#7309) misuse of R-bugs (was inappropriate definition

Folks, I'm a bit puzzled by the fact that if I generate 100,000 standard normal variates using rnorm() and perform the Jarque-Bera on the resulting vector, I get p-values that vary drastically from run to run. Is this expected? Surely the p-val should be close to 1 for each test? Are 100,000 variates sufficient for this test? Or is it that rnorm() is not a robust random number generator?

HI Philippos, Try this: dat1<- read.csv("Validation_data_set3.csv",sep=",",stringsAsFactors=FALSE) #converted to csv str(dat1) #'data.frame':??? 12573 obs. of? 17 variables: # $ Removed.AGC????????????????????????????? : num? 65.67 46.17 41.26 14.09 5.38 ... # $ Removed.SST????????????????????????????? : chr? "" "46.1658" "41.2566"

Hi, I would appreciate if someone could help me on track with this problem. I want to compute some parameters from a system of equations given a number of sample observations. The system looks like this: sum_i( A+b_i>0 & A+b_i>C+d_i) = x sum_i( C+d_i>0 & C+d_i>A+b_i) = y sum_i( exp(E+f_i) * ( A+b_i>0 & A+b_i>C+d_i) = z A, C, E are free variables while the other

Hello everybody, this is a bit off topic, so maybe you can just reply to me personally... What is the typical notation for 'skewness', 'kurtosis' and maybe 'excess kurtosis'? Thank you, Martina -------------------------------------------------------------------------- Department of Statistics Office Phone: (614) 292-1567 1958 Neil Avenue, 304E Cockins Hall

Hi, I'm a newcomer to R, having previously used SPSS. One problem I have run into is computing kurtosis. A test dataset is here: http://www.whinlatter.ukfsn.org/2401.dat > library(moments) > data <- read.table("2401.dat", header=T) > attach(data) > loglen <- log10(Length) With SPSS, I get Skewness -0.320 Kurtosis -1.138 With R: > skewness(loglen) [1]

Dear R-Users and Developpers, Currently R does not include functions to compute the skewness and kurtosis. I programmed it myself in the following way, but probably *real* programmers/statisticians can do that better: mykurtosis <- function(x) { m4 <- mean((x-mean(x))^4) kurt <- m4/(sd(x)^4)-3 kurt } myskewness <- function(x) { m3 <- mean((x-mean(x))^3) skew <-

Greetings. For R gurus this may be a no brainer, but I could not find pointers to efficient computation of this beast in past help files. Background - I wish to implement a Cramer-von Mises type test statistic which involves double sums of max(X_i,Y_j) where X and Y are vectors of differing length. I am currently using ifelse pointwise in a vector, but have a nagging suspicion that there is a

Dear useRs, I was searching CRAN for implementation of kurtosis and skewness tests, and found that there is some kind of lack on it. So, I have written three functions: 1. Anscombe-Glynn test for kurtosis 2. Bonett-Seier test based on Geary's kurtosis (which is not widely known, but I was inspired by original paper describing it, found coincidentally in Elsevier database) 3.

Greetings Everybody: I generated a 1.2MB dta file based on the general social survey with Stata8 for linux. The file can be re-opened with Stata, but when I bring it into R, it says all the values are missing for most of the variables. This dataset is called "morgen.dta" and I dropped a copy online in case you are interested http://www.ku.edu/~pauljohn/R/morgen.dta looks like this

Hello I do not know very much about statistics (and English language too :-( ), then I come in search of a clarification (explanation): I found two distinct results on KURTOSIS and I do not know which of them is the correct one. Any aid will be welcome! klebyn ################ CODE rnorm(1000) -> x library(moments) kurtosis(x) skewness(x) detach("package:moments")

We use the tar option of smbclient (Samba 1.9.16p9, SunOS 4.1.3) to backup Windows 95 PCs: smbclient "\\${server}\backup" -U COMPUTING -d 1 -D "\" -Tc tarfile >! logfile backup is set as a share of c:\ The backups are incomplete (as shown by the logfile and gzip -tf tarfile). Typically, only very roughly 5700 of 6600 files get copied. I expect I've got a bad

Hi I'm becoming a bit confused. I have been trying out the Samba PDC code from BRANCH_NTDOM on Solaris 2.5.1 (with great success !). I now want to download the latest version and messages on the samba-ntdom list say to download from the main branch. What exactly do I put in my cvs command ? The NTDOM FAQ says 'cvs -d :pserver:cvs@samba.anu.edu.au:/cvsroot co samba '. When I do

Hi Samba people I saw this as a problem in the archives about a month ago but no solution ? I'm running samba 2.2.10 on Solaris 8. We have a script to replace /usr/bin/passwd which changes the NIS passwd and then uses smbpasswd to change the user's smb password to the same thing. Some users are getting : NIS password changed for diane on cauchy machine 127.0.0.1 rejected the password

Hi I have successfully joined my Sun NAS device to my Samba PDC (version 2.2.10, security=user, encrypt passwords=yes) However, when I try to map a drive to \\nas\home it wont accept my (correct) password and the NAS log says : NetrSamlogon[BIOSS\janet]: SAMBA_ACCESS_DENIED (Samba PDC) Have I missed something ? Do I have to set up a trust account or is it to do with encrypted passwords ? --

Hi I have a Samba server running 1.9.18p7 on a Solaris 2.5 machine. It is set up with a netlogon share. I have 4 Win95 PCs for which I have created identical <machinename>.bat files in the netlogon directory. However, only one of those PCs is running the 'root preexec' script and running the commands in %m.bat. Anyone have any idea why ? I've appended the global and netlogon

Hi If I add a specific machine to my hosts allow list in the global section of my smb.conf, then that machine can access any of the shares on that server. If I only add the machine to a hosts allow list in a specific share, when I try and connect to that share from that machine I get an error : 'Specified computer did not recieve your request' and in log.smb : 'Connection denied from

I wonder whether the functions for skewness and kurtosis in the e1071 package are based on correct formulas. The functions in the package e1071 are: # -------------------------------------------- skewness <- function (x, na.rm = FALSE) { if (na.rm) x <- x[!is.na(x)] sum((x - mean(x))^3)/(length(x) * sd(x)^3) } # -------------------------------------------- and #

Hi I managed to get my Sun NAS filer to join my Samba domain (2.2.10). However, I cannot mount any of the shares from the NAS to a PC in the same domain, getting an error : 06/21/05 14:15:53 NetrSamLogon[BIOSS\janet]: SAMBA ACCESS_DENIED (Samba PDC) When I asked Sun about it they responded : "Further reseach has shown that a connection to SAMBA domain controller as is not supported by the

Hi, I have the following function: > kurtosis <-function(x) (mean((x-mean(x))^4))/(sd(x)^4) #x is a vector and data > print(mydata) V1 V2 V3 V4 V5 1 1007_s_at DDR1 2865.1 2901.3 1978.3 2 1053_at RFC2 103.6 81.6 108.0 3

Thanks to all who response my problem. Here are my summary : 1. from Dirk Eddelbuettel <edd at debian.org> We could try a mixture of normals -- ie flip a coin (use a uniform with some cutoff c where 0 < c < 1 ) to choose between N(0, sigma_1) and N(0, sigma_2). 2. from Michaell Taylor <michaell.taylor at reis.com> We could use the gld library to specify the lambdas of