Hi
The function loess takes very long time if the dataset is very huge
I have around 100 records
and used only one independent variable. still it takes very long time
Any suggestion to reduce the time
-
Thanks in Advance
Arun
--
View this message in context:
http://r.789695.n4.
On Thu, 12 Apr 2012, jpm miao wrote:
It seems to only works for zoo or ts objects?
It was tested for zoo, ts, and plain numeric vectors.
I tried to run it for xts objects, and error message occurs.
I had a quick look if this can be avoided. But xts' merge() method is not
fully compatible
On Thu, 12 Apr 2012, jpm miao wrote:
Example:
Will ts objects be obsolete or modified?
No, I don't think so.
a[,1]
1983 Q1 2.747365190
1983 Q2 2.791594762
1983 Q3 -0.009953715
1983 Q4 -0.015059485
1984 Q1 -1.190061246
1984 Q2 -0.553031799
1984 Q3 0.686874720
1984 Q4 0.95
On Wed, Apr 11, 2012 at 8:12 PM, Gabor Grothendieck
wrote:
> On Wed, Apr 11, 2012 at 10:17 PM, Whit Armstrong
> wrote:
>> I must admit I'm a little ashamed to have been using R for so long,
>> and still lack a sound understanding of deferred calls, eval, deparse,
>> substitute, and friends.
>>
>>
At 12:03 12/04/2012, you wrote:
>I had the same problem! So, as I'm a linux user,
>I prefer use linux terminal. On terminal I type this to compile
>
>R CMD Sweave --encoding=utf-8 myfile.Rnw
>
>and the compilation is successful. Try to set the encoding option in Sweave().
>
>Bests.
>Walmes.
>
>=
It seems to only works for zoo or ts objects?
I tried to run it for xts objects, and error message occurs. Once I
coerce it to be an zoo object (by as.zoo), it works.
Error message:
Error in model.frame.default(formula = dynformula(PIh - PI ~ L(X, 0:i) + :
variable lengths differ (found for
If I calculate the C-index after 10-fold cross-validation based on Cox
model, which statistic test
should I used to compare two C-statistc values.I means I use
cross-validation to obtain a predictive probability for each individual,
then combine all the probabilities to calculate the C-index.If I
Is there some way to remove superscripts from objects returned by
html/xmlParse (XML package)?
h <- "CataDog"
doc <- htmlParse(h)
xpathSApply(doc, "//p", xmlValue)
[1] "Cata" "Dog"
I could probably remove the tags from the "h" object above, but I'd
rather just work with the results from htmlPa
Hi,
I am working on an implementation of R within some web application (R is
just a part of a larger project, so R has to be incorporated within an
existing setup)
Here is the scenario, PHP grabs a csv dataset from some server and turns it
into a 2 dimensional array (with each row being an arra
Hi, I'm trying to calculate a ratematrix for a RNA aligment (U instead of T)
in order to use it as a ratematrix in Profidst (a phylogenetic program which
takes into account both the primary sequence and the secondary structure of
the RNA, in my case rRNA). The sequence-structure aligment has been m
I am trying to schedule my R script using cron, but it is not working. It
seems R can not find packages in cron. Anyone can help me? Thanks.
The following is my bash script
# source my porfile
. /home/winie/.profile
# script.R will load packages
R CMD BATCH /home/script.R
--
View this message
What makes you think those are the only two options?
I happen to prefer the xts convention, but there is a lot of code out there
that successfully uses ts just as it is, and I can't see breaking all of that
to meet an arbitrary preference of sign convention. (It isn't my decision
anyway...)
---
Yes, this is as documented. See ? lag.xts under details for the
justification and how to change the default if desired.
Michael
On Wed, Apr 11, 2012 at 11:13 PM, jpm miao wrote:
> BTW, zoo is like ts in the application of lag.
> In other words, zoo and xts are opposite in this issue.
>
> 2012/4/
BTW, zoo is like ts in the application of lag.
In other words, zoo and xts are opposite in this issue.
2012/4/12 jpm miao
> Example:
> Will ts objects be obsolete or modified?
>
> > a[,1]
> 1983 Q1 2.747365190
> 1983 Q2 2.791594762
> 1983 Q3 -0.009953715
> 1983 Q4 -0.015059485
On Wed, Apr 11, 2012 at 10:17 PM, Whit Armstrong
wrote:
> I must admit I'm a little ashamed to have been using R for so long,
> and still lack a sound understanding of deferred calls, eval, deparse,
> substitute, and friends.
>
> I'm attempting to make a deferred call to a function which has defau
Example:
Will ts objects be obsolete or modified?
> a[,1]
1983 Q1 2.747365190
1983 Q2 2.791594762
1983 Q3 -0.009953715
1983 Q4 -0.015059485
1984 Q1 -1.190061246
1984 Q2 -0.553031799
1984 Q3 0.686874720
1984 Q4 0.953911035> lag(a,4)[,1]
1983 Q1 NA
1983
I must admit I'm a little ashamed to have been using R for so long,
and still lack a sound understanding of deferred calls, eval, deparse,
substitute, and friends.
I'm attempting to make a deferred call to a function which has default
arguments in the following way:
call.foo <- function(f) {
I had the same problem! So, as I'm a linux user, I prefer use linux
terminal. On terminal I type this to compile
R CMD Sweave --encoding=utf-8 myfile.Rnw
and the compilation is successful. Try to set the encoding option in
Sweave().
Bests.
Walmes.
===
Hi
I ran the following sweave file in R2.14.1 and upgraded to R2.15
yesterday with the same setup
I got the following error message when I rand the following Sweave file
> setwd("D:/Cic/Sweave/Parasite/Comb/12")
> Sweave("D:/Cic/Sweave/Parasite/Comb/12/ParasiteComb12.Rnw")
Error: c("'ParasiteC
Thanks again.
By the way, any idea why I get this error:
merged3 <- merge_all(list_of_files , by = "Name")
Error in `[.data.frame`(df, , match(names(dfs[[1]]), names(df))) :
undefined columns selected
> sessionInfo()
R version 2.14.2 (2012-02-29)
Platform: i386-pc-mingw32/i386 (32-bit)
loca
I repeat myself: Please reply all to the list and don't send HTML.
I'm not aware of facilities to do that automatically (though others
might now them) -- you can draw arrows manually with the arrows()
function.
Michael
On Wed, Apr 11, 2012 at 7:08 PM, John Kim wrote:
> thank you so much..
>
> i
Alternatively, you can use
curve(x^3, from = -5, to = 5); abline(h = 0, v = 0, lty = 2)
which will work even if the axes aren't in the middle of the image.
Michael
On Wed, Apr 11, 2012 at 6:54 PM, R. Michael Weylandt
wrote:
> Please reply all to the list and don't send HTML.
>
> curve(x^3, fro
Please reply all to the list and don't send HTML.
curve(x^3, from = -5, to = 5); grid()
If you only want dotted lines at x = 0 & y = 0 use
grid(2,2)
instead.
Michael
On Wed, Apr 11, 2012 at 6:49 PM, John Kim wrote:
> thanks for the reply..
>
> i tried the code below.. it shows the graph.. bu
On Wed, Apr 11, 2012 at 06:04:28AM -0700, capy_bara wrote:
> Dear all,
>
> I want to explore the nullspace of a matrix S: I currently use the function
> Null from the MASS package to get a basis for the null space:
> > S = matrix(nrow=3, ncol=5, c(1,0,0,-1,1,1,1,-1,-1,0,-1,0,0,0,-1)); S
> > MASS:
On 11-Apr-2012 16:28:31 Christopher Kelvin wrote:
> Hello,
> can i implement this as 10% censored data where t gives me
> failure and x censored.
> Thank you
>
> p=2;b=120
> n=50
>
> set.seed(132);
> r<-sample(1:50,45)
> t<-rweibull(r,shape=p,scale=b)
> t
> set.seed(123);_
> cens <- sample(1:50,
On Wed, Apr 11, 2012 at 06:04:28AM -0700, capy_bara wrote:
> Dear all,
>
> I want to explore the nullspace of a matrix S: I currently use the function
> Null from the MASS package to get a basis for the null space:
> > S = matrix(nrow=3, ncol=5, c(1,0,0,-1,1,1,1,-1,-1,0,-1,0,0,0,-1)); S
> > MASS:
On Apr 11, 2012, at 4:41 PM, John Kim wrote:
can anybody tell me how i can draw x- y- axis and draw x^3 graph using
R graph??
i need nice coordinate system with legends and coordinate
numberings..
Look at:
?curve
and nice graph of x^3 on it.. it will be nice if you tell me how i
can ce
The easiest way is to just use ?curve (type ?curve at the prompt to
get documentation for curve): e.g., curve(x^3, from = -5, to = 5)
You could also build the plot yourself like:
x <- seq(-5, 5, length.out = 200)
y <- x^3
plot(x,y)
Michael
On Wed, Apr 11, 2012 at 4:41 PM, John Kim wrote:
> ca
can anybody tell me how i can draw x- y- axis and draw x^3 graph using
R graph??
i need nice coordinate system with legends and coordinate
numberings..
and nice graph of x^3 on it.. it will be nice if you tell me how i
can center the graph..
i want the origin (0,0) to be right in the middle of
On Sun, Apr 1, 2012 at 9:16 AM, Till Bayer wrote:
> Hi all!
>
> I want to generate a heat map from an all-vs-all comparison. I have the
> data, already scaled to 0-1. However, I have the values only for the
> comparisons in one way, and not for the comparisons between the same group
> (which are a
On 11/04/2012 19:33, Křištof Želechovski wrote:
Dnia środa, 11 kwietnia 2012 13:26:43 Prof Brian Ripley pisze:
On 11/04/2012 12:42, Duncan Murdoch wrote:
On 12-04-11 12:43 AM, Křištof Želechovski wrote:
http://finzi.psych.upenn.edu/R/library/datasets/html/morley.html>
"The classical data of M
> What I would have liked is something like a cloud of lines, similar to what
> I get when I convert the data into a matrix (why do I not just use a matrix?
> I come from MATLAB and this seems natural, however, my data is large and a
> data frame seems to be an advantageous way to handle that).
It
On Wed, Apr 11, 2012 at 7:16 AM, David Winsemius wrote:
>
> On Apr 11, 2012, at 9:03 AM, David Winsemius wrote:
>
>>
>> On Apr 11, 2012, at 6:28 AM, maxbre wrote:
>>
>>> hi, I just realised I want to go a little further in the control of the
>>> chart
>>> appearance and I would like to have the sa
And it's now fixed in the dev version.
Hadley
On Tue, Mar 13, 2012 at 11:37 AM, Helios de Rosario
wrote:
> Michael,
>
> Thanks for the pointer to the discussion in the ggplot list. It seems
> that the reason of this behaviour of facet_grid() is already known and
> being discussed by the developer
Thank you Ilai. Problem solved. There is a small detail, alpha affects the
rug and the curve line opacity too. Is possible to specify it just to
polygon?
Bests.
==
Walmes Marques Zeviani
LEG (Laboratório de Estatística e Geoi
Dnia środa, 11 kwietnia 2012 13:26:43 Prof Brian Ripley pisze:
> On 11/04/2012 12:42, Duncan Murdoch wrote:
> > On 12-04-11 12:43 AM, Křištof Želechovski wrote:
> >> >> http://finzi.psych.upenn.edu/R/library/datasets/html/morley.html>
> >>
> >> "The classical data of Michaelson and Morley on the
Hello all~
I am interested in clarifying something more conceptual, so I won't be
providing any data or code here.
>From what I understand, partial dependence plots can help you understand the
relative dependence on a variable, and the subsequent values of that
variable, after "averaging out th
densityplot(~y|B, groups=A, data=dt,
plot.points="rug",
col=trellis.par.get("superpose.polygon")$col, alpha=.5,
panel=panel.superpose,
panel.groups=my.panel.densityplot)
Worked for me (i.e. semi-transparent superpose.polygon colors). Is
that not wha
On Apr 11, 2012, at 16:51 , John C Nash wrote:
> nls() often gives this message, which is misleading in that it is the
> Jacobian that is not
> of full rank in the solution of J * delta ~ - residuals or in more
> conventional
> Gauss-Newton J' * J delta = -g = - J' * residuals. My view is
David,
Right you are! Thanks for pointing that out.
strsplit <- 1:10
strsplit("With spaces", NULL)
strsplit
Jean
David Winsemius wrote on 04/11/2012 01:17:07 PM:
> [image removed]
>
> Re: [R] strsplit help
>
> David Winsemius
>
> to:
>
> Jean V Adams
>
> 04/11/2012 01:19 PM
>
> Cc:
On Apr 11, 2012, at 2:01 PM, Jean V Adams wrote:
Alison,
Your code works fine on the first six lines of the data that you
provided.
Rumino_Reps_agreeWalign <- data.frame(
geneid = c("657313.locus_tag:RTO_08940",
"457412.251848018",
"657314.locus_tag:CK5
FYI,
whenever getting an error, run traceback() *immediately after* (i.e.
before any other commands) and include that in your error report.
Also include the output of sessionInfo(). It helps tremendously and
spares lots of second guessing.
/Henrik
On Wed, Apr 11, 2012 at 8:28 AM, krtek wrote:
Hello,
I'm doing some graphics for a paper and a need customize such with filled
region above the density curve. My attempts I get something very near what
I need, but I don't solve the problem of use semitransparent filled. Below
a minimal reproducible code. Someone has any idea?
require(lattice
Alison,
Your code works fine on the first six lines of the data that you provided.
Rumino_Reps_agreeWalign <- data.frame(
geneid = c("657313.locus_tag:RTO_08940",
"457412.251848018",
"657314.locus_tag:CK5_20630",
"657323.locus_tag:CK1_330
Thanks a ton, Michael! Everything is running much faster now!!
Wyatt
-Original Message-
From: R. Michael Weylandt [mailto:michael.weyla...@gmail.com]
Sent: Wednesday, April 11, 2012 1:41 PM
To: Wyatt McMahon
Cc: R-help@r-project.org
Subject: Re: [R] multicore/mcparallel error
Parallel i
Parallel is one of those packages (like tools or grid) that is not
loaded by default, but comes with the standard installation.
Running library(parallel) will make mclapply() available. Then (at
least for me) it's easy to parallelize some code:
system.time(lapply(1:5, function(x) max(rnorm(50
Hello,
I've got a small dataset on box turtle shell measurements that I would like to
perform a detrended correspondence analysis on. I thought that it would be
interesting to examine the morphometrics for each species in the area of
overlap and in areas where neither species occurs.
I've tak
Can you give the line of code that gives the error? That'd make it
much easier to see what's there that should/shouldn't be.
Michael
On Wed, Apr 11, 2012 at 11:28 AM, krtek wrote:
> Thank you! Updating to R-patched really helped me.
>
> My problem were not been into my code. I've tried to run s
Hello:
I have looked at the bnlearn and deal packages for infering bayesian
network. Can anyone suggest any other suitable package for constructing
bayesian gene regulatory network using gene expression data?
Thanks!
John
[[alternative HTML version deleted]]
Thank you! Updating to R-patched really helped me.
My problem were not been into my code. I've tried to run source() and the
error had occurred again. I didn't use a mixture of different versions of R,
version 2.13.2 has been deleted.
--
View this message in context:
http://r.789695.n4.nabble.
Oh yes, I see now the problem...
thank you
max
--
View this message in context:
http://r.789695.n4.nabble.com/xyplot-lattice-fine-control-of-axes-limits-and-thick-marks-with-log-scale-tp4511897p4549180.html
Sent from the R help mailing list archive at Nabble.com.
___
Unfortunately, this won't help for expression arrays. Last time I checked,
those IDATs appeared to be encrypted.
crlmm, methylumi, and minfi can all read IDAT files... *if* they are
"version 3" or later (e.g. genotyping, methylation, etc).
Otherwise you are probably stuck with GenomeStudio if it
Hi Ekta,
You can use the bioconductor package methylumi to read in .idat files. The
function methylumIDAT() will read in .idat files and provide you with beta
values along with the M and U values as a MethyLumiSet object. You can then use
the functions methylumi.bgcorr() and normalizeMethyLumiS
Hello,
can i implement this as 10% censored data where t gives me failure and x
censored.
Thank you
p=2;b=120
n=50
set.seed(132);
r<-sample(1:50,45)
t<-rweibull(r,shape=p,scale=b)
t
set.seed(123);
cens <- sample(1:50, 5)
x<-runif(cens,shape=p,scale=b)
x
Chris Guure
Researcher,
Institute for
Hi,
A line can be added to an existing plot using the abline function. For
example, if a is the intercept and b is the slope, the command would be
abline(a=a, b=b)
To overlay a new plot on an existing one, use the command: par(new=TRUE).
For example:
plot(1:10, 1:10)
par(new=TRUE)
plot(log
Hello,
>
> I would like to know for each variable what is date format?
> Is it “ %y-%d-%m “ or “ %m-%d-%y “
>
It's 4 digits, with uppercase 'Y'
Try
Data <- data.frame(var1=c("12-15-2011", "12-15-2011", "12-17-2001"),
var2=c("2001-15-12", "2001-15-01", "2001-15-01")
Dear all,
I want to use string split to parse column names, however, I am having
some errors that I don't understand.
I see a problem when I try to rbind the output from strsplit.
please let me know if I'm missing something obvious,
thanks,
alison
here are my commands:
>strsplit<-strsplit(as
nls() often gives this message, which is misleading in that it is the Jacobian
that is not
of full rank in the solution of J * delta ~ - residuals or in more
conventional
Gauss-Newton J' * J delta = -g = - J' * residuals. My view is that the
gradient itself
cannot be "singular". It's just
Hello R users,
I am analizing survival data (mostly uncensored) and want to extract the
most out of it.
Since I have more than one factor, I´ve read that the survival regression
can help to test the interactions between factors, and then decide how to do
the comparisons using the Log-rank test (sur
This has nothing to do with your question, but instead of using
>class=c(rep(1,3),rep(2,3),rep(3,3))
>
It's probably easier to use class = rep(1:3, each =3)
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat
Michael,
Thanks for your help!
I may not be understanding you. I've upgraded to 2.15. From your text
below, I thought the parallel package was included with R 2.15, but I
cannot find an mclapply function there. One exists through the multicore
package, which I've now installed. I also went loo
Hi Michael,
Thank you for your help!
I did some googling and researching... Reading the following article,
http://www.ecd.bnl.gov/pubs/BNL-79819-2008-JA.pdf
It seems that once we estimate the parameters of the bivariate normal
distribution,
then we can plug into the formula of conditional distrib
One of the products of Project MOSAIC (funded by an NSF CCLI grant)
has been the development of an R package with the goal of making it
easier to use R, especially in teaching situations. We're not quite
ready to declare that we've reached version 1.0, but version 0.4 does
represent a fairly large
Hi all,
I appreciate that if someone please help me with an R solution/code/library
for scraping Google Insight data. It seems difficult to get hold of the URL
corresponding to the resulting CSV file that is appropriate for
*read.csv*or another similar function.
Best, R
[[alternative HTML
In addition to what Jeff and David have said...
If you really want to create a separate data frame for each subgroup then
you need use the assign function, and also data.split[[i]] instead of
data.split[1].
for (i in 1:3) {
assign( paste('sub',i,sep='') , data.split[[i]]
}
Jeff showed you how
Given the caveats Ted describes here:
http://tolstoy.newcastle.edu.au/R/help/05/06/5992.html it seems that
bootstrapping might be the only way to get (somewhat) credible
prediction intervals: the boot package on CRAN can help to facilitate
getting these. Here's some documentation for CI:
http://ww
I use ssh -Y in this situation.
You can test with other X client applications, such as xclock. Try running
xclock after your login. If you get the same "unable to open" message then
it's not an R problem.
-Don
--
Don MacQueen
Lawrence Livermore National Laboratory
7000 East Ave., L-627
Livermo
On 11.04.2012 15:27, Prof Brian Ripley wrote:
On 11/04/2012 14:08, R. Michael Weylandt wrote:
Well, obviously the interface to .Internal(identical) changed between
2.13.x and 2.15 -- as, ?.Internal says: Only true R wizards should
even consider using this function...
Anyways, find where you u
Hi all,
Are there functions in R that could help me do the following?
We have a special type of regression which is called Geometric Mean
Regression.
We have done some search and found the following:
https://stat.ethz.ch/pipermail/r-help/2011-July/285022.html
The question is: how to do the sta
> -Original Message-
> This would do it in your example:
>
> > levels(animals$V1) <- c("cat","tiger","dog","fish")
> > table(animals)
>
> cattigerdog fish
> 2 200
But be very wary of
levels(animals$V2)<- c("cat","tiger","dog","fish")
table(animals
> I would like to get the frequency
> (counts) of all the variables in a single column (that is
> easy), but I would also like to return the value 0 for the
> absence of variables defined in another column.
If you use factor() on your columns and include all the animals in the factor
levels, yo
This would do it in your example:
> levels(animals$V1) <- c("cat","tiger","dog","fish")
> table(animals)
cattigerdog fish
2 200
HTH
David
cat tiger dog fish
2 2 0 0
On 11 April 2012 14:21, Daniel Gabrieli wrote:
> animals = as.data
On Apr 11, 2012, at 9:21 AM, Daniel Gabrieli wrote:
Hi,
I hope this is not too trivial, but I've had this recurring problem
and I think there is super easy solution, just not sure what it is.
Please see short example below. I would like to get the frequency
(counts) of all the variables in a
On 11/04/2012 14:08, R. Michael Weylandt wrote:
Well, obviously the interface to .Internal(identical) changed between
2.13.x and 2.15 -- as, ?.Internal says: Only true R wizards should
even consider using this function...
Anyways, find where you use it in your script and then we can help --
ther
Hi,
I hope this is not too trivial, but I've had this recurring problem
and I think there is super easy solution, just not sure what it is.
Please see short example below. I would like to get the frequency
(counts) of all the variables in a single column (that is easy), but I
would also like to r
On Apr 11, 2012, at 9:03 AM, David Winsemius wrote:
On Apr 11, 2012, at 6:28 AM, maxbre wrote:
hi, I just realised I want to go a little further in the control of
the chart
appearance and I would like to have the same number of ticks
displayed in
both axes of all panels
given this code.
Simply pass all = FALSE to merge_all
merge_all(list_of_files, by = "Name", all = FALSE)
Michael
On Wed, Apr 11, 2012 at 1:09 AM, Chintanu wrote:
> Thanks to David and Michael
>
> Michael:
>
> That works, however with a glitch. Each of my 24 files has got two columns:
> "Name", and "Rank/score".
Well, obviously the interface to .Internal(identical) changed between
2.13.x and 2.15 -- as, ?.Internal says: Only true R wizards should
even consider using this function...
Anyways, find where you use it in your script and then we can help --
there's not much we can do without seeing the relevant
Dear all,
I want to explore the nullspace of a matrix S: I currently use the function
Null from the MASS package to get a basis for the null space:
> S = matrix(nrow=3, ncol=5, c(1,0,0,-1,1,1,1,-1,-1,0,-1,0,0,0,-1)); S
> MASS::Null(t(S))
My problem is that I actually need a nonnegative basis for
Here's one:
x <- xts(1:5, Sys.Date() + 1:5)
y <- xts(rep(NA, length(x)), time(x))
or another, less direct but shorter:
y <- cbind(x, NA)[,2]
Michael
On Wed, Apr 11, 2012 at 4:16 AM, jpm miao wrote:
> Hello,
>
> Thank you very much.
>
> How can I create an empty time series with the same
On Apr 11, 2012, at 6:28 AM, maxbre wrote:
hi, I just realised I want to go a little further in the control of
the chart
appearance and I would like to have the same number of ticks
displayed in
both axes of all panels
given this code
xyplot(tv ~ ms | sub_family, data=tm,
#as.tab
hi, I just realised I want to go a little further in the control of the chart
appearance and I would like to have the same number of ticks displayed in
both axes of all panels
given this code
xyplot(tv ~ ms | sub_family, data=tm,
#as.table=TRUE,
aspect="xy",
xlab = expres
Hello, I've just installed R-2.15.0 (I've had R 2.13.2 before it and I've
deleted everything before I started to install new version).
When I've tried to run my script by command "source()" I received this
message:
Error in source("script.R") :
7 arguments passed to .Internal(identical) which r
Dear List,
If I calculate the C-index after 10-fold cross-validation based on Cox
model, which statistic test
should I used to compare two C-statistc values.I means I use
cross-validation to obtain a predictive
probability for each individual, then combine all the probabilities to
calculate th
Hi Pooja,
You must use options command, something like this
options(max.print=5.5E5)
For more information type? ?options
--
Bernardo Rangel Tura, M.D,MPH,Ph.D
National Institute of Cardiology
Brazil
Thanks Bernardo;
cut and paste above line done..
--
View this message in context:
http://r.
On 11/04/2012 12:42, Duncan Murdoch wrote:
On 12-04-11 12:43 AM, Křištof Želechovski wrote:
http://finzi.psych.upenn.edu/R/library/datasets/html/morley.html>
"The classical data of Michaelson and Morley on the speed of light"
Can you provide more information about the data? How were they obtai
On 12-04-11 12:43 AM, Křištof Želechovski wrote:
http://finzi.psych.upenn.edu/R/library/datasets/html/morley.html>
"The classical data of Michaelson and Morley on the speed of light"
Can you provide more information about the data? How were they obtained,
etc.? I do not have the book "Genstat
Dear Michael,
>From ?effect (under "Warnings and Limitations"):
"Binomial generalized linear models cannot have a matrix of successes and
failures on the left-hand side of the model formula; instead specify the
proportion of successes (i.e., successes/(successes + failures)) as the
response, a
On Tue, Apr 10, 2012 at 11:21 PM, jpm miao wrote:
> Hello,
>
> I am writing codes for time series computation but encountering some
> problems
>
> Given the quarterly data from 1983Q1 to 1984Q2
>
> PI1<-ts(c(2.747365190,2.791594762, -0.009953715, -0.015059485,
> -1.190061246, -0.553031799, 0.
> sessionInfo()
R version 2.15.0 (2012-03-30)
Platform: x86_64-apple-darwin9.8.0/x86_64 (64-bit)
locale:
[1] en_US.UTF-8/en_US.UTF-8/en_US.UTF-8/C/en_US.UTF-8/en_US.UTF-8
attached base packages:
[1] grid datasets splines utils stats
[6] graphics grDevices methods base
oth
Dear R users,
I need to fit and then simulate the fitted 2 state Markov-switching VAR
model. My google search shows that there is a package called MSVAR for
problems of this kind. But it seems the package is removed from the CRAN
repository.
May I get a help from someone here please?
Thanks
Mam
Dear Bioc and R List Users,
I am having trouble analysing illumine data generated from BeadScan. I have
.idat files and JPEG images. I realise that i need bead-level summary data to
be able to begin quality control followed by normalization. Is there a way i
can read .idat files for expression a
Hi Petr and Michael,
Thank you very much for the response.
It worked! Now I can use multiple R sessions simultaneously, I just need
to save their workspace separately. And If I want those objects in the
single session, I just load them all in this session, and save the
workspace again togethe
Hi Petr and Michael,
Thank you very much for the response.
It worked! Now I can use multiple R sessions simultaneously, I just need
to save their workspace separately. And If I want those objects in the
single session, I just load them all in this session, and save the
workspace again togethe
Hello,
Thank you very much.
How can I create an empty time series with the same dates as existing
time series?
If x is a ts object, then it would be easy:
y<-ts(NA, start=start(x), end=end(x),frequency=frequency(x))
What can I do if x is a zoo or xts object?
I come up with a c
On 04/11/2012 12:59 AM, Nick Fankhauser wrote:
I've got the strange problem with clock24.plot that only the first data
point (phase = 23.38, size = 0.44) from the phases/sizes numeric vectors
is plotted.
Does anyone have an idea why this could be?
library(plotrix)
phases<- c(23.38, 22.29, 22.71)
On 04/10/2012 11:50 PM, anaraster wrote:
Is there a way to access the numeric results (standard deviation and
correlation) obtained with the taylor.diagram ?
Hi anaraster,
That wouldn't be too difficult, just alter the code to return a list
like this:
return(list(oldpar,R,sd.r,sd.f))
instea
97 matches
Mail list logo
