Displaying 20 results from an estimated 10000 matches similar to: "foo.RData or foo.r?"
2003 Jul 28
3
data manipulation: getting mean value every 5 rows
Dear All,
I would like to ask you how to accomplish a little tricky data
manipulation. I have a large dataset, looking something like:
temp line cage number
18 18 1 6678.63
18 18 1 7774.458
18 18 1 7845.902
18 18 1 9483.578
18 18 1 8983.555
18 18 1 9181.052
18 18 1 9458.696
18 18 1 8138.616
18 18 1 7981.994
18 18 1 7556.491
18 18 1 7672.137
18 18 1 6607.776
18 18 1 8383.65
18 18 1 7129.852
18 18
2003 Aug 05
2
bwplot colours
Dear All,
is any way I can change the colours of the box and umbrella of a bwplot
without having to go to the length of:
box.rectangle<-trellis.par.get("box.rectangle")
box.rectangle$col<-"black"
trellis.par.set("box.rectangle", box.rectangle)
etc...
but straight from the call:
bwplot(y ~ x | z, mydata) ?
Regards,
Federico Calboli
2003 Jun 17
1
books
Dear All,
using R is noticed that I would benefit by going back to the books and
"refresh" my linear algebra and calculus.
Has anybody suggestions on which books are better as reference on the two
these two topics?
I was thinking to get Lang's "A First Course in Calculus", but I I have no
idea what could be a good one for linear algebra. Something detailed and
clear
2003 Jun 26
3
degrees of freedom in a LME model
Dear All,
I am analysing some data for a colleague (not my data, gotta be published
so I cannot divulge).
My response variable is the number of matings observed per day for some
fruitlies.
My factors are:
Day: the observations were taken on 9 days
Regime: 3 selection regimes
Line: 3 replicates per selection regime.
I have 81 observations in total
The lines are coded A to I, so I do not need
2003 Jul 14
1
methods help and glmmPQL
Dear All,
I would like to ask you to help me with my memeory. I remember using some
function that would list all the possible methods I could apply to an
object. Say, if I had an object of
class=lme,
it would tell me that that I could do stuff like
qqnorm(myobjct), or VarCorr(myobject). In general, a very complete list.
I though this list of all possible methods would pop out by typing
2003 Sep 23
3
problems installing Design and Hmisc libs
Dear All,
when I try to:
install.packages("Design"); install.packages("Hmisc")
I get the following error messages:
* Installing *source* package 'Design' ...
** libs
g77 -mieee-fp -O2 -fomit-frame-pointer -pipe -march=i586
-mcpu=pentiumpro -O2 -fomit-frame-pointer -pipe -march=i586
-mcpu=pentiumpro -c lrmfit.f -o lrmfit.o
make: g77: Command not found
make: ***
2005 Mar 08
5
removing message: [Previously saved workspace restored]
Dear All,
I saved by mistake the environment I was working in after typing q(),
and now I get the annoying message:
[Previously saved workspace restored]
I have already deleted all the objects in the environment, saving it as
an empty environment, so it's just a matter of nitpicking I suppose. The
message does not appear if I start R from any other place in the
directory tree.
I am
2003 Jul 17
2
line colors in lattice.xyplot with png device.
Hi,
R is very new for me, so excuse if my questions are too basic...
BTW, are there any forum where new R users could get help without
annoying this huge mailing list ?
In following code, I'd like to choose the color for each of the curve
diplayed.
png(filename = filename, width = 950, height = 600, pointsize = 10, bg =
"white")
xyplot(HITS+MS1*3+FREQ~TIME
2003 Jul 20
1
Problem with packages...
Hi, I just installed the last version of R taken from the CRAN mirron in
University of Vi?osa, Brasil. Since I use Mandrake 9.1 I get the proper .rpm
file. I installed fine and all the test runs OK.
I tryied to install sevela packages ussing the
install.packages("package.name")
and the installation runs OK... apparently.
However, when I try to use the new package I get the following
2003 Jul 27
2
continuous independent variable in lme
Dear All,
I am writing to ask a clarification on what R, and in particular lme, is
doing.
I have an experiment where fly wing area was measured in 4 selection lines,
measured at 18 and 25 degrees. I am using a lme model because I have three
replicated per line (coded 1:12 so I need not use getGroups to creat an
orederd factor).
The lines are called: "18"; "25";
2005 Jul 01
2
loop over large dataset
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
2005 Jul 21
3
vectorising ifelse()
Hi All,
is there any chance of vectorising the two ifelse() statements in the
following code:
for(i in gp){
new[i,1] = ifelse(srow[i]>0, new[srow[i],zippo[i]], sample(1:100, 1,
prob =Y1, rep = T))
new[i,2] = ifelse(drow[i]>0, new[drow[i]>0,zappo[i]], sample(1:100,
1, prob =Y1, rep = T))
}
Where I am forced to check if the value of drow and srow are >0 for each
line... in
2006 Sep 26
5
putting stuff into bins...
Hi All,
I have a vector of data, a vector of bin breakpoints and I want to put my data
in the bins and then extract fanciful informations like the mean value of each bin.
I know I can write my own function, but I would have thought that R should have
somewhere a function that took as arguments something like (data, breaks, what
to do with the data in the bins). I surey could not find it
2003 Jun 19
1
GLME
Hi All,
does anyone know if the package GLME by J. Pinheiro is available anywhere
in any form? checking on the archive I got that it was at some point, as as
a beta version (for S-Plus only, alas)...
Cheers,
Federico
=========================
Federico C.F. Calboli
Department of Biology
University College London
Room 327
Darwin Building
Gower Street
London
WClE 6BT
Tel: (+44) 020 7679 4395
2006 Feb 28
4
subsetting a list of matrices
Hi All,
I have a list of matrices:
> x
[,1] [,2]
[1,] 1 4
[2,] 2 5
[3,] 3 6
> y
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] 18 21 24 27 30 33
[2,] 19 22 25 28 31 34
[3,] 20 23 26 29 32 35
> z =list(x,y)
I want to create a second list that is has a subset each matrix in the
list subsetting so I get the 2nd and 3rd row of each (and
2007 Jun 26
2
fisher information matrix
Hi All,
a colleague wants to calculate the Fisher information matrix for a model he
wrote (not in R). He can easily get the neg-log-likelihood and the best fit
parameters at the minimum. He can also get negLLs for other parameter values too.
Given these data, is there a way in R to calculate the Fisher information matrix?
Best,
Federico
--
Federico C. F. Calboli
Department of Epidemiology
2003 Jan 16
1
graphics
Dear R community,
I need to plot the results of some simulations I did using QTL
Cartographer. I am plotting LOD scores over three chromosomes. The three
plot have to be one next to the other.
The procedure I am using is:
par(mfrow=c(1,3))
plot(x$x, x$y, ylim=c(0,35), type="l", col="blue", las=1, xaxs="i",
yaxs="i", xlab="X Chromosome",
2005 May 24
2
R unable to run on Mac OS 10.4 Tiger
Hello,
I'm running a PB G4 with Mac OS 10.4.1. I have downloaded the latest
version R-2.1.0a.dmg. It appears that R does not work. It launches
itself, but the window never gets ready, there is written "Loading R..."
and a small progress wheel keeps turning indefinitely.
Could someone help or suggest something?
THANKS !!
Guillaume
2006 Feb 08
2
logical condition in vector operation
HI All,
I have a data frame such as:
> test
x y p d
[1,] 1 0 10 21 0
[2,] 2 3 11 12 0
[3,] 3 4 12 23 0
[4,] 3 5 13 24 0
and I want to perfor some operations on the first two coulums,
conditional on the uneqaulity values on the 3rd and 4th columns.
For instance:
j = 3
test[test[,1] == j, 5] = test[test[,1] == j,2] + test[test[,2] == j,1]
gives me the result:
test:
x y p d
2006 Oct 26
2
pairs matchning
Hi
You could try to find an equivalent representation as a string and try to
match those.
> (A <- cbind(sample(1:2, 10, rep=TRUE), sample(1:2, 10, rep=TRUE)))
[,1] [,2]
[1,] 1 2
[2,] 1 2
[3,] 1 2
[4,] 2 2
[5,] 1 1
[6,] 1 2
[7,] 1 2
[8,] 1 1
[9,] 1 2
[10,] 1 1
> (B <- unique(A))
[,1] [,2]
[1,] 1 2