similar to: how to figure out "skewness"

I think there is a question in R. I donot know the reason. This is my data about comulative percentage figure. The result is not right. The first point (no=1,leiji=26.94350) in the plot figure was showen in a lower location. Why? Thanks ! > fre no leiji 1 1 26.94350 2 2 46.86401 3 3 60.59032 4 4 72.17355 5 5 77.85521 6 6 82.05853 7 7 85.56495 8 8 87.64378 9 9

Hello. I came across your response in an R forum and could use your help. I have a data set with 472 rows. I want to delete rows 416 through 472. The name of my data set is MERGE. I am an extreme R novice. How do I write a script to accomplish this? Thank you. --- Christopher H. Porter, M.A., M.Ed. Director, Undergraduate Recruitment College of Engineering and Science Clemson University

R-users, Has anyone implemented (coded) the Rai and Van Ryzin dose-response model from their 1985 Biometrics paper. I'm trying to avoid re-inventing the wheel. If not any hints on coding this would be greatly appreciated. Thanks, -- ==================== Calvin L. Williams, Ph.D. ==================== Department of Mathematical Sciences, Clemson University

Dear R users, I hope to explain the concepts of skewness and kurtosis by generating series of distributions with same skewness and different kurtosis or with same kurtosis and different skewness, but it seems that i cannot find the right functions. I have searched the mailing list, but no answers were found. Is it possible to do that in R? Which function could be used? Thanks a lot. --

Hello everyone, Does anyone know what would cause the skewness() function (from e1071), as well as skew() from psych, to return a value of NaN? I have a vector of positively-skewed data (https://docs.google.com/file/d/0B6-m45Jvl3ZmYzlHRVRHRURzbVk/edit?usp=sharing) which these functions return a value for like normal: > skewness( data ) # returns 1.400405 but when I instead give those

Now I have the forest plot data with x, y locations, and I measured the DBH for every indivicuals. The data looks like that: x=runif(100,0,100) y=runif(100,0,100) dbh=runif(100,1,100) rdata=data.frame(x,y,dbh) > rdata[1:5,] x y dbh 1 99.5354145 1.412844 34.10112 2 0.8259361 87.737036 39.12914 3 6.5678613 65.699032 22.55990 4 67.2987881 72.053877 45.83978 5 2.2491372

R-Users, I'm in the process of developing some notes on optimizing likelihoods for a class I'm teaching next term. Can anyone offer examples showing the optimization of a complex likelihood function that they wouldn't mind sharing using optim(). I've adapted some of the S-plus likelihood examplesfrom the ``Guide to Statistics" and the Zip example from ``Statistical Models in

Hi, is there a test for the H0 skewness=0 (or with skewness as test statistic and normality as H0) implemented in R? Thank you, Christian *********************************************************************** Christian Hennig Fachbereich Mathematik-SPST/ZMS, Universitaet Hamburg hennig at math.uni-hamburg.de, http://www.math.uni-hamburg.de/home/hennig/

I've been using Samba with OpenLDAP with great success on normal servers. Recently however, it appeared to us that for remote locations it is more economically viable to replace Samba servers with Samba running on little routers like ASUS WL-500g with openwrt firmware/software. It has a broadcom/mipsel CPU, and thanks to openwrt (http://openwrt.org), it is possible to run lots of software

I have a dataset giving traffic between pairs of ports, and the lat/lon of each port, roughly as follows: sample data as follows > df2VisitTrips[1:4,] Origin Destination NumberOfTrips OriginLatitude OriginLongitude DestinationLatitude DestinationLongitude 1 P1 P16 1 39.45965 -80.15633 40.76111 -79.54583 2 P1 P3

Hi, in which R-package I could find skewness and kurtosis measures for a distribution? I built some functions: gamma1<-function(x) { m=mean(x) n=length(x) s=sqrt(var(x)) m3=sum((x-m)^3)/n g1=m3/(s^3) return(g1) } skewness<-function(x) { m=mean(x) me=median(x) s=sqrt(var(x)) sk=(m-me)/s return(sk) } bowley<-function(x) { q<-as.vector(quantile(x,prob=c(.25,.50,.75)))

I am trying to reorder a factor data type so that when I plot stats associated with the factor, the ordering makes sense. For instance, if I have a factor entered as follows ... A = as.factor(c("1", "10", "3", "3", "10", "10")) levels(A) ... the ordering does not really make sense (assuming I want the factor ordered by integer

I checked the help and the mailing list archives, but I can find no mention of a routine that calculates higher moments like skewness and kurtosis. Of course, these are easy enough to write myself, but I was thinking that they MUST be in here. Am I wrong? Thanks. -Frank

Hi, Is there a function in R to calculate the coefficient of skewness of some data? I had expected there to be one, but can find no information about it. Thanks for any pointers. Robin Williams Met Office summer intern - Health Forecasting robin.williams@metoffice.gov.uk [[alternative HTML version deleted]]

I have a question about the D'Agostino skewness test and the Anscombe-Glynn kurtosis test. agostino.test(x, alternative = c("two.sided", "less", "greater")) anscombe.test(x, alternative = c("two.sided", "less", "greater")) The option "alternative" in those two functions seems to be the null hypothesis. In the output, the

Is there some reasonable way to generate random data from a distribution that has some degree of skewness and/or kurtosis, but would otherwise be normal? thanks, -------------- next part -------------- A non-text attachment was scrubbed... Name: greiff.vcf Type: text/x-vcard Size: 398 bytes Desc: Card for Warren R. Greiff Url :

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, Can anyone please tell me how to get ".probe" function in "struct xenbus_driver" called? I want to do one simple test of event channel communication between peer modules of domU and dom0. I wrote one backend and one frontend driver (I tried registered as misc/blk/input device) where I put in all my xenbus routines, but cannot get .probe called after insmod my module, I

Hello. There are so many sound files in /var/lib/asterisk/en. Is there an easy way to let me play back all of them one by one while I am watching CLI to see the current file name? Thanks for help. -- Jian Gao IT Technician SJ Geophysics Ltd. <http://www.sjgeophysics.com> jian.gao at sjgeophysics.com <mailto:jian.gao at sjgeophysics.com> Tel: (604)582-1100

Dear all, Right now, I'm doing research about outlier in statistical data (univariate and multivariate) and I want to simulate its behavior. My problem is : How to generate random data from distribution with zero skewness and some kurtosis values in R ? A. Kudus ===================== Dept. of Statistics Bandung Islamic University I n d o n e s i a ==========================