similar to: using nls to fit a curve to data

Hi, I'm programming in R and below is a summary of a generalized linear model: ************************************************** *** Call: glm(formula = offspring ~ degdays, family = quasi(link = "log", variance = "mu"), data = fecundity) Deviance Residuals: Min 1Q Median 3Q Max -0.76674 -0.29117 -0.09664 0.15668 1.00800 Coefficients: Estimate Std. Error t value

If I have a model line = lm(y~x1) and I want to use a for loop to change the number of explanatory variables, how would I do this? So for example I want to store the model objects in a list. model1 = lm(y~x1) model2 = lm(y~x1+x2) model3 = lm(y~x1+x2+x3) model4 = lm(y~x1+x2+x3+x4) model5 = lm(y~x1+x2+x3+x4+x5)... model10. model_function = function(x){ for(i in 1:x) { } If x =1, then the list

Hi I'm trying to fit arima models with the arima() function and I have two questions. ###### ##1. ## ###### I have n observations for my time series. Now, no matter what arima(p,d,q)- model I fit, I always get n residuals. How is that possible? For example: If I try this out myself on an AR(1) and calculate the fitted values from the estimated coefficients I can calculate n-1 residuals.

Hi all, I am reading a huge data set(12M rows) that contains family information, Offspring, Parent1 and Parent2 Parent1 and parent2 should be in the first column as an offspring before their offspring information. Their parent information (parent1 and parent2) should be set to zero, if unknown. Also the first column should be unique. Here is my sample data set and desired output. fam

Sorry that I was unclear. For an individual to qualify for my analysis I want both of the following two criteria to be fulfilled: First, I want to select measurement taken at a certain age: for the focal individual the year of measurement (year) should be the same as year.hatch Second, I want the focal individual to be born by a mother that reproduces for the first time. So the /parents /of

Hi people, I'd like to understand valid users and file permissions better. I have a share which is not behaving as I expect. [family] path = /home/shares/family create mask = 0664 directory mask = 0775 force group = parental guest ok = No valid users = @parental, @family writeable = Yes in Group parental are mum & dad; in group family are mum, dad and offspring. With file permissions

Dear Rlers, in the following dataset I would like to insert a new column that refers to the data of the previous row. My question is whether the probability of a female (Id) changes if she had given birth to a pup in the previous year. So my dataframe consists of the column Id, year (2003-2007 for each Id) and offspring (=of; 1-0): Id year of 1 2003 1 1 2004 0 1

All, Sometimes when I do an rsync, it times out and leaves on orphaned process on the remote machine. I am running rsync over ssh using public key authentication. I synchronize about twenty-five machines, and this problem only appears on one of them. It worked until recently on this machine. I don't know what changed or what is different. I run rsync with a command line like this

Hi, All I have to do some job with PV-driver, who could give me some advice, thank you very much! 2010-10-19 xudongxen _______________________________________________ Xen-users mailing list Xen-users@lists.xensource.com http://lists.xensource.com/xen-users

http://bugzilla.mindrot.org/show_bug.cgi?id=797 Summary: ignores setting of CONSOLE var in /etc/default/login Product: Portable OpenSSH Version: 3.7.1p2 Platform: UltraSparc OS/Version: Solaris Status: NEW Severity: normal Priority: P2 Component: sshd AssignedTo: openssh-bugs at mindrot.org

Hi, My Samba server (Solaris 7, 64-bit mode) with AFS option has been running for a year. I just changed the "Max Log Size" global parameter few times today. Each time I committed the change, I got the following same set of error messages in my Samba log (the log.smb file): [2002/01/03 17:40:56, 0] lib/util_sock.c:set_socket_options(151) Failed to set socket option

I am getting "Chi-squared approximation may be incorrect in: chisq.test(x)" with the data bleow. Frequency distribution of number of male offspring in families of size 5. Number of Male Offspring N 0 518 1 2245 2 4621 3 4753 4 2476 5

Dear R users, I encounter a problem regarding number of significant digits on y-axis. Below is my basic code: myplotkid<-xyplot(expected_offspringnumber~afr|decade,groups=SES,data1, auto.key=list(space="right"),layout=c(9,1),xlab="",ylab="Offspring number", aspect="fill",scales=list(x=list(draw=F)),strip=T) >From this code, you can see there are 9

Is there a way I can have R automatically execute the commands in a source file without ever having to use R interactively? If so, what arguments should I pass to the UNIX call to R? I need to do this to run several R jobs in parallel. An alternative may be to have R and an R package behave as a stand-alone application that can be called from the UNIX command line. Is there any documentation on

R-Users, Hopefully someone can shed some light on these questions as I had little luck searching the archives (although I probably missed something in my search due to the search phrase). I estimated multinomial probabilities for some count data (number successful offspring) ranging from 0 to 8 (9 possible response categories). I constructed a barplot (using barplot2) and I want to

Hi all, I wonder if anyone can help, I have a dataframe with columns for... 'I.D' 'age' 'mothers I.D' 01 5 03 02 6 06 03 16 NA 04 8 06 05 3 NA 06 17 NA I need to create a new column for 'mothers age' which puts the age of the individual with 'mothers i.d' into the row for her offspring (so

Hi, I am sampling two random columns from females and two random columns from males to produce tetraploid offspring. For every female I am sampling a random male. In the end I want to write out a a matrix with all the offspring, but that does not work. I get always only the offspring from the last females. There must be a mistake in my script: moms<-read.delim("females.txt",

Hi. Thank you very much in advance for your help. I have generated data from two simple linear models and used k-means clustering (stats4) to identify two clusters in the generated data. Next, I would like to do simple linear regression for each separate cluster. I can do this if I first use the cluster labels to define two separate data frames with the subset function. However, I would

Dear R users, I encounter a problem regarding space control in xyplot. Basically, I want to control spaces between label, tick and axis. I remember there is a function called mgp in general plot. Is there a similar function for xyplot or xYplot? Below is my basic code: myplotkid<-xyplot(expected_offspringnumber~afr|decade,groups=SES,data1,

Hi All, I'd like to ask for a few clarifications. I am doing some calculations over some biggish datasets. One has ~ 23000 rows, and 6 columns, the other has ~620000 rows and 6 columns. I am using these datasets to perform a simulation of of haplotype coalescence over a pedigree (the datestes themselves are pedigree information). I created a new dataset (same number of rows as the pedigree