similar to: number of rescaling cycles in decorana

Is there a way to apply a function, say cor(), to each combination of some number of variables, and this, without using loops? For example, I have day, hour, var1 and var2. How could I compute cor(var1,var2) for each day*hour combination and obtain a matrix with day, hour and the cor value for each combination? Thanks for your time, Marc =================== Marc Bélisle Professeur adjoint

Howdee, My question appears at #6 below: 1. I want to model the growth of each of a large number of individuals using a 4-parameter logistic growth curve. 2. nlme does not converge with the random structure that I want to use. 3. nlsList does not converge for some individuals. 4. I decided to go around nlsList using: t(sapply(split(data, list(data$id)), function(subd){coef(nls(mass ~

Howdee, *** I know that the lmList() function exists, yet I don't want to use it. *** Would anyone be kind enough to tell how I can apply the function lm() to each level of a given factor so to obtain the intercept and slope for each factor level within a matrix? For instance, suppose a dataframe containing 3 variables: id, x and y. I want to compute the function lm() for each level

Hi there, I have 3 numeric variables: day (e.g., 05), month (e.g., 06), year (e.g., 2009). I would like to create a (string) variable of the following form: month/day/year (e.g., 06/05/2009). I would be grateful to anyone who could point me toward a solution. Sincerely, Marc =================== Marc B?lisle Professeur agr?g? Chaire de recherche du Canada en ?cologie spatiale et en ?cologie

Howdee, I am able to fit a 4-parameter logistic growth curve to a dataset which comprise many individuals (using R v. 2.3.1). Yet, if I want to obtain the parameters for each individual (i.e., for each 'id') using nlsList, then I obtain an Error message which I have trouble interpreting. Any advice as to how I can solve this problem? Thanks for your time, Marc > reg <-nls(mass ~

Howdee, One of my student spotted something I can't explain: a probability >1 vs a normal probability density function. > dnorm(x=1, mean=1, sd=0.4) [1] 0.9973557 > dnorm(x=1, mean=1, sd=0.39) [1] 1.022929 > dnorm(x=1, mean=1, sd=0.3) [1] 1.329808 > dnorm(x=1, mean=1, sd=0.1) [1] 3.989423 > dnorm(x=1, mean=1, sd=0.01) [1] 39.89423 > dnorm(x=1, mean=1, sd=0.001) [1]

Hi After trying a simple decorana analysis (from the vegan package) on a simple data frame which contains no NA's the following error was returned: > tt_decorana(covN) Error in decorana(covN) : NA/NaN/Inf in foreign function call (arg 1) Have any vegan users come across this error and know what can be done about it? Cheers, J

Hi I’m working in behavioural ecology I’m interesting in doing a multi-factorial MANOVA with a random effect as of independent variable. Do you know if it’s possible to do that and how I can do it? Thanks Julien Julien Martin Laboratoire de Denis Réale Chaire de Recherche du Canada en écologie comportementale Canadian Research Chair in behavioural ecology Département des sciences

In the general sense you may get some help by looking at the control dependence graph. - dibyendu ----- Original Message ----- From: Jianfei Hu [mailto:hujianfei258 at gmail.com] Sent: Saturday, August 25, 2012 07:43 AM To: LLVMdev at cs.uiuc.edu <LLVMdev at cs.uiuc.edu> Subject: [LLVMdev] How to Check whether BasicBlock resides in a conditional branch Hello All, I want to dertermine

dear R-users, I have to model counts where all counts above some threshold have been censored. In the same dataset I have too many zeroes for a Poisson or even a negative binomial distribution to make sense, so I would need a zero-inflated-censored negative binomial family for use in glmer (or glmmADMB?). That seems not to exist. my question is : how could I add a custom-built family of

Hi Jianfei Hu, the GVN pass does something like this in the logic around GVN::propagateEquality. If in your example it was if a == 2 // BasicBlock A then then it replaces all occurrences of a with 2 in BasicBlock A. For this it needs to understand which basic blocks can only be reached via this conditional edge "a == 2". Ciao, Duncan. > Hello All, > > I want to

Hello All, I want to dertermine whether a basicblock is in a conditional branch. such as, //============================= if a > 2 // BasicBlock A then BasicBlock B endif BasicBlock C //============================= What I want to identify is BasicBlock B, which is in a condtional block. but C is not. In other words, I want to distinguish BasicBlocks that * must * be executed and that

Can't you do it by performing some analysis on CFG? You can traverse that structure with BFS. And after that for all the BB you have visited more than once, you try to find a parent that has a branch instruction as a terminator. Additionally you ensure that there are no BB with branches as terminators on your way. If such parent exist, you mark that there is exist a direct connection between

2012/8/25 Iaroslav Markov <ymarkov at cs.stonybrook.edu>: > Can't you do it by performing some analysis on CFG? You can traverse that structure with BFS. And after that for all the BB you have visited more than once, you try to find a parent that has a branch instruction as a terminator. Additionally you ensure that there are no BB with branches as terminators on your way. If such

Dear all, Does someone know an R function implementing the method of Sison and Glaz (1995) (see full ref below) for computing Simultaneous confidence intervals for multinomial proportions? As alternative method, I think to boostrap the mean of each proportion and get in that way confidence interval of the mean. I observed 21 times a response that could be one out of 8 categories

Hello, with Rails 3.0.3 "Café Noir ".strip => "Café noir" but "Café ".strip => "Caf\303\251" In fact, strip() doesn''t works if the last printable character is accentuated. Surprisingly " écologie".strip works fine. I''ve tried to dig deeper in active_support multibyte source code but didn''t found any solution.

Hello, I apologize if this seems like an obvious question, but I have been looking everywhere and have yet to find an answer. I am doing a multinomial regression with multinom() in the nnet package. I have a 3 level ordered response (ordered()) variable and 4 predictors, 3 of which are numerical and one which is an ordered factor (also ordered()) with 5 levels (a, b, c, d, e). My question is in

I am using R version 2.12.1 on Win XP. I have done a heatmap with dendrogram using the heatmap.2 function. The heatmap basically looks like I want it to be, but the labels of the columns are cut off. I.e. the textual labels of the columns, although they are not very long (less than 12 characters), do not fit into the window and can not be read entirely. If I manually resize the graphics

When I export a graphic into a windows metafile, import that file into a Word document, the width of the character of my title isn’t uniform anymore. The only way I found for solving this problem is to use the courier font, but I would like to use a font like (font.main=1). Does anybody know about that problem? obs<-matrix(scan("D:\\.....) win.metafile(filename =

Hi, Let say I have this matrix: > mat<-matrix(cbind(rnorm(20),rnorm(20)), ncol = 2) And I want to rescale values of column [,1] and [,2] using values from row 1 to 5, such that the values of row 1:5 should be rescale to the same amplitude (kinda like take the z-score of population from row 1:5). column 1 and 2 are two different samples need to be compared. I hope I made myself clear