similar to: Non-linear random variables

Hi All (first timer here, be gentle!!) I'm trying to extract the residual variance out of a glm with poisson data to conduct power and sample size to detect a change over years. I have equal reps within year (for now). Is the residual variance just the grand mean of the original count data, or something trickier? glm.Male_glm(M~year*strata,family=poisson,offset=dur) levels(year) [1]

> Hi All (first timer here again) > I've sent this message once before but I'd like to add a bit more. > I'm trying to conduct power and sample size to detect a change over years > from poisson data > Is there some simulation routine in R that I could use (or that someone has used eslwhere) on a set of generated data that has a mean, say 10% less than current

Hi, I'm new to R and am trying to implement a state-space model, and do some simulations looking at how well it works. Does anyone have an example I could use as a template? I'm having trouble working out where to start, and haven't been able to find anything similar in the help files or on the web. My example is a surplus production fisheries model, so something similar to that

Dear R-list members, Cynthia M. Jones wrote a paper (Fitting growth curves to retrospective size-at-age data, Fisheries Research 46(2000):123-129; abstract at http://www.elsevier.nl/gej-ng/10/19/44/70/24/37/abstract.html)where the SAS procedure MIXED, Macro NLINMIX (Littell et. al., 1996)was used to estimate the von Bertalanffy growth function parameters assuming that data from the same fish are

I would like to calculate a first and second derivative and am having problems finding a simple solution. My syntax may be off as I am not a mathematician, so pardon ahead of time. data: t<-seq(0,4, by=.1) The function is: H(t) = t^3-6*t^2+5*t + 30 from here I plot the curve: plot(x,y ,xlab="x-values", ylab="f(x)", type="l") But would like to similarly plot the

Hello I am trying to fit a REML to some soil mineral data which has been collected over the time period 1999 - 2004. I want to know if the 19 different treatments imposed, differ in terms of their soil mineral content. A tree model of the data has shown differences between the treatments can be attributed to the Magnesium, Potassium and organic matter content of the soil, with Magnesium being the

Hola a tod en s, estoy trabajando con modelos de regresión mixta con el paquete "nlme", concretamente con la función "lme". Una vez que hago mi modelo y el summary del mismo obtengo por ejemplo la significación (valor p) de la variable independiente, así como el valor del AIC para comparar mi modelo con otro similar, pero no sé cómo obtener un valor del grado de ajuste de la

I went to RedHats' page and found hardware requirements and it said that the Dell Poweredge 6300 is a supported hardware. Then I went to linuxcare.com and it does not list the machine. Before I look like a buffoon, or continue to, Redhat will run on this machine, right? I appreciate the help. Thanks, --------------------------------------------------- Tony Ricker Technology Coordinator

All, I have the following in my smb.comf for printing commands... # NOTE: If you have a BSD-style print system there is no need to # specifically define each individual printer [printers] comment = All Printers path = /var/spool/samba browseable = no print ok = yes print command = /usr/bin/lpr -P%p -r %s lpq command = /usr/bin/lpq -P%p lprm command = /usr/bin/lprm -P%p %j

I am using Rsync between a Redhat Linux box and an AIX RS600. We have a about 30gb of database we need to sync to a backup server. Sounds good, right? The problem is that Rsync is so slow when we do the initial dump. We have files that are 1 - 5gb. It takes around 14-20 hours to Rsync the file structure, which means our daily backup will miss that days data because we are waiting for the data

hai, I want to perform a simulation like this仯伜 Suppose that I have one population ,it's size is 500, is composed of x,y and z. The probability of x, y and z is respectively is 0.3, 0.5, 0.2. I wan to simulate a new same size population based ratio of x, y and z, how can I get and assess the number of x, y and z. ? luan Key Laboratory for Sustainable Utilization of Marine Fisheries

Dear Alberto, There are fisheries people also in Europe using AD Model Builder (Denmark and England for instance), but you are probably right that it is more widespread in North America. There is also effort going on where people try to make assessment models written in ADMB callable from R. best regards, hans > think AD Model Builder is mainly used for fisheries assessment in North

There has recently been some discussion on the list about AD Model builder and the suitability of R for constructing the types of models used in fisheries management. https://stat.ethz.ch/pipermail/r-help/2006-January/086841.html https://stat.ethz.ch/pipermail/r-help/2006-January/086858.html I think that many R users understimate the numerical challenges that some of the typical

Greetings! I am trying to compare simulated environmental conditions from a model against a recruitment time series for a species of crab by first dropping 5 data points, and then using the remainder to attempt to simulate the missing data as a measure of best fit and using the following code: all.mat<-as.matrix(comb,ncol=ncol(comb),nrow=nrow(comb)) obs<-as.matrix(R2,24,1) mod<-all.mat

I am running R version 2.9.2 on Windows XP OS with RODBC version Version: 1.3-0. Has anyone out there in the R user community successfully appended records to a DB2 table on a remote database using the sqlSave function in the RODBC package? (or by any other means from R?) I posed a similar question a few months ago and unfortunately, did not receive a response. I was hoping recent upgrades to

Hi all and thanks in advance. I am regressing Time and Weight, and then predicting Weight at different Time. The format of the Time data is day/month/year. How can I get R to use time series data such as this? Keith -- M. Keith Cox, Ph.D. Alaska NOAA Fisheries, National Marine Fisheries Service Auke Bay Laboratories 17109 Pt. Lena Loop Rd. Juneau, AK 99801 Keith.Cox at noaa.gov marlinkcox at

Hi, member-list i'm searching for packages or routines to analise fisheries acoustics data. Something make import data of scientific ecosound (EK500) of paralel port, filter signal, analyse pattern of pixel fish shoals... Is there someone work with this ? If no, why dont make colective effort to do this?? (package fish_acoustic). Thanks Marcelo

Hi. I just installed Centos 4. I''m pretty sure that I chose to have it install postgresql but when the system came up, it wasn''t there. No worries. I installed it from the net with ''yum''. Unfortunately, when I started it up and it tried to init the database, I got a bunch of SELinux errors: Mar 3 13:24:22 dirty kernel: audit(1109874262.006:0): avc:

Hi, am very new to R and I've written an optim function, but can't get it to work least.squares.fitter<-function(start.params,gr,low.constraints,high.constraints,model.one.stepper,data,scale,ploton=F) { result<-optim(par=start.params,method=c('Nelder-Mead'),fn=least.squares.fit,lower=low.constraints,upper=high.constraints,data=data,scale=scale,ploton=ploton)

Dear R user: I am studying the allele data of two populations. the following is the data: a1 a2 a3 a4 a5 a6 a7 a8 a9 a10 a11 a12 a13 a14 a15 a16 a17 pop1 0.0217 0.0000 0.0109 0.0435 0.0435 0.0000 0.0109 0.0543 0.1739 0.0761 0.1413 0.1522 0.1087 0.0870 0.0435 0.0217 0.0109 pop2 0.0213 0.0213 0.0000 0.0000 0.0000 0.0426 0.1702 0.2128 0.1596 0.1809 0.0957 0.0745 0.0106