similar to: estimate alpha beta for gamma

I apologize that I've asked a similar question before, but being new to R I don't think I did a very good job of formating the question. I've included a text file since the date set is somewhat large. What I have is a huge string of numbers in a text file. The numbers are all separated by comma's and the groups are separated by a semicolon. What I would like to do is read each

Hi! I have 2 data.frames: "fish" and "popn": >fish xloc yloc id birth size weight energy gonad 20 15 15 54 -60 107.9 63.0 15952.9 8.0 21 15 15 32 -60 105.1 61.4 15538.8 7.8 91 4 43 96 -60 118.9 69.4 17573.2 8.8 71 32 4 64 -60 121.6 71.0 17976.0 9.0 34 2 64 20 -60 116.2 67.9 17173.0 8.6 95 6 20 58 -60 106.5

Dear r-users, I have dailly rainfall data from year 1971 to 2000. I use aggregate to form monthly rainfall data.  What I don't understand is that the data for the year 2000 become on the top, instead of year 1971.  Here are some codes and output: agg_dt1     <- aggregate(x=dt1[,4],by=dt1[,c(1,2)],FUN=sum) > head(agg_dt1,20); tail(agg_dt1,20)    Tahun Bulan     x 1      0     1 398.6

Hi all, I have this figure: http://docs.google.com/Doc?id=df5zfsj4_103rjt2v4d5 created with the following steps: > x [1] 90.4 57.8 77.0 103.7 55.4 217.5 68.1 85.3 152.0 113.0 97.1 89.9 [13] 68.1 83.7 77.4 34.5 104.9 170.3 88.6 88.1 108.8 77.4 85.6 82.7 [25] 81.3 108.0 49.5 71.0 85.7 99.3 203.5 275.9 51.1 84.8 16.5 72.6 [37] 160.5 158.3 136.7 140.0 98.4 116.1

HI All R Gurus, Please reply directly to this email if you have a suggestion, as I am not currently an R-help member. I would like to complete a non-linear regression with groups. Having gone down the path of completing such things under GENSTAT, I would prefer to do it with R. My example data set is... TimeTreat ExtrVol2 SorbedT ET1 39.5 384.5 ET1

Hello Everybody, I am running R 2.2.0 with Windows XP i am trying to fit nonlinear differential equation to data sets which looks like this: Week N C 0 1 1 1 5 6 2 6.2 12.2 3 59 71.2 4 39 110.2 5 38 148.2 6 44 192.2 7 20.4 212.6 8 19.4 232 9 34.2 266.2 10 35.4 301.6 and i need to fit these data to the following diff equation: dNdt=a*N-b*N*C, dCdt=N^2, Where a=birth rate, b=death rate and N=

I'm using lmer to fit mixed-effect logistic regression models. This is for a small data set. First, I fit a constant: Generalized linear mixed model fit using Laplace Formula: propm ~ (1 | study) Data: inducedSR71507.dat Family: binomial(logit link) AIC BIC logLik deviance 183.7 189.4 -89.84 179.7 Random effects: Groups Name Variance Std.Dev. study (Intercept)

Dear all, I am trying to fit a generalized linear mixed model to deal with both temporal and spatial pseudorreplication. I have repeated seasonal measurements (3 seasons during 3 years and 2 season during the last year, factor named percod, with 11 levels) of a bird presence/absence (pres.f) on 14 artificial ponds (charca.f). The ponds are integrated in 4 greater spatial units (zepa). I have

Hi r-users,   I have a data below but I would like to rearrange my data so that it will be in the following format: year, month, day, rain.  I also would like the data to be arranged so that it will be continues for example year 1939 it will go from january 1 to december 31 and continue with year 1940 from january 1st to 31st December and so on.  What I have now is all january for all the years

`grDevices::convertColor` performance can be improved by 30-300x with small changes to the code. `colorRamp(space='Lab')` uses `convertColor` so it too benefits from substantial performance gains. `convertColor` vectorizes explicitly over the rows of the input color matrix using `apply`. The level-1 patch [1] illustrates a possible minimal set of changes to achieve this with just R

