similar to: Running glm in batch and exporting results (AIC) to HTML

Dear all, I get an error message when I run my model and I am not sure what to do about it. I try to determine what factors influence the survival of voles. I use a mixed-model because I have several voles per site (varying from 2 to 19 voles). Here is the model: ### fm5 <-lmer(data=cdrgsaou2, alive~factor(pacut)+factor(agecamp)+factor(sex)+ResCondCorp+(1|factor(cdrgsa ou2$ids)),

Hello alls, I found in the literature a technique that has been evaluated as one of the more robust to assess statistically the significance of the loadings in a PCA: bootstrapping the eigenvector (Jackson, Ecology 1993, 74: 2204-2214; Peres-Neto and al. 2003. Ecology 84:2347-2363). However, I'm not able to transform by myself the following steps into a R program, yet? Can someone could help

Dear all, I'm trying to combine in R the two following datasets. I'm would be very grateful if someone could help me with this problem. Dataset 1: X Y Z 1 27 48 2 25 50 3 27 40 4 28 56 Where X is a unique number and Y and Z some variables Dataset 2: X I J 1 10 12 1 12 12 1 23 30 1 40 46 2 7 8 2 8 9 2 9 8 3 2 0 4 98 87 4 78 89 Where X IS NOT this time a unique number. I and J are some

Dear alls, Could someone tell me how to select a subset of string observations (e.g. "females" in a sex column) with sqlQuery in the RODBC library? Indeed, I'm trying to select a subset of observations on my access database with: female<-sqlQuery(mychannel,"SELECT Micromammiferes.sex FROM Micromammiferes WHERE (((Micromammiferes.sex)="females"));") The sql

Dear all, I ran the following model without any problem previously to the update of lme4: fm2<-lmer(data=NGud,family=poisson, seed~hab*seedtray +(1|site)+(1|site:hab)) I have 25 sites, 2 habitats ("hab") per site, 8 seedtrays per habitat (4 as control, and 4 as treatment), and I'm interested in comparing the number of seed in seedtrays as a function of the treatment and its

Dear Rusers, Trying to reduce my for loops addiction, could somebody tell me if there are ways to simplify (and perhaps accelerate ?) the following line for (i in 1:N) for (j in 1:N) m[i,j] = b[i]-b[j]; (where m is a NxN matrix and b a vector of length N) Thanks for any hint.

Hello Listers, I am trying to sample a vector to create a new one of sample length, witha sum equal to the sum of the initial vector: initial = 10, 30, 10 (sum=50) sample example = 5, 35, 10 (sum=50) or 25, 15, 10 (sum=50), etc ... My problem is to control the sum, so it stays constant. Any suggestions would be very helpful ! Thank you in advance, Eric

Dear all, I'm building binomial mixed-model using lme4 package. I'm able to obtain outputs properly except when I include two particular variables: date (from 23 to 34; 1 being to first sampling day) and Latitude (UTM/100 000, from 55.42 to 56.53). No "NA" is any of those variables. In those cases, I get the warning message: "nlminb failed to converge" I tried to

hi, I have a problem. The R shows that Error in 1/ue : non-numeric argument to binary operator. Here is the code: # simulation of tempered stable processes by compound Poisson approximation tsp<-function(n,e,a,c,lema){ x<-numeric(n) for (i in 1:n){ repeat{ w<-runif(1) v<-runif(1) x<-e*(w^(-1/a)) fex<-function(x1){ if (x1>=e) {a1<-a*(e^(-a))/(x^a-1)}}

Hello, It seems I need to restart webrick everytime I modify a module in order for my modifications to be taken into account. Is this normal and is there a workaround? Regards, -- ---------------------------------------------------------------------- Yannick Majoros http://www.inma.ucl.ac.be/~majoros Informaticien UCL/INMA-MEMA 4, avenue G. Lema?tre B-1348 Louvain-la-Neuve Tel:

Hello, How can I pass a hash as parameter from a template with url_for? e.g. : if I have test = { :a => "a", :b => "b" } How can I pass that with my request? Is it possible to pass nested hash''es to? Right now, url_for :action => :test, test doesn''t seem to work. Thank you, --

Hello all, I have a vector of character strings, in which I have letters, numbers, and symbols. What I wish to do is obtain a vector of the same length with just the numbers. A quick example - extract of the original vector : "lema, rb 2%" "rb 2%" "rb 3%" "rb 4%" "rb 3%" "rb 2%,mineuse" "rb" "rb" "rb 12"

Hello, Rather new to RoR, so I''m not sure about the terminology and such :-/ . My question is: how could I know the model which is associated to a controller? For example: - controller class is TestController, which is a subclass of ApplicationController - associated model class is Test How could I write some code in a method of ApplicationController to dermine the current

Hello, I always have the same kind of design problem, which I am not sure how to solve correctly: - say I have a variable, "status", which I want to calculate once in my controller, and access a couple of times in my view - "status" variable initialization could happen in a subclass, overriding the regular method for that What I do now is using an instance variable

Hello, Here is the description of a problem we habe been having for some time. Originally posted on the postfix-users mailing list, but no solution found to date. >>>> We have been having problems with our mail system for a couple of >>>> months. Every now and then, some mailboxes get corrupted: a >>>> half-truncated message, usually without headers, is at

One more question: Is there some way to set the human_name of a column? e.g.: human_name for column address1 shouldn''t be Address1 but "Address, line 1". If not, should I make a hash with my custom names? Best regards, -- ---------------------------------------------------------------------- Yannick Majoros http://www.inma.ucl.ac.be/~majoros Informaticien UCL/INMA-MEMA

Hello, I have a People table and a Addresses table. A person can have one or more addresses, but should at least have one, so there is a address_id field in People. Now, I would like to have a form to fill the name of a new person and its address from the same place. I could use person.address.country, for example (it works), but I would like to simply use person.country for some reasons

Dear Rusers, I would like some comments about the following results (under R-2.2.0) > m = matrix(1:6 , 2 , 3) > m [,1] [,2] [,3] [1,] 1 3 5 [2,] 2 4 6 > z1 = m[(m[,1]==2),] > z1 [1] 2 4 6 > is.matrix(z1) [1] FALSE > z2 = m[(m[,1]==0),] > z2 [,1] [,2] [,3] > is.matrix(z2) [1] TRUE Considered together, I'm a bit surprised about

Basically, i want to create a boxplot from population data, with a solid line for the median, and then a dotted line for a specific persons value. Is there a way to artificially introduce another line or mark of some kind in a box plot? Thanks again for all of yoru help. Hopefully soon I wont have to keep bugging you all. ~Erithid :-p

Dear R users, First of all sorry for this question, surely quite naive. (I searched on the R site but was unable to find by myself). I have a table, called infile : 1 2 3 4 5 6 7 8 9 I would like to read it and ignore the rows with 1st element > 3 I do it now with a for loop, and it's ok, but I was expecting something simpler, like : intable = read.table(infile); newtable =