similar to: Zeta and Zipf distribution

Dear R users, first I take this opportunity to greet all the R community for your continuous efforts. I wrote a function to calculate the pdf and cdf of a custom distribution (mixed gamma model). The function is the following: pmixedgamma3 <- function(y, shape1, rate1, shape2, rate2, prev) { density.function<-function(x) { shape1=shape1 rate1=rate1 shape2=shape2

Dear R user, I'm using qda (quadratic discriminant analysis) function (package MASS) to classify 58 explanatory variables (numeric type with different ranges) using a grouping variable (factor 2 levels "0" "1"). I'm using the qda method for class 'data.frame' (in this way I don't need to specify a formula). Using the function:

Using R, I plotted a log-log plot of the frequencies in the Brown Corpus using plot(sort(file.tfl$f, decreasing=TRUE), xlab="rank", ylab="frequency", log="x,y") However, I would also like to add lines showing the curves for a Zipfian distribution and for Zipf-Mandelbrot. I have seen these in many articles that used R in creating graphs. Thank you! [[alternative HTML

Dear R-users, I am new to R and would like to use it for fitting the zipf distribution to some numeric data that I have. Here's the snippet that I use: library(VGAM) X <- read.table(file("~\\mydata.txt", encoding="latin1")) w <- as.vector(t((X[2]))) w <- w/sum(w) y <- (1:length(w)) fit = vglm (y ~ 1, zipf, tra=TRUE, weight=w) zipf(N=NULL,

Hi, Is there a Zipf-like distribution RNG in R? Thanks, Weiguang

Hi We have a couple of FSC Computers here which seem to have a very strange bios. I am trying to boot a MS-DOS disk to flash an Scsi Raid controller. Since the computer doesn't have a floppy anymore I tried syslinux from an usb stick. The menu loads fine but I cannot boot disc images with memdisk on this computer. It works fine on others, of course. But now comes the weird part, if I use the

Dear all, Using: library(vegan) data(BCI) mod <- radfit(BCI[1,]) mod RAD models, family poisson No. of species 93, total abundance 448 par1 par2 par3 Deviance AIC BIC Null 39.5261 315.4362 315.4362 Preemption 0.042797 21.8939 299.8041 302.3367 Lognormal 1.0687 1.0186 25.1528 305.0629 310.1281

An embedded and charset-unspecified text was scrubbed... Name: not available URL: <https://stat.ethz.ch/pipermail/r-help/attachments/20091126/7134fc17/attachment-0001.pl>

Hi, I have a question on generating random variables based on zipf-mandelbrot distribution. So when I execute the following lines: ZM = lnre ("zm", alpha = 2/3, B=0.1) zmsample = rlnre (ZM, n =100) zmsample It generates 100 random values based on a zipf-mandelbrot distribution as below. But how do I make sure the generated random number is within the range of 1 - 6000 only? Can I

Dear all, How can I fit a truncated power law to a vector? I can't find a function to do that. If the function provides an AIC, even better. John

You are correct, ubyte 2 simply means 3rd element of the structure. This quantity is always constant (it specifies the field number), so we could have used any signed or unsigned integer type for it instead of ubyte. UByte means we are limited to at most 256 structure fields but it also makes the bytecode representation more compact. (To remove this restriction, we may migrate to UInt in the

Hi everyone. Newbie to statistics. I have 40 matrices of ~400 values. how may I determine whether the distribution follows zips law? response <-sample (1:20,400*4, replace= TRUE) Thank you vry much. -- View this message in context: http://r.789695.n4.nabble.com/A-introductory-question-about-Zips-law-Newbie-to-statistics-tp4521190p4521190.html Sent from the R help mailing list archive at

Dear R users, I would like to employ count data as covariates while fitting a logistic regression model. My question is: do I violate any assumption of the logistic (and, more in general, of the generalized linear) models by employing count, non-negative integer variables as independent variables? I found a lot of references in the literature regarding hot to use count data as outcome,

Sorry for the dumb question, but I cant work out how to do this. Quick version, How can I re-bin a given frequency distribution using new breaks without reference to the original data? Given distribution has integer valued bins. Long version, I am loading a frequency table into R from a file. The original data is very large, and it is a very simple process to get a frequency distribution

Hi all, I would like install 389 Directory Server on a Centos 6.8 server. I followed the instructions reported in : https://wiki.centos.org/HowTos/DirectoryServerSetup?highlight=%28directory%29%7C%28server%29 and installed with yum the following packages : 389-ds-base-1.2.11.15-75.el6_8.x86_64 389-ds-console-1.2.6-1.el6.noarch 389-ds-base-libs-1.2.11.15-75.el6_8.x86_64

Hi I am having problems loading bcdw (http://bootcd.narod.ru) from the isolinux bootloader. I use this to load the acronis programs. It worked fine up till version 3.51 but does not work anymore with 3.52. Here is the code to load the second bootmanager label l kernel /bcdw/bcdwboot.bin Any insights to what might be wrong here? With 3.52 I get Error! Cannot load

Ok guys I get this errors: Network UPS Tools ? UPS driver controller 2.4.1 Network UPS Tools ? Generic HID driver 0.34 (2.4.1) USB communication driver 0.31 No matching HID UPS found my ups.conf is: [ap1000] driver = usbhid-ups port = auto I?ve also tried to set a specific usb port and change driver. My ups is a AMACROX 1000 http://www.amacrox.com//english/01_product/02_01detail.aspx?pdid=46

I have a question about GetElementPtr. Suppose I have an GetElementPtr Instruction GI: %reg = getelementptr %ST* %s, uint 1, ubyte 2, ubyte 1, uint 5, uint 13 I want to check if this is the reference of a component of a structure, how can I do that? Should I check which operand of this instruction is 'ubyte' type? How can I do that in code? should I use ubyte *a =

A little patch that makes the nv50 mesa driver use the correct GPU values for the matching depth formats. There doesn't seem to be Z16, but I added Z32_FLOAT to show the tile_flags the various formats need in the page tables. There'd also be Z32_FLOAT_S8_UNORM (0x19, tile flags 0xe000), but that's not yet defined in the pipe formats header, and I actually don't know if it's

Hi everyone, I am plotting a boxplot with main title as main = bquote(paste(.(ts.ind[s]), ": ", bar(zeta), " Boxplot from 2001 to 2009", sep = "")) but it doesn't work.  The program said they cannot find the function "bar".  Does anyone know how to do it correctly?  Thanks. tin [[alternative HTML version deleted]]