Hola a todos: Necesito pegar 4 bases de datos de N variables y N columnas, EL PROBLEMA ES QUE 3 CONSTAN DE 13 COLUMNAS, y LA CUARTA CON 11 COLUMNAS. He usado rbind y smartbind, sin embargo el resultado no es el esperado, ya que coloca las columnas a un costado y no hacia abajo como debería Como resultado necesito algo como LO QUE PRESENTO ABAJO, RELLENANDO LOS ESPACIOS VACIOS CON "NA"

Dear R users, I have more of a statistical/econometric question than straight R one. I have a data set with the discrete hazard rate of small firms survival on 400 counties over a period of 9 years. This data was generated using census information from the VAT registration number of each one of these business. I would like to analyze the effect of regional factors (deprivation index, real

Hi, I have been doing some tests using hping2 and TCP SYN pings targeting local and remote hosts from two FreeBSD 4.7-RELEASE-p10 and one Linux 2.4.18 host. The three machines have the same hardware configuration and have been running for 6 months now. The average load isnt too high (usually 0.01 to 0.15) on the FreeBSD machines. Here is the output from hping2 (excuse me the line wrap),

https://bugs.freedesktop.org/show_bug.cgi?id=38016 Summary: nouveau driver crashes xorg when starting KDE Product: xorg Version: git Platform: Other OS/Version: All Status: NEW Severity: normal Priority: medium Component: Driver/nouveau AssignedTo: nouveau at lists.freedesktop.org

Dear Rusers, I have met a difficult problem on explaining the differences of principal component analysis(PCA) between R,S-PLUS and SAS/STATA/SPSS, which wasn't met before. Althought they have got the same eigenvalues, their coeffiecients were different. First, I list my results from R,S-PLUS and SAS/STATA/SPSS, and then show the original dataset, hoping sb. to try and explain it.

Hola Carlos: Tanto rbind(base) como rbind(data.table) generan el mismo resultado que no busco. Que solución se le podría dar específicamente a la data de 11 columnas, que es la que me da problema. Saludos. [Descripción: FIRMA2] De: Carlos Ortega [mailto:cof en qualityexcellence.es] Enviado el: martes, 14 de febrero de 2017 14:21 Para: Javier Valdes Cantallopts (DGA) CC: r-help-es en

De esa forma me sigue entregando lo mismo, multiplica las variables por dos (80 obs 26 variables) debería ser (80 obs 13 variables).En concreto pega la información pero la costado, no hacia abajo. Solucionara el tema si le agrego dos columnas ficticias para que queden todas las bases con 13 variables? Saludos. [Descripción: FIRMA2] De: Carlos Ortega [mailto:cof en qualityexcellence.es] Enviado

Both get.hist.quote, and its derivative priceIts, rely on download.file() to fetch financial data series from Yahoo! in .csv format. They allow for nice interactive demonstrations of what one can do with R. Unfortunately, both are currently broken as Yahoo! decided to add a somewhat useless html comment at the end of the csv 'stream', breaking the regular format of n rows with k columns.

Dear R-Users, I can't understand the behaviour of quasibinomial in lmer. It doesn't appear to be calculating a scaling parameter, and looks to be reducing the standard errors of fixed effects estimates when overdispersion is present (and when it is not present also)! A simple demo of what I'm seeing is given below. Comments appreciated? Thanks, Russell Millar Dept of Stat U.

Hello zfs-discuss, One of a disk started to behave strangely. Apr 11 16:07:42 thumper-9.srv sata: [ID 801593 kern.notice] NOTICE: /pci at 1,0/pci1022,7458 at 3/pci11ab,11ab at 1: Apr 11 16:07:42 thumper-9.srv port 6: device reset Apr 11 16:07:42 thumper-9.srv scsi: [ID 107833 kern.warning] WARNING: /pci at 1,0/pci1022,7458 at 3/pci11ab,11ab at 1/disk at 6,0 (sd27): Apr 11 16:07:42 thumper-9.srv