5
[19] scales_0.4.0 geosphere_1.5-5 colorspace_1.2-6 labeling_0.3
stringi_1.1.1 lazyeval_0.2.0munsell_0.4.3 rjson_0.2.15
This is not an urgent issue, I just think it's curious, so it would be nice to
understand why it happens.
Thanks,
Colin
___
Dear Duncan,
Thank you for your quick reply. I've got the basic version of what I'm
looking for now (see below). My next step will be your rgl::persp3d suggestion
for the hidden lines control.
Best,
Colin
library(MASS)
X <- mvrnorm(1000,mu=c(0,0),Sigma=matrix(c(1,0,0,
0,xlab="x_b",ylab="x_a",zlab="f") ->res
Any suggestions are very appreciated.
Thank you,
Colin
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/
Hi all,
are there any packages to perform a market simulation with the conjoint
analysis' results?
Thanks,
Colin
--
View this message in context:
http://r.789695.n4.nabble.com/Conjoint-Analysis-in-R-tp842239p4651095.html
Sent from the R help mailing list archive at Nabbl
unning the function against each
row of the table it fails
mydataframe["TotalDif"] <- difcount(surfA,surfB)
Error in file(fname, "r") : invalid 'description' argument
This error seems to be coming from the readAsciiGrid function, but its
being given the
the labels are incorrect. The labels cycle through CA-I,
CAIII, CA-II, in that order.
The plot file is attached.
Colin
On Mon, Feb 13, 2012 at 1:31 AM, Jim Lemon wrote:
> On 02/13/2012 09:51 AM, Colin Wahl wrote:
>>
>> Does anyone have any recommendations for producing dotplots wit
ve double checked
that all the numeric columns in the data frame are not factors, and
the error persists.
I'd really prefer not just calling it quits and resorting to
old-school sloppy bar charts, but if thats what I need to do to finish
this in a timely manner, then so be it.
Thank you,
Col
a why R cant do vector functions for these data?
Solution?
Thanks,
Colin Wahl
M.S. Biology candidate
Western Washington University
fit<-lm(sandcomb ~ CCEC25)
z<-predict(fit, data.frame(CCEC25=range(CCEC25)))
lines(range(CCEC25), z, lty=2, lwd=0.75, col="grey51")
> is
value?
A description of the data and the R code and output is attached for
anyone interested.
Best wishes,
Colin Aitken
---
David Winsemius wrote:
On Nov 7, 2011, at 12:59 PM, Colin Aitken wrote:
How does R estimate the intercept term \alpha in a loglinear
model with
+ \beta), exp(\alpha + \gamma), exp(\alpha +
\beta + \gamma). The MLE of \alpha appears to be log(n_{.1} *
n_{1.}/n_{..}), but this is not equal to the intercept given by R in the
example I tried.
With thanks in anticipation,
Colin Aitken
--
Professor Colin Aitken,
Professor of Forensic Statistics
on samples(rows))?
Thank you,
Colin Wahl
Graduate student,
Western Washington University
code & background:
I am looking at the ordination of abiotic stream variables between different
sampling locations.
abiot.pca=prcomp(all24[, c(10, 13:18)], retx=TRUE, center=TRUE, scale=
TRUE
Yes that sounds very much like it. I'm using R 13.1 on Windows 7 32 bit. I did
do a search but didn't find anything - probably not using the correct search
terms. Glad its not just me.
Thank you,
Col.
From: Michael Sumner
Cc: Peter Morgan ; "r-help@r-project
also just tried the image command and
that was giving the same result as well.
What I'm wanting are blocks together without the white lines going through but
just cant seem to find a way to do that?
Best regards,
Col.
From: Peter Morgan
To: Colin Ford
Se
Hello,
I must start by saying that I am an R novice and am sorry if this is a
no-brainer...
I have an csv file that contains a grid of 300x300 data points. Each point
represents
a 1Km square on a map. Each point is either a floating point number or NA. I
load the
data in with:
data <- read.
might come from the fact that I do have DD format in the LatLong
file. Maybe I should be able to specify that it is zone 11 north. What do
you think?
Colin
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.e
"inside" in the
barplot function would be the way to go but I am lacking some time to
acquire the skills required to implement it.
Colin
On Fri, Jul 29, 2011 at 11:05 AM, Jannis wrote:
> Why dont you plot a stacked barplot with the two values for each depth (or
> whatever is
not yet implemented.
Colin
On Fri, Jul 29, 2011 at 9:02 AM, Jannis wrote:
> Dear Collin,
>
>
> as always, a reproducible example code would help us to understand what you
> want to do. This way I can only guess
>
> And my guess would be that it is much easier to use:
>
; would allow not to plot the line
dividing two bars and using a null color within these bars would allow
superposition with the actual bar plot and indicate the sample depth.
Do you have any suggestions or experience to share?
Colin
[[alternative
ut$N.upper, as.numeric(y), nut$N.lower, as.numeric(y), lty=1
, col=1)
})
I have no experience with panel functions and would very much appreciate
advice.
Thank you,
Colin Wahl
Graduate Student
Dept. of Biology
Western Washington University
[[alternative HTML version deleted]]
I would rather use cleveland dot plots than bar charts to display my study
results. I have not been able to find (or figure out) an R package that is
capable of producing the publication quality dot charts Im looking for. I
have either not been able to get error bars (lattice), cannot order the dat
Any advice on how to do this?
Thanks,
Colin Curry
--
Colin Curry
Ph.D. Candidate
Canadian Rivers Institute
Department of Biology
University of New Brunswick, Fredericton
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
h
me.
Many thanks again,
Colin.
-Original Message-
From: Uwe Ligges [mailto:lig...@statistik.tu-dortmund.de]
Sent: 14 January 2010 14:06
To: Colin Millar
Cc: Gabor Grothendieck; r-help@r-project.org
Subject: Re: [R] R package dependencies
For the original question:
> what are a packa
ect.org/bin/windows/contrib/2.10> "
destdir <- "C:/colin/SCOTS-user-group/pkgs"
lab.pkgs <-
c("gamair","mgcv","survival","lme4","nlme","lattice","MASS","nnet","splines","stats4&quo
have listed as buisness needs, i would like to include all
reverse dependencies of this collection that do not have dlls.
I hope this makes sense (the question not the reason).
Kind regards,
Colin.
[[alternative HTML version deleted]]
__
R
I've written a very simple bit of code to plot a trajectory using rgl:
x <-
(c(0,-5.947,-11.496,-16.665,-21.474,-25.947,-30.116,-34.017,-37.684,-41.148,-44.435,-47.568,-50.567,-53.448,-56.223,-58.906))
y <-
(c(0,33.729,65.198,94.514,121.784,147.151,170.797,192.920,213.717,233.360,252.003,269
@r-project.org [mailto:r-help-boun...@r-project.org]
On Behalf Of Colin Millar
Sent: 27 November 2009 16:41
To: Anastasia; r-help@r-project.org
Subject: Re: [R] Simple Function doesn't work?
Hi,
You would also make your code more efficient and possible more readable
by doing
ReturnsGrid
])
/ m)
drop(out)
}
(xx <- ReturnsGrid(0, 9, 3))
#[1] 0 3 6
(xx <- ReturnsGrid(0:2, 9:11, 3))
#[1,]012
#[2,]345
#[3,]678
But it seems like you could also do it using sequence ...
seq(x, y-1, by = m)
HTH,
Colin
-Original Message-
From: r-help-b
n", "dx", "df")) {
command <- match.arg(command)
switch(command,
scrn = cat("scrn :Screening","\n"),
dx = cat("dx:Diagnosis","\n"),
df = cat("df:Don't Forget","\n")
)
}
Colin.
I think you just missed some commas out...
aar <-
function(command = c("scrn", "dx", "df"))
{
command <- match.arg(command)
switch(command,
scrn = cat("scrn :Screening","\n"),
dx = cat("dx:Diagnosis","
Hi,
try ?table
# for example
(s3 <- table(s))
# and if you want a single value
s3["a"]
# or
s3[1]
HTH,
Colin.
From: r-help-boun...@r-project.org on behalf of Marc Giombetti
Sent: Tue 17/11/2009 22:55
To: r-help@r-project.org
Subje
-help-boun...@r-project.org]
On Behalf Of Colin Millar
Sent: 17 November 2009 16:10
To: Zhiyuan Jason ZHENG; r-h...@stat.math.ethz.ch
Subject: Re: [R] Calculating the power of a negative number
Hi,
Look at
?NumericConstants
At the bottom of the details section you will find:
"Note t
this helps,
Colin.
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
On Behalf Of Zhiyuan Jason ZHENG
Sent: 17 November 2009 14:00
To: r-h...@stat.math.ethz.ch
Subject: [R] Calculating the power of a negative number
Hello,
I use R a lot, one thi
ylab = "Frequency", xlab = "Score",
xlim = c(1,5), ylim = c(0,100),
main = colnames(data)[i]
)
}
print(plots[[1]])
# or export
for (i in seq_along(plots))
{
png(paste("hist", i, ".png", sep = ""))
print(plots[[i]])
dev.off()
* sum_i
[ (y_i - (a + b*x_i))^2 ]
You can simplyfiy this function if you condition on the variance and intercept.
You are then left with a simple function coresponding to the log-likelihood
curve.
This function is trivial to minimise.
Hope this helps.
Colin
Thanks! The formattingOut.odt document had exactly what I needed. The style
"highlight" wasnt defined in my environment.
Cheers,
Colin
On Wed, Nov 4, 2009 at 4:59 AM, Max Kuhn wrote:
> It's hard to say without a reproducible example or the output form
> sessionInfo().
et as I want. What are valid values that
tableStyles$cell can be set to? Does this have to be a custom style defined
in the ODT file beforehand?
Thanks in advance,
Colin Robertson
PhD Candidate
Dept of Geography
University of Victoria
[[alternative HTML version deleted]]
___
A quick google on 'normality test' (no quotes) gives
http://en.wikipedia.org/wiki/Normality_test. This gives you a few more tests
than the KS test.
Cheers,
Colin.
Steve Lianoglou wrote:
> Hi,
>
> I think you can also use a qq-plot to do the same, no? You won't
Hi,
This should work - head is quite a usefull summary function
head(df.mydata[df.mydata$A=="X" & df.mydata$C < 2, ],3)
Colin.
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
On Behalf Of Noah Silverman
Sent: 26 Aug
Or perhaps use a temporary vector might be neater?
tmp <- with(df.mydata, B[A=="X" & C < 2])
df.mydata[order(tmp) %in% 1:3,] # gives df with highest three values of B
or
head(df.mydata[order(tmp),],3) # gives first 3 rows of df sorted by B
Colin.
-Original Message---
ge, but for
various reasons I don't want to (I need to run the process many, many times on
very large datasets). Thanks for any pointers - I'm being very slow today...
Colin
PS sessionInfo()
R version 2.8.1 (2008-12-22)
i386-pc-mingw32
locale:
LC_COLLATE=English_United Kingd
ase
other attached packages:
[1] rgl_0.83-3 misc3d_0.6-1 debug_1.1.0mvbutils_1.1.1 svSocket_0.9-5
svIO_0.9-5 R2HTML_1.58svMisc_0.9-5 svIDE_0.9-5
loaded via a namespace (and not attached):
[1] tools_2.8.1
Thanks!
Colin
--
Please note that the views expressed in this
;JostD' for all combinations
of columns, perhaps returning a matrix of distances. Thanks for any help or
suggestions.
Cheers
Colin
Function:
JostD <- function(DF, X, Y) {
Ni1 <- DF[,X]
Ni2 <- DF[,Y]
N1 <- sum(Ni1)
N2 <- sum(Ni2)
pi1 <-Ni1/N1
pi2
I have one other question regarding the nnet package in R:
What is the value of the step size (often referred to as "eta") for
the gradient descent function?
I dont see anywhere in the API where this value can be modified.
Thanks,
-Colin
___
4) Will this package ever be updated to allow for multiple layers
instead of just one? (just out of curiousity)
I have to present this paper on Friday May 2nd so I would greatly
appreciate a timely response.
Thanks,
-Colin
__
R-help@r-project.or
FactorA2:FactorB1 FactorA2:FactorB3
1 1 0 0
2 1 0 0
3 1 0 0
4 1 1 0
5 1 0 0
6 1 0 1
If anyone has any advice o
to calculate, for example, the appropriate 95% confidence
intervals for the predicted values of autosome=ore + mtDNA=ore,
autosome=ore + mtDNA=simw501, etc.
Sincerely,
Kristi Montooth and Colin Meiklejohn
Ecology and Evolutionary Biology
Brown University
__
odel
results.
Any help / clarification appreciated,
Colin
Colin Robertson
University of Victoria
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE
class.
Any help appreciated,
Thanks
Colin
Colin Robertson
Dept of Geography
University of Victoria
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r
48 matches
Mail list logo
