peter dalgaard gmail.com> writes:
>
>
> > On 25 Jul 2015, at 21:49 , Atte Tenkanen
utu.fi> wrote:
> >
> > Hi,
> >
> > I have tried to find a way to find opposite
or complementary colors in R.
> >
> > I would like to form a color circle with R
like this one: http://nobetty.net/dandls/
col
csiro.au> writes:
> I'm trying to fit a binomial GLM with user defined
link
function (negative exponential), however I seem to
> be unable to find the correct starting values to
initialise such a model. I've tried taking starting
> values from a logistic and log models fit to the
same data an
Marc Girondot yahoo.fr> writes:
> > outer(0:1, 0:1, FUN=function(x, y) {1})
> Erreur dans outer(0:1, 0:1, FUN = function(x, y) { :
>dims [produit 4] ne correspond pas à la longueur de l'objet [1]
Because whatever the dimensions of your 2 input vectors,
this function simply returns the value
Doran, Harold air.org> writes:
> I am trying to generate a binary matrix where
row
in the matrix is guaranteed to have at
least one 1.
> Ideally, I would like most rowSums to be equal
2 or 3
with some 1s and some 4s. But,
rowSums cannot be equal
> to 0.
>
> I can tinker with the vector o
Chen, George roswellpark.org> writes:
> My apologies for asking this question that may have
been asked before. I am trying to plot activity
> dependent on time conditioned by the subject.
Code for sample data below.
> So I have something like this
> xyplot(Activity~Time|Subject).
> This works fi
Roland Deutsch tuwien.ac.at> writes:
> in my research I frequently work with binomial
response models, which
> are of course part of the generalized linear
models. While I do use
> common link functions such as the logit, probit
and cloglog, I often
> have the need of invoking the lesser-kno
Alaios yahoo.com> writes:
> I have in my code some vectors that are not of equal size.
I would like to be able for each of these vectors
> select 6 elements that are (almost) equally spaced.
So the first one would be at (or close) to the beginning
> the last one at (or close) to the end and the
Martin Batholdy googlemail.com> writes:
> I would like to colour different areas of a plot.
> But I don't know how to do this efficiently.
>
> here an example:
> (t = time)
>
> t <- 1:100
> bg_colors <- c(rep('green',20), rep('yellow',10),
rep('green',20), rep('red',5),
> rep('yellow',45))
>
Ken Knoblauch inserm.fr> writes:
>
> Bos, Roger rothschild.com> writes:
> > I am using a sum of squared differences in the
> objective function of an optimization problem I am
> doing and I
> > have managed to speed it up using the
> outer function ve
Bos, Roger rothschild.com> writes:
> I am using a sum of squared differences in the
objective function of an optimization problem I am
doing and I
> have managed to speed it up using the
outer function versus the nested for loops, but my
suspicion is that
> the calculation could be done even q
Martin Batholdy googlemail.com> writes:
> I have to import multiple R-files.
> Each file consists of several functions with the same
function name across the R-files.
> When I import all files one by one (with source())
I overwrite the function definition of the previous file
> until only the ve
On Jun 12, 2013, at 4:36 PM, Ken Knoblauch wrote:
You seem to treating the input values as xyY when they should be XYZ
(case matters).
So, I would do something like this
D65 <- c(0.3127, 0.329, 0.3583)
X <- 100 * D65[1]
Y <- 100 * D65[2]
Z <- 100 * D65[3]
XYZ <- data.frame
You seem to treating the input values as xyY when they should be XYZ
(case matters).
So, I would do something like this
D65 <- c(0.3127, 0.329, 0.3583)
X <- 100 * D65[1]
Y <- 100 * D65[2]
Z <- 100 * D65[3]
XYZ <- data.frame(X = X, Y = Y, Z = Z)
convertColor(XYZ, from = "XYZ", to = "sRGB")
I'd have to look it up and I'm not home at the moment. Can see later on. I
would have thought that it would be normalized to have a jnd equal to 1 but I'm
not sure.
Ken
Sent from my iPhone
___
Ken Knoblauch
Inserm U846
Stem-Cell and Brain Research Institute
18 av du Doyen Lé
Hi John,
Out of curiosity and if it is not much trouble, I would be curious if Luv
worked any better than Lab. I think that Luv is supposed to be preferred for
monitors and Lab for surfaces but they are generally pretty similar.
Best,
Ken
Sent from my iPhone
___
Ken Knoblauch
Inserm U846
John Fox mcmaster.ca> writes:
> I'm interested in locating the named colour that's "closest" to an
arbitrary RGB colour. The best that I've
> been able to come up is the following, which uses
HSV colours for the comparison:
>
> r2c <- function(){
> hexnumerals <- 0:15
> names(hexnumeral
t Ken points
out (and I do appreciate him making these points). One can readily
demonstrate the gamut limitations by printing the diagram Ishida
links to on different devices. My hope is to get something close
and include a disclaimer. Bryan
On Mar 18, 2013, at 7:08 AM, Ken Knoblauch
ishi soichi gmail.com> writes:
> Has anyone plotted or is it possible to plot
>
> CIE *xy* chromaticity diagram
>
> http://en.wikipedia.org/wiki/File:CIE1931xy_blank.svg
>
> I need this plot in color.
>
> ishida
And following up on my previous mail (diatribe), after
having contemplated th
ishi soichi gmail.com> writes:
> Has anyone plotted or is it possible to plot
>
> CIE *xy* chromaticity diagram
>
> http://en.wikipedia.org/wiki/File:CIE1931xy_blank.svg
>
> I need this plot in color.
>
> ishida
>
I think that plotting the spectral locus and the line of purples is trivial,
Michael Friendly yorku.ca> writes:
>
> To partly answer my own question: It wasn't that hard to hack the
> result of model.matrix() to remove the intercept,
>
> remove.intercept <- function(x) {
> if (colnames(x)[1] == "(Intercept)") {
> x <- x[,-1]
> attr(x,
14.3 Hmisc_3.10-1 R2HTML_2.2 svMisc_0.9-65
TinnR_1.0-5tools_2.15.2
-- Bert
On Sun, Nov 18, 2012 at 8:09 AM, Ken Knoblauch
wrote:
Tom Roche pobox.com> writes:
As described @
<<< clipped >>>
However I will need to before-and-after compare this to
the r
lattice with the _value_ of
the level (an
> atmospheric pressure), rather than the name or
index of the level.
> How to do that?
>
> TIA, Tom Roche pobox.com>
maybe, see ?strip.custom in lattice
--
Ken Knoblauch
Inserm U846
Stem-cell and Brain Research Institute
Department of Integ
vice?
Try looking at
http://R.research.att.com/libs/
>
> TIA -- Christian
>
--
Ken Knoblauch
Inserm U846
Stem-cell and Brain Research Institute
Department of Integrative Neurosciences
18 avenue du Doyen Lépine
69500 Bron
France
tel: +33 (0)4 72 91 34 77
fax: +33 (0)4 72 91 34 61
portable: +33
Jessica da Silva gmail.com> writes:
> I am trying to create a scatterplot, coding each point to
one of 5
> populations. I was successful when I did this for one
set of data, yet
> when I try plotting other data a blank plot appears
(although the axes are
> labelled and I can fit the regression
Christof Kluß email.uni-kiel.de> writes:
> But now I search a self-starter for y = a + b * c^x
>
> Do you know a package that implement such self-starter or is there a
> easy way to do it myself?
You might be able to adapt the SS.calib function
from the psyphy package for this.
It was design
es several data sets from published psychophysical
experiments using detection and rating scale measures for
estimating signal detection parameters, psychometric functions,
classification images, etc.
Best,
Ken
--
Ken Knoblauch
Inserm U846
Stem-cell and Brain Research Institute
Department of I
rstand what you mean by "You
seem to be getting complete separation on X5 "?
Can you please be more elaborate?
Thanks,
On Thu, Apr 12, 2012 at 4:06 PM, ken knoblauch
wrote:
Christofer Bogaso gmail.com> writes:
Dear all, I am fitting a LOGIT model on this Data...
Christofer Bogaso gmail.com> writes:
> Dear all, I am fitting a LOGIT model on this Data...
<< snip >>---
> glm(Data[,1] ~ Data[,-1], binomial(link = logit))
>
> Call: glm(formula = Data[, 1] ~ Data[, -1], family = binomial(link = logit))
>
> Coefficients:
> (Intercept) Data[, -
ken knoblauch inserm.fr> writes:
>
> Michael Bach gmail.com> writes:
> > how do I e.g. square each second element of a
> vector with an even
> > number of elements? Or more generally to
> apply a function to every
> > 'nth' element of a vec
example:
> v <- c(1, 2, 3, 4)
> mysquare <- function (x) { return (x*x) }
> w <- applyfun(v, mysquare, 2)
> then w should be c(1, 4, 3, 16)
> Michael Bach
Hi Michael,
v^(2 - seq_along(v) %% 2)
[1] 1 4 3 16
Ken
--
Ken Knoblauch
Inserm U846
Stem-cell and Brain Research
t; Is there any way to tell glm() to add this
> parameter in the estimation or do I have to write my own estimator with
> optim()?
If the parameter cannot be made into a coefficient of the linear predictor,
then I'm afraid that you will have to roll your own.
> Thanks,
>
&g
How can I do then?
>
> What I know is only for 2 vectors via "intersect" function,
but don't know how to
deal with multiple
vectors.
>
Reduce(intersect, list(v1 = c("a","b","c","d"),
v2 = c("a","b&quo
Christof Kluß email.uni-kiel.de> writes:
> Am 02-01-2012 10:54, schrieb ken knoblauch:
> > Christof Kluß email.uni-kiel.de> writes:
> >> lme<- lme(conc ~ name/time - 1,
> >> random=conc~time|nr,method="ML",data=measurements)
> > see pl
me) and
> seperate colors for the measurements (nr).
>
> How would you do that?
>
> thx
> Christof
>
see plot.augPred in the nlme package
--
Ken Knoblauch
Inserm U846
Stem-cell and Brain Research Institute
Department of Integrative Neurosciences
18 avenue du Doyen Lépine
ibm.com/dx/proceedings/
pravda/truevis.htm
who was (is) quite concerned with this issue, as well,
as the excellent article by Zeileis, Hornik and Murrell
http://statmath.wu.ac.at/~zeileis/papers/
Zeileis+Hornik+Murrell-2009.pdf
HTH,
Ken
--
Ken Knoblauch
Inserm U846
Stem-cell and Brain Research
ot;http://r.789695.n4.nabble.com/file/n3870788/111004_Lode_Outlines.csv";,
header = TRUE,sep = ",",)
par(mfrow = c(1, 2), pty = "s")
plot(z ~ y, Data_poly, type = "l")
fh <- with(Data_poly, which(z > 240))
D_poly <- rbind(Data_poly[fh, ], Data_poly[-rev(fh), ])
D_
you describe, the
subjects are nested in this, i.e., some had a high effort of 5
and others of 3. Perhaps, the following would work then
glmer(y ~ EffortLevel/(effort + costs + scr) + (1 | id), family = binomial)
I think that if each observer has a unique id, that the nesting
will be automatic for
in the documentation of the R.matlab package.
--
Ken Knoblauch
Inserm U846
Stem-cell and Brain Research Institute
Department of Integrative Neurosciences
18 avenue du Doyen Lépine
69500 Bron
France
tel: +33 (0)4 72 91 34 77
fax: +33 (0)4 72 91 34 61
portable: +33 (0)6 84 10 64 10
http:
ou all! This list is pleasure!!!
>
> Marc
>
But, try
all.equal(tt, t)
[1] TRUE
and see the R FAQ 7.31
--
Ken Knoblauch
Inserm U846
Stem-cell and Brain Research Institute
Department of Integrative Neurosciences
18 avenue du Doyen Lépine
69500 Bron
France
tel: +33 (0)4 72 91 34 77
fax: +33 (0
> How can I do that?
>
> I would like to thank you in advance for your help
> Best Regards
> Alex
>
>
> [[alternative HTML version deleted]]
Read section 2.7 of An Introduction to R that comes with the
distribution.
--
Ken Knoblauch
Inserm U846
Stem-cell a
eturn(c(v,w))
> }
>
> # The following for loop works
> result<-data.frame()
> for (i in 1:length(df1[,1])) {
> result<-rbind(result,fcttest(df1[i,1],df1[i,2],df1[i,3]))
why bother with lapply when you can just do this
with(df1, cbind(df1[[1]] * df1[[2]], df1[[2]] + df1[[3]
s to do this but how about
DTA <- cbind(day = dta$day, stack(dta[, -1]))
xyplot(values ~ day | ind, DTA, type = "b", layout = c(2, 6))
for which you can add additional annotations as desired.
By the way, do you realize that you have repeated column
names in your data frame?
Peng, C gmail.com> writes:
>
>
> what is ESP package? Thanks.
I've heard that It's only available over from a repository
accessible through a next-generation
wifi system call oui-ja.
(Beware humor travels poorly over the internet
and across linguistic differences!).
n verify for yourself that a factor yields FALSE here
x <- db1[[1]]
is.vector(x)
[1] FALSE
so I think that this at least explains why it doesn't work as
you expected.
> Thank you for your help.
>
> Kenneth
--
Ken Knoblauch
Inserm U846
Stem-cell and Brain Research Institute
Department o
- function(x, logY = TRUE, ...) {
xy <- ecdflt(x)
if (logY) xy[, 2] <- log10(xy[, 2])
panel.xyplot(xy[, 1], xy[, 2], ...)
}
xyplot(X ~ X | F, panel = function(x, y = NULL, ...){
panel.ecdflt(x, ...)
})
>
> Many thanks,
> Jeff
>
> --
> Jeff Stevens
Roger Koenker uiuc.edu> writes:
> I'm trying to redo an old plot with:
> > sessionInfo()
> R version 2.11.0 Under development (unstable) (2010-02-09 r51113)
> x86_64-apple-darwin9.8.0
> When I do:
>
> pdf("lty.pdf",height = 6, width = 8)
> u <- 1:100/100
> y <- matrix(rep(1:10,each = 100),100)
lass(x)=c('c1','c2')
> test(x)
It works fine for me if you add a default method,
which I think is what it is looking for.
test.default <- function(x){
cat("default")
x
}
test(x)
c1
c2
default[1] 1
attr(,"class")
[1] "c1" "c2"
--
ame output which looks ok.
Just for the record...
Z
R> sessionInfo()
R version 2.10.1 (2009-12-14)
i486-pc-linux-gnu
locale:
[1] C
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] fortunes_1.3-7
On Thu, 11 Feb 2010
ld
the same output which looks ok.
Just for the record...
Z
R> sessionInfo()
R version 2.10.1 (2009-12-14)
i486-pc-linux-gnu
locale:
[1] C
attached base packages:
[1] stats graphics grDevices utils datasets methods base
other attached packages:
[1] fortunes_1.3-7
On Thu, 11 Feb
Hadley Wickham rice.edu> writes:
>
> Hi all,
>
> Is there a fast way to determine the number of lines in a file? I'm
> looking for something like count.lines analogous to count.fields.
>
> Hadley
How about something like
length(readLines(fname))
Ken
_
and in R, it is some uncalibrated combination
of frame buffer values that is being used.
> Best,
>
> baptiste
>
Ken
--
Ken Knoblauch
Inserm U846
Stem-cell and Brain Research Institute
Department of Integrative Neurosciences
18 avenue du Doyen Lépine
69500 Bron
France
tel: +33 (0)
st^ex/(Contrast^fx + sig^fx)),
start = list(Rm = 30, sig = 0.05, ex = 3,
fx = 3.1))
-Peter Ehlers
Ken Knoblauch wrote:
Hi,
I'm getting an error that I don't understand when updating an nls
object. Here is a toy example.
dd <- structure(list(Contrast = c(0.0
onInfo()
R version 2.10.1 Patched (2010-01-25 r51051)
i386-apple-darwin9.8.0
locale:
[1] en_US.UTF-8/en_US.UTF-8/C/C/en_US.UTF-8/en_US.UTF-8
attached base packages:
[1] stats graphics grDevices utils datasets methods
[7] base
Thanks, in advance, for any help.
Ken
--
Ken Knoblauch
Ins
3
823
933
or just vectors
expand.grid(1:3, 1:3)
Var1 Var2
111
221
331
412
522
632
713
823
933
--
Ken Knoblauch
Inserm U846
Stem-cell and Brain Research Institute
Department of Integrative Neurosciences
18
.frame") so that they would still
inherit other data frame methods. My rbind.mlds.df works fine
with them, and I document it accordingly.
HTH.
Ken
--
Ken Knoblauch
Inserm U846
Stem-cell and Brain Research Institute
Department of Integrative Neurosciences
18 avenue du Doyen Lépine
69500 Bro
ons/R.app
to open as many copies of the app as you like. Be careful though,
because these inherit environment variables from the terminal
session, not necessarily the same as those when running the
app from the Finder. I was bitten by that the first time I
tried this.
Ken
--
Ken Knoblauch
Inser
0 0
5 0 0
6 0 0
7 0 1
8 0 1
attr(,"assign")
[1] 1 1
attr(,"contrasts")
attr(,"contrasts")$z
[1] "contr.treatment"
> > thanks
> > Ben Bolker
> >
> >
--
Ken Knoblauch
Inserm U846
Stem-cell
aking, the dichromat package is of great value in avoiding color
choices that about 1% of the population would have trouble discriminating.
Ken
--
Ken Knoblauch
Inserm U846
Stem-cell and Brain Research Institute
Department of Integrative Neurosciences
18 avenue du Doyen Lépine
69500 Bron
Fran
of a function that does this?
>
> Thanks,
> Derek McCrae Norton
Try Rnews Volume 5/1, May 2005,
The Programmer's Niche by John Fox
How Do You Spell That Number?
--
Ken Knoblauch
Inserm U846
Stem-cell and Brain Research Institute
Department of Integrative Neurosciences
18 avenue d
desirability. This package contains
tools to estimate the additive contribution of the n scales
to the judgment by a maximum likelihood method under several
hypotheses of how the perceptual dimensions interact.
--
Ken Knoblauch
Inserm U846
Stem-cell and Brain Research Institute
Department of Integrative
default hypothesis when fitting with lm. Also, the
default link function with the poisson family is log.
So, these are things to take into account in any potential
comparison.
Ken
--
Ken Knoblauch
Inserm U846
Stem-cell and Brain Research Institute
Department of Integrative Neurosciences
18 avenue du
that the responsible person is on vacation...
--
Ken Knoblauch
Inserm U846
Stem-cell and Brain Research Institute
Department of Integrative Neurosciences
18 avenue du Doyen Lépine
69500 Bron
France
tel: +33 (0)4 72 91 34 77
fax: +33 (0)4 72 91 34 61
portable: +33 (0)6 84 10 64 10
http://www.sbri.
Ben Bolker ufl.edu> writes:
> >> > Here is a toy example that illustrates the overshoot of the formula
> >> > \documentclass[12pt]{article}
> >> > \usepackage{geometry}
> >> > \geometry{left=2in,right=2in}
> >> > \begin{document}
> >> > <>=
> >> > op <- options(width = 65, digits = 3)
> >> > ddata
Ben Bolker ufl.edu> writes:
> > In the Sweave output for summary for several types
> > of model objects and also for the comparison of models
> > with anova, I find that that the display of the call(s)
> > or formula does not obey the width option, even with
> > keep.source=TRUE set, so that a lon
B = factor(letters[1:2]),
C = factor(LETTERS[1:5]),
S = factor(paste("S", 1:10, sep = "")),
R = rnorm(10))
mod1 <- lm(R ~ A + B + C + S,
ddataframe)
summary(mod1)
@
\end{document}
I have attached the pdf output.
Thanks,
lx <- levels(x)
ly <- levels(y)
lxy <- union(lx, ly)
xy <- cbind(levels(x)[x], levels(y)[y])
xy <- t(xy)
dim(xy) <- NULL
xy <- factor(xy, levels = lxy)
xy
}
> splice.factor(factor(1:3), factor(4:6))
[1] 1 4 2 5
Doran, Harold air.org> writes:
>
> Ista
>
> There are several functions in the MiscPsycho package that can be sued
> for classical item analysis.
>
Since when is classical item analysis a crime?
No wonder the USA is considered such a litigious society!
Ken
--
Ken
eave spaces around the = than it would be for <-.
>
> the reason being ...?
>
> vQ
>
ergonomy!
--
Ken Knoblauch
Inserm U846
Stem-cell and Brain Research Institute
Department of Integrative Neurosciences
18 avenue du Doyen Lépine
69500 Bron
France
tel: +33 (0)4 72 91 34 77
fax:
8224 .
[4,] . . -1.007848357 -0.1117796 -0.5555834
[5,] . . . -0.6816979 -0.6052127
Ken
--
Ken Knoblauch
Inserm U846
Stem-cell and Brain Research Institute
Department of Integrative Neurosciences
18 avenue du Doyen Lépine
69500 Bron
France
te
Ken Knoblauch inserm.fr> writes:
>
> venkata kirankumar gmail.com> writes:
> > I am trying to parse a vector for caliculating minimum in that vector the
> > vector having values like
> >
> > 1Kontrolle
> > 2 Placebo
> >
olle is being treated as a factor so you
are seeing only the codes of the levels. There is
probably something more elegant, but you need
something like,
as.numeric(sapply(with(dd, strsplit(levels(Placebo)[Placebo], "m")), "[[", 1))
--
Ken Knoblauch
Inserm U846
Institut Cellul
nnectivity
in neural systems.
There are also packages for analysing psychophysical data which
are relevant for behavioral neuroscience,
psyphy, MLDS, sdtalt, etc.
Would there be enough for CRAN TASK VIEW?
Ken
--
Ken Knoblauch
Inserm U846
Institut Cellule Souche et Cerveau
Département Neurosciences
te
a user-specified link and the source of the make.link
function can be useful to.
Ken
--
Ken Knoblauch
Inserm U846
Institut Cellule Souche et Cerveau
Département Neurosciences Intégratives
18 avenue du Doyen Lépine
69500 Bron
France
tel: +33 (0)4 72 91 34 77
fax: +33 (0)4 72 91 34 61
por
Ken Knoblauch inserm.fr> writes:
> Robbert Langenberg gmail.com> writes:
> > I am trying to get a prediction of my GAM on a response type. So that I
> > eventually get plots with the correct values on my ylab.
> > The problem I am encountering now is that I cannot seem
and the by term, something like
model.matrix(~ day:mapID - 1, data = mergeday)
in your case.
I added the appropriate columns into my data frame and
also to the newdata for predict. You can see an example
in the appendix of
http://www.journalofvision.org/8/16/10/
HTH,
Ken
--
Ken
abel1",100), rep("label2",50), rep("label3",150))
> df <- data.frame(as.factor(l), x)
> plot(df)
Just to complete my response,
the documentation for plot.data.frame indicates
For a two-column data frame it plots the second column
against the first by the most a
Hi,
Antje yahoo.de> writes:
>
> Hi folks,
>
> I've just discovered that the following code leads to boxplot
> (surprisingly to me).
> Can anybody explain to me why? Is this documented somewhere? I've never
> consider this option before.
>
> x <- rnorm(300)
> l <- c(rep("label1",100), rep("lab
might try the link mafc.logit(m = 2) defined in the
psyphy package. Continuing with your example,
library(psyphy)
fit <- glm(y ~ x, binomial(mafc.logit(2)),
control = glm.control(maxit = 100)) # default didn't converge
x.ord <- order(x)
lines(x[x.ord], fitted(fit)[x.ord], col = &q
Ken Knoblauch inserm.fr> writes:
>
> Susana Zuloaga hotmail.com> writes:
>
> >
> > Hi all
> >
> > I am one recent user of R and have a few doubts
> > I did a binomial GLM with 3 - factor and now I have to test contrasts to
> > iden
is not incorrect to use them in GLM? there is a way to do
> contrasts between treatments for GLM as a Tukey for the ANOVA?
>
> Susana
see
https://stat.ethz.ch/pipermail/r-help/2003-November/041559.html
--
Ken Knoblauch
Inserm U846
Institut Cellule Souche et Cerveau
Département Neuros
Andrew Barr gmail.com> writes:
> This maybe a newbie question. I have a dataframe
that looks like the sample
> at the bottom of the email. I have monthly
precipitation data from several
> sites over several years. For each site,
I need to extract years that have
> a complete series of 12 mon
Anh Tran ucla.edu> writes:
> I always open more than 1 R console in Windows. I can't figure out a way to
> do this with OS X yet. I need that to utilize the duo core on my desktop.
> How would I do that?
>
Have a look here
https://stat.ethz.ch/pipermail/r-sig-mac/2008-April/004814.html
__
Daren Tan hotmail.com> writes:
> Any better solution than this ?
> sum(strsplit("TCGACAATCGGTAACCCGTCT", "")[[1]] == "G")
Try
table(strsplit("TCGACAATCGGTAACCCGTCT", ""))
A C G T
5 7 8 5
and get all 4 at once.
HTH
--
K
mple, my own psyphy and MLDS).
Ken
--
Ken Knoblauch
Inserm U846
Institut Cellule Souche et Cerveau
Département Neurosciences Intégratives
18 avenue du Doyen Lépine
69500 Bron
France
tel: +33 (0)4 72 91 34 77
fax: +33 (0)4 72 91 34 61
portable: +33 (0)6 84
Megh Dal yahoo.com> writes:
> I have one question on expand.grid() function.
> When I write following syntax :expand.grid(c("u", "l"),
>c("u", "l"), c("u", "l")) I get following as
> desired :
> Var1 Var2 Var3
> 1uuu
> 2luu
> 3ulu
> 4llu
> 5
Ken Knoblauch inserm.fr> writes:
> Daren Tan hotmail.com> writes:
> > I tried aggregate, apply etc, but can't get the right result.
> do.call("expand.grid", rep(list(c("u", "l")), 3))
> Var1 Var2 Var3
> 1uuu
> 2l
t;[3,] "cc" "C|E"
>
How about
do.call("expand.grid", rep(list(c("u", "l")), 3))
Var1 Var2 Var3
1uuu
2luu
3ulu
4llu
5uul
6lul
7ul
Armin Goralczyk gmail.com> writes:
> In a function I have a plot and want to add symbols/text only when
> indicated by a logical vector (which was generated by the function
> before, not manually like in the following example):
>
> plot(1:10, 1:10)
> lv <- c(T,T,T,F,F,F,T,T,T,F)
> text(1:10, 1:10
Hi,
Charles Annis, P.E. StatisticalEngineering.com> writes:
> logit.FC <- function(POD.floor = 0, POD.ceiling =1)
> { if (POD.floor < 0 | POD.floor > 1) stop ("POD.floor must be between zero
> and one.")
> if (POD.ceiling < 0 | POD.ceiling > 1) stop ("POD.ceiling must be
> between zero and
andy dsl.pipex.com> writes:
> I am trying to import an *.xls spreadsheet into R. I am doing this as
> follows:
> > read.table(file("A5_DL.xls"))
> So I copied it all over to a text document and tried to import that, thus:
> > read.table("A5.txt")
> The error I got then was:
> Error in scan(file
Johannes Hüsing huesing.name> writes:
>
> Am 02.03.2008 um 17:44 schrieb Gabor Csardi:
>
> > I'm not a statistician, but do i remember well that among all
> > distributions with a given mean and variance, the normal distribution
> > has the highest entropy? This is good enough for me to call it
csiro.au> writes:
>
> y <- sort(rnorm(20)) # say...
>
> m <- s <- numeric(19)
>
> for(i in 2:20) {
> m[i-1] <- mean(y[1:i])
> s[i-1] <- sd(y[1:i])
> }
> -Original Message-
> On Behalf Of cvandy
> Subject: [R] Loop with variable index
> I have a list of 20 values. The fir
hadley wickham gmail.com> writes:
<<>>
>For example, I have this alpha function in
> ggplot:
> alpha <- function(colour, alpha) {
> col <- col2rgb(colour, TRUE) / 255
> col[4, ] <- rep(alpha, length(colour))
> new_col <- rgb(col[1,], col[2,], col[3,], col[4,])
> new_col[is.na(colour)] <- N
en what about
http://finzi.psych.upenn.edu/R/doc/html/packages.html
which allows you to go in a little deeper and see the help pages of each
function
and
http://cged.genes.nig.ac.jp/RGM2/index.php
HTH
ken
--
Ken Knoblauch
Inserm U846
Institut Cellule Souche et Cerveau
Département Neu
Dimitris Rizopoulos med.kuleuven.be> writes:
>
> the 'col' argument of plot() may also be a vector, e.g.,
>
> x <- sort(rnorm(15))
> y <- rnorm(15)
> plot(x, y,
> col = rep(c("red", "blue", "magenta"), c(5, 6, 4)))
> Dimitris
>
>
> Dimitris Rizopoulos
>
> - Original Message
uses the quantile function
and then looking at fivenum to see that it did not.
Looking at the help for fivenum
led me to boxplot.stats where I was that it w
as not necessarily doing the same thing.
HTH
--
Ken Knoblauch
Inserm U846
Institut Cellule Souche et Cerveau
Dép
ts for most measurements, so one day you may be able to
> do that in whatever passes for par() by then.
>
> Duncan Murdoch
For anyone interested, there is a larger selection of units online here,
as well as other places,
http://en.wikipedia.org/wiki/SI_prefix
upon seeing which, I wo
Julien Barnier no-log.org> writes:
>
> Hi,
> When I try this example at work, I got an error message during Sweave
> which is close from something like "warning : semitransparency not
> supported by this device". When I try at home (with the given
> sessionInfo), I got no warning but in both the
Aleksi Lehtonen gmail.com> writes:
> I am trying to estimate several non-linear models simultaneously. I don't
> want to use non-linear mixed model, but non-linear model with same form, but
> it should be estimated separately according to variable group (I have lots
> of groups that have lots of
99 matches
Mail list logo