similar to: Help w/ matrix calc

*Thanks* all those who took the time to help me (even if the "question" was not related to - the use of - R). Now I think I can soundly make my point w/ the referee (can I use your replies? If so I intend to properly cite its use?!?). Regards, Eduardo Esteves ps - Sorry for not explaining the "biological details" of my posting: RNA/DNA is the ratio of RNA content to

Dear All, I'm confortable with xyplot(...) and panel.lmline(...) statements (at least I thought I did :). I've used the following code to plot the decline in log-abundance of fish larvae (no.larvae) with age (age.cls, 4 to 27 days-old) for specific dates of sampling (day, 9 dates). I further plotted data with different colors and regression lines for ages 5-14 d and 17-23 d in a 7-by-1

Dear All, I'me having (much) trouble understanding why it happened and answering a referee's comment to part of a submitted manuscript. I've tried to google for help but... I'm really confident that although this is a R-Help list someone can help me! I used R to do an ANCOVA w/ RNA/DNA as the dep var, sl as the indep var and gut (a factor w/ levels: prey and empty) as the

Dear All, A few weeks ago I posted a question to this list but unfortunately got no answer! A friend warned me of my english. Again, the problem is: A 2-level 5-factors completely randomized design was used to investigate the potential effects of those factors on a solution's characters (several response variables). Each response-variable was measured repeatedly (13 times) during a

Dear All, I'm a Marine Biologist using R (in a most informal and applied way) for a couple of years. Actually, I've been using R to analyse the results supporting my thesis. I would like to put together (in the same file) the data (.txt files) and the scripts files (.R files) for each chapter. How can I do this? Thanks in advance, Eduardo Esteves [[alternative HTML version deleted]]

Dear All, I'd like to fit a "kind" of logistic model to small data-set using nonlinear least-squares regression. A transcript of R-script are reproduced below. Estimated B and T (the model's coeff, herein B=-8,50 and T=5,46) seem appropriate (at least visually) but are quite diff from those obtained w/ SPSS (Levenberg-Marquardt): B=-19,56 and T=2,37. Am I doing something wrong in

Dear All, I'm trying to calculate the following: g.rate<-(SL-int)/NO where SL and NO are individual length and counts for the same subject, and the int value is a day-specific value i.e. for each DAY=88, 101..., 172 I'd like to use a different values of int=9.32, 8.43, ..., 9.81. All variables are compiled in a larger data.frame of the form: CODE SL NO DAY (...) 123 12.5 14

Dear All, I'm having trouble tweaking a forest plot made using the R meta-analysis package metafor. I did the analysis based upon the correlation coeff from studies and plotted the corresponding forest plot easily > q2<-rma(yi,vi,mods=cbind(grupo),data=qim) > q2 > forest (q2,transf=transf.ztor,digits=3, ... ,alim=c(0,1),refline=.5) > text(-1.55,42,"Esp?cie

Hi, I have a dataset like this, > testdata Grouped Data: expr ~ visit | subject expr visit subject 1 6.502782 V1 A 2 6.354506 V1 B 3 6.349184 V1 C 4 6.386301 V2 A 5 6.376405 V2 B 6 6.758640 V2 C 7 6.414142 V3 A 8 6.354521 V3 B 9 6.396636 V3 C I tried the command >

Dear all I have made some lattice graphs (bwplot(dead ~ treat | group,...). 'group' is a factor with three levels (Artemia larvae, Abalone larvae and Abalone spat), and the result is a graph with three panels corresponding to the levels in 'group', ordered in alphabetical order from bottom to top, as expected. How does one re-order the order in which the levels of

Dear all, How do I make R show the R-squared (deviance explained by the model) in a logistic regression? Below is how I write my syntax. Basically I want to investigate density-dependence in parasitism of larvae. Note that in the end I perform a F-test because the dispersion factor (residual deviance / residual df) is significantly higher than 1. But how do I make R show the

> Background: > OS: Linux Ubuntu Dapper > release: R 2.3.0 > editor: GNU Emacs 21.4.1 > front-end: ESS 5.2.3 > ----------------------------- > Colleagues > I have two spatial datasets (latitude, longitude, fish eggs) and (latitude, longitude, fish larvae) at the same 280 stations (i.e. 280 cases). I want to determine if the 2 datasets are spatially correlated. In other

Background: OS: Linux Ubuntu Dapper R : Version 2.3.1 (2006-06-01) Emacs 21.4.1 ESS Colleagues This is neither a bug report, nor a complaint, so I think this is the appropriate list. I am getting an odd font problem with a pdf generated using pdf(file = "../figures/PCNM.models.larvae.fluor.zoo.stability.2001.and.2002.pdf", title="", family="serif",

Background: OS: Linux Mandrake 10.1 release: R 2.0.0 editor: GNU Emacs 21.3.2 front-end: ESS 5.2.3 --------------------------------- Colleagues I have a set of lattice plots, and want to plot 4 of them on the page. I am having trouble with the layout. grid.newpage() pushViewport(viewport(layout = grid.layout(2,2))) pushviewport(viewport(layout.pos.col = 1, layout.pos.row = 1)) working trellis

I'm using the image.plot() function (fields package), but I want to enlarge the characters of the legend (as they are too small to be read in a combined figure), but there is no way I can find a command to do this. I can enlarge the legend bar (with legend.witdh), axis character size (cex.axis) or the total legend size (legend.shrink), but not the character size of the legend characters

Hi all I am not sure this question is better suitted for the R-devel or R-help, but here it goes anyway I am not a software developer but I have produced a personal set of functions that I would like now to combine in a package. I am finding this really hard to do only by following the manuals provided (the language is for software developers I guess which I am far to be). Even so I could

Full_Name: Fernando Henrique Ferraz Pereira da Rosa Version: 1.8.0 OS: Linux 2.4.21 Submission from: (NULL) (200.206.211.169) After installing R 1.8.0, the R DBI interface stopped working. I tracked it down as a problem in the 'methods' package, that comes in the default installation. Somehow the function '.valueClassTest' which is defined on package 'methods',

Hello, I have a matrix with values, with columns c1..cn. I need the values to be normalized between 0 and 1 by column. Therefor, the 0 should correspond to the minimum value in the column c1 and 1 should correspond to the maximum value in the column c1. The remaining columns should be organized in the same way. Does a function in R exists for this purpose? Thanks, Rui [[alternative HTML

Perhaps this is a common question but I haven't been able to find the answer. I have data with many factors, each taking many values. However, only relatively few combinations appear in the data, ie have nonzero counts, in other words the resulting table is sparse. Say we have 10 factors each with 10 levels. The result of table() would exceed the memory space (on a 32bit machine). Is there

Dear all. I am searching for KMeans ++ for R. I cannot find it. Do you know any package with it? Best regards, Rui