similar to: Formulating anova for partially nested model

Hi, I'm using R to plot some data and associated linear models, and I was wondering if there was a parameter to fix the line width in the graphs (for axis, lines, ...) because I explored the associated help pages and could not find it. Thx in advance, -- Ir. Yves Brostaux - Statistics and Computer Science Dpt. Gembloux Agricultural University 8, avenue de la Facult? B-5030 Gembloux

Dear members, I'm not very used with Unix systems, but I got an old Sun/UltraSparc workstation running Solaris 2.5.1 and I tried to install R on it to able able to do some batch R processing while working on my pc. I downloaded and installed following required and recommanded programs before installing R : perl 5.8.0, readline 4.3, gzip 1.3.5, zlib 1.1.4, make 3.80, jpeg 6b, libpng 1.2.4

Hello, I'm trying to plot a set of 3 dependant variables (y) against 4 predictors (x) in a matrix-like plot, sharing x- an y-axis for all the plot on the same column/line : y1/x1 y1/x2 y1/x3 y1/x4 y2/x1 y2/x2 y2/x3 y2/x4 y3/x1 y3/x2 y3/x3 y3/x4 In fact, this plot is a rectangular selection of the result of pairs(), limited to the relations between x's and y's

Hello everybody, I'm testing the randomForest package in order to do some simulations and I get some trouble with the prediction of new values. The random forest computation is fine but each time I try to predict values with the newly created object, I get an error message. I thought I was because NA values in the dataframe, but I cleaned them and still got the same error. What am I

Hello ! I'm struggling again against lattice graprhics. ;) I'm trying to produce a conditionnal xyplot with two conditionning factors (let's say A and B). I want the levels of those factors (A1, A2, etc) to show in the margins of the lattice plot, not in the strips between the panels. A1 A2 A3 plot11 plot12 plot13 B1 plot21 plot22 plot23 B2 I managed to remove the

Dear list, Here's a little problem I already solved with my own coding style, but I feel there is a more efficient and cleaner way to write it, but had no success finding the "clever" solution. I want to produce a factor from a subset of the combination of two vectors. I have the vectors a et b in a data-frame : > df <- expand.grid(a=c(0, 5, 10, 25, 50), b=c(0, 25, 50,

Dear list, I'm stuck with a little graphical problem. I'm generating several lattice plots which are printed in a single postcript device opened by > trellis.device(postscript, theme=canonical.theme("postscript", color=F), file="an_phase2_graph.ps", paper="a4", pointsize = 10, onefile=TRUE, horizontal=TRUE) Everything works fine,but some of these

Hello again, First I want to thank all the people who answered my question about line width in graphs. I promise I will learn the 'par' help page by heart for the end of the month ! I now want to trace some ellipses to emphasize groups of data. I found how to trace circles with 'symbols()', but no ellipse. I'm planning on writing my own function based on

Hi, a questionaire data-set with more than one row for one Person make me problems. (1) i want create a function which count the used row for the relevant id ! example: id c.row 1 1 1 2 1 3 2 1 2 2 4 1 4 2 4 3 4 4 4 5 (2) i have got a problem to filter & analyse , because the filter criteria is unlucky a row above or under the ratings which i would like to

What is the recommended/optimal way to perform aggregates on data frames with chron objects? Here is an example: >raw.data 1 07/09/01 4000 2 07/09/01 2000 3 07/11/01 1000 4 07/13/01 800 5 07/13/01 700 6 07/16/01 600 7 07/17/01 500 I'm trying to construct a function that would first aggregate the data (second column) by day (grouping by the first column) according to a

Hello I have a vector gd of 279 values of diameter of limpets. Using 'hist(gd, breaks=seq(10,60,5),plot=F)' I can obtain a vector of 10 centers of size-class ($mids). But I'd like to work on a vector of 279 values being the mids (centers of size-classes) corresponding to the respective 279 values in gd. I looked in the different manuals and in archives of list but didn't find

So in fact : I have a data wich colnames are : "id","VEMS5A","VEMS6A","VEMS7A","VEMS8A","VEMS9A","VEMS10A","VEMS11A","VEMS12A","VEMS13A","VEMS14A","VEMS15A","VEMS16A","VEMS17A","VEMS18A","cohort") With "VEMS5A" is the value when

Hello I want to compare the slope (let's say 'b') of a linear model obtained with lm() to a theoretical value (let's say 'th'). To do so, I think I should compute a 't value' using something like : (b - 'th')/standard.deviation(b) and then look at the p-value of this computed t. I don't understand how to do this in a simple way, just using lm()

Excuse me. I don't understant how to use the function "stack". I'm not sure to know the good syntax. I've got a data witch column names are the time when the mesures was done and an other who indicates the age of the subjetcs. I would like to do a new matrix with patients older than a value and the new matrix will contain in the first column the time of the mesures and in

this is a clearer (I hope) version of an earlier post - My problem is formulating the random = argument to give estimates of all 9 random components for this kind of setup where there are (I think) 9 variance/covariance components ... Study.1 Study.2 ... Study.5 Treatment T1: subject: 1 2 3 4 5 6 ... 13 14 15 Treatment T2: subject: 16 17

Hi, I have this table (BTW, published by FBI and representing age distribution of the WTC tragedy victims): "LABEL" "FREQ" "MIDPOINT" "1" "Infant (under 1)" 0 "2" "1 to 4" 5 "3" "5 to 8" 1 "4" "9 to 12" 3 "5" "13 to 16" 0 "6" "17 to

On Thu, Dec 29, 2016 at 05:44:50PM +0100, Giulio Teslano wrote: > a. What other types of extended wild card(s) options are there ? > > or is this still currently limited to these two characters '*?' ? As I said, the branch "adds support for arbitrary glob-style wildcard patterns (where * matches 0 or more characters and ? a single

Dear list, I'm migrating a project from Matlab to R, and I'm facing a relatively complicated problem for nls. My objective function is below: >> objFun <- function(yEx,xEx,tEx,gamma,theta,kappa){ yTh <- pdfDY(xEx,tEx,gamma,theta,kappa) sum(log(yEx/yTh)^2) } The equation is yTh=P(xEx,tEx) + noise. I collect my data in: >> data <-

I made an ANOVA with repeated mesures (aov(Mesure~Distance*Genre*Correct+Error(Sujet/(Distance*Genre*Correct)), data)) and I would like to extract the pvalues. The output is: ----------------------------------------------------------- Error: Sujet Df Sum Sq Mean Sq F value Pr(>F) Residuals 21 97.082 4.623 Error: Sujet:Distance Df Sum Sq Mean Sq F value Pr(>F) Distance

Hello, We have Xapian-Omega installed (Linux) and working in default mode and have browsed several interesting pages on the main site, at trac.xapian.org (the wiki) and in the mailing list. Having tested various search options (up to now only for Omega) we would like to ask a couple of questions. 1. Is it possible to search for Patterns of Text with Omega and/or Xapian Queries ?