hi,
I am a novice in R, therefore apologies if my question is too basic, and
thanks for your time.
I used the line below and got the following error:
> getSymbols("DGS1MO", src= "FRED", from="2013-01-01", to=Sys.Date())
Error in read.table(file = file, header = header, sep = sep, quote = quo
Hello,
The Rstudio support forum is at http://support.rstudio.com.
Regards,
Pascal
2013/9/10 Bembi Prima
> Hi All,
>
> Anyone here using RStudio Server?
> I want to ask how can I put .Rprofile that can be accessed by all user.
> So if someone want to add general function, he can just edit one
Hello,
1. You don't say which library you used (I found by myself)
2. The second structure is wrong
3. What specific way is it supposed to follow?
Hereafter is how I understand it:
> plot(layout, t='n')
> text(layout, unique(c(t(data
Regards,
Pascal
2013/9/10 Luca Gaegauf
> Hey everyon
Hi All,
Anyone here using RStudio Server?
I want to ask how can I put .Rprofile that can be accessed by all user.
So if someone want to add general function, he can just edit one .Rprofile
without distribution difficulties.
If I just update .Rprofile in my home directory, other user's .Rprofile
wi
Hello,
Are you sure of the alpha value?
Regards,
Pascal
2013/9/10 Folkes, Michael
> Hello all,
> I've been using binconf (package Hmisc) at a range of alpha values and
> noticed that using the 'Wilson' method when alpha is larger (i.e. narrow
> CI), results in the upper value being smaller th
Hey everyone,
I have a small problem understanding how matrices are accessed.
I created 2 matrices:
data <- structure(c("A", "B", "C", "D", "D", "E", "H", "H", "H", "I",
"F", "G", "F", "F", "G", "G", "F", "G", "I", "J"),
.Dim = c(10L, 2L), .Dimnames = l
Hey everyone,
I have a small problem understanding how matrices are accessed.
I created 2 matrices:
data <- structure(c("A", "B", "C", "D", "D", "E", "H", "H", "H", "I",
"F", "G", "F", "F", "G", "G", "F", "G", "I", "J"),
.Dim = c(10L, 2L), .Dimnames = lis
Please advise
--
Barrett
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
BTW, replying to an old thread here but I ran into this issue today with
the CvM2SL2Test package. The problem is on line 51 of src/XCvMTest.cc in
that package. If you get a stack trace from that memory not mapped
location in your code you'll see its from std::sort, because line 51 is
violating so
Dear R Users:
I would like to use symbolic integration from rSymPy and pass in the actual
function.
So I have
> sympy("x=Symbol('x')")
and if I use
> sympy("integrate(cos(x),x)")
it works fine.
But what I want to do is something like
> y <- "cos(x)"
> sympy("integrate(y,x)")
"cos(x)*x"
Has an
On Sep 9, 2013, at 9:00 AM, Cech, Christian wrote:
> Dear all,
>
> I want to create a line-plot with two lines and some additional scatter-plots.
> However, adding a line to the plot does not work when I specify the
> y-argument in the plot command, while it does work when y is not specified.
>
Hello everyone,
I was trying to fit an arima model to a daily historical data, but, for
some reason, havent been able to.
I basically have 212 observations (from 12/1/2012 to 06/30/2013) containing
the number of transits for a particular vessel.
The following messages are produced by R:
dailytr
Is there a package or a command that does window aggregation like
select
sum(col1) over
(partition by col2, col3 order by col4
rows between unbounded preceding and current row) as sum1
from table1 ;
the above is Netezza syntax, but Postgre has same capability.
Stephen B
[[alternative HT
Hello all,
I've been using binconf (package Hmisc) at a range of alpha values and
noticed that using the 'Wilson' method when alpha is larger (i.e. narrow
CI), results in the upper value being smaller than the lower value. The
'exact' and 'asymptotic' methods give results in the realm I'd expect.
B
Hi David:
Try:
unlist(lapply(sara.data[,4:length(sara.data)],function(x){indx1<-
1:(length(x)/2); indx2<- ((length(x)/2)+1):length(x);
x1<-data.frame(FH=x[indx1],LH=x[indx2],Group1=sara.data$Group[indx1],Group2=sara.data$Group[indx2]);
x2<- x1[!(is.na(x1$FH) | is.na(x1$LH)),]; if(nrow(x2)==0)
>> On Sep 9, 2013, at 9:00 AM, Cech, Christian wrote:
>>
>>> Dear all,
>>>
>>> I want to create a line-plot with two lines and some additional
>>> scatter-plots.
>>> However, adding a line to the plot does not work when I specify the
>>> y-argument in the plot command, while it does work when y i
Hello Arun,
Thanks so much--while I haven't tried it yet, this seems as though it
will be an excellent way to skip the categories (Actb, etc) that have
missing values (NAs).
The second part of my question, which I didn't ask before: instead of
skipping, is there a way to continue with the wilco
Hi,
May be you can try this:
dat1New<- dat1[!(duplicated(dat1$gene)|duplicated(dat1$gene,fromLast=TRUE)),]
dat2<-dat1[duplicated(dat1$gene)|duplicated(dat1$gene,fromLast=TRUE),]
lst1<-split(dat2,dat2$gene)
dat3<-unsplit(lapply(lst1,function(x) {x1<- sum(apply(x[,6:32],2,function(y)
y[1]>=y[2]
On Sep 9, 2013, at 12:27 PM, Cech, Christian wrote:
>>> On Sep 9, 2013, at 9:00 AM, Cech, Christian wrote:
>>>
Dear all,
I want to create a line-plot with two lines and some additional
scatter-plots.
However, adding a line to the plot does not work when I specify the
On Mon, 09 Sep 2013, "Cech, Christian" writes:
> Dear all,
>
> I want to create a line-plot with two lines and some additional scatter-plots.
> However, adding a line to the plot does not work when I specify the
> y-argument in the plot command, while it does work when y is not specified.
>
> I
On 09/09/2013 20:36, Paul Bernal wrote:
Hello everyone,
I was trying to fit an arima model to a daily historical data, but, for
some reason, havent been able to.
I basically have 212 observations (from 12/1/2012 to 06/30/2013) containing
Those dates are not in a standard format: ISO 8601 is p
HI Ahmed,
No problem.
You got the error because one of the IDs had all NAs for AUC24.
Also, In your dataset, there are NAs in BLADDERWALL, BLADDERWALLC, AUC12d,
Cmaxd, in addition to WT, and AUC24. Do you want to follow any particular
methodology for filling the NAs in these columns? or t
HI Ahmed,
No problem.
You got the error because one of the IDs had all NAs for AUC24.
Also, In your dataset, there are NAs in BLADDERWALL, BLADDERWALLC, AUC12d,
Cmaxd, in addition to WT, and AUC24. How do you want to fill those NAs?
res<-unsplit(lapply(split(u3s,u3s$ID),function(x) { x$AUC
Hi R user,
I was trying to put four figures in a row, in which value of Y axis are similar
but different in X axis. Therefore I wanted to put four figures in row.
To remove y axis for second, thrid and forth figures, I used [yaxt="n"]. But, I
saw that tick mark of Y was found to be missing. I d
Revolution Analytics staff write about R every weekday at the Revolutions blog:
http://blog.revolutionanalytics.com
and every month I post a summary of articles from the previous month
of particular interest to readers of r-help.
In case you missed them, here are some articles related to R from t
On 13-09-09 6:14 PM, Paul Bernal wrote:
Dear Mr. Brian,
What is the ISO 8601 standard for dates?
http://lmgtfy.com/?q=ISO+8601
Duncan Murdoch
Excuse my ignorance and best regards,
Paul
El 09/09/2013 16:41, "Prof Brian Ripley" escribió:
On 09/09/2013 20:36, Paul Bernal wrote:
Hello
Thierry: I thought the post asked how to evaluate the arguments in the same
way as the formula ... that is, without changing the paradigm (e.g. by
requiring character arguments). So I think that Harold needs to explicitly
specify in his code where to evaluate semMat, something like:
semMat <- eval
Hi,
I want to obtain the F test associated an ADF test
In tseries, I can obtain t t stat of Dickey fuller, it is posible to obtain
the related F test?
Many thanks in advance
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing
Dear Mr. Brian,
What is the ISO 8601 standard for dates?
Excuse my ignorance and best regards,
Paul
El 09/09/2013 16:41, "Prof Brian Ripley" escribió:
> On 09/09/2013 20:36, Paul Bernal wrote:
>
>> Hello everyone,
>>
>> I was trying to fit an arima model to a daily historical data, but, for
>
Hi,
Please ?dput() your example dataset. Not sure this helps or not.
u3s<- read.table(text="Current-ID visit AUC Wight ID1
101 3 . . 1
101 4 10 13 2
101 5 . . 3
102 3 . . 4
102 4 4 10 5
102 5 . . 6
103 3 . . 7
103 4 6 9 8
103 5 . . 9",sep="",header=TRUE,na.strings=".",check.names=FALSE)
u3d<-
Hi,
You may try:
unlist(lapply(sara.data[,4:length(sara.data)],function(x)
{x1<-tapply(is.na(x),list(sara.data$Groups),FUN=sum); if(x1[1]!=x1[2]) NULL
else
wilcox.test(x~sara.data$Groups,paired=TRUE,alternative="two.sided")$p.value}))
# Bcl2 Ccl5 Cd27 Cd28
#0.1250 0.1875 0.8125 0.8125
In the second example, in the three functions that create the plot:
plot()
lines()
points()
all must have two arguments, the first of which is the x-values (dates),
and the second of which is the y-values associated with the dates.
For example, in your second version, where you have
points
I'm sure I'm missing something really obvious in the "for loop"...
Here is simplified data for 3 patients, we need filling in Na's with same WT
for each patient, AUC halved for visit 3, doubled for visit 5 for the same
patient, based on visit 4
for(i in unique(u3s$ID)){
Hi all,
Imagine I hava e a simulated variable
original<-matrix(rnorm(1),100,100)
If I install tseries and I run adf.test(original[,2],k=0)
the result is:
Augmented Dickey-Fuller Test
data: original[, 2]
Dickey-Fuller = -11.5645, Lag order = 0, p-value = 0.01
alternative hypothesis: stati
Please don't post in HTML. (Read the Posting Guide.)
Consider using the sqldf package.
---
Jeff NewmillerThe . . Go Live...
DCN:Basics: ##.#. ##.#. Live Go...
Dear all,
I am facing a problem with running a script that i have written in R and i hope
you can help me spot which exactly is the problem and how i can solve it.
The error I get is the following:
Error in .jcheck(silent = FALSE) :
Java Exception .jcall(row[[ir]],
"Lorg/apache/poi/ss/us
Hi
I have 6 rather big data sets (between 40 and 80 lines) on transport
data (times, distances and travelers between nodes). They all have a common
index (start-end nodes).
I want to aggregate this data, but for that I have to merge them.
I tried to use "merge" with the result that R (3.0
On 09/09/2013 12:55 PM, Kristi Glover wrote:
Hi R user,
I was trying to put four figures in a row, in which value of Y axis are similar
but different in X axis. Therefore I wanted to put four figures in row.
To remove y axis for second, thrid and forth figures, I used [yaxt="n"]. But, I
saw th
Hi!
I am currently working with a project where I want to plot the regression
line in a plot using ggplot.
The problem occurs when I want to add the second variable, i.e. the z in the
source code:
p = ggplot(data = dat, aes_string(x = "sd", y = "mean", z = "corr"))
p = p + stat_smooth(method = l
Hi!
For example if "data" is the complete dataset with both x and y values:
tempdata = data[complete.cases(data[,1:2]),] # Regression data
model = lm(y~x, data = tempdata) # Linear model
>From this you can calculate the regression value of the missing values.
Hope this helped!
Reg
Your sample code does not run and you didn't show the format of
your inputs. If your 'ttoevP' is a table without duplicate Start/End
pairs that maps Start/End to OEV.T then using match and subscripting
can be quicker than merge. E.g.,
f0 <- function (dataP, ttoevP)
{
encode <- func
On Mon, 9 Sep 2013, bgnumis wrote:
Hi all,
Imagine I hava e a simulated variable
original<-matrix(rnorm(1),100,100)
If I install tseries and I run adf.test(original[,2],k=0)
the result is:
Augmented Dickey-Fuller Test
data: original[, 2]
Dickey-Fuller = -11.5645, Lag order = 0, p-value
Dear all,
I would like to know if there is any R package that can deal with longitudinal
GWAS analysis. Briefly, I have phenotypic data for some individuals (measured
at different time points; not necessarily the same number of measures, so one
individual can have 2 measures over time, while an
Hi John
> 1). Assume now I can calculate these "adjusted" standard deviation from
> sparse PCA, should the percent variation explained by each sparse PC be
> calculated using the sum of all these "adjusted" variance (i.e. square of the
> "adjusted" standard deviation) as the denominator (then t
Dear Harold,
An easy work-around would be to pass the names of the variables as a character
vector.
fm <- lm.eiv(y ~ x1 + x2, dat, ind = c(2,3), semDep = 0, semMat = c("sem1",
"sem2"))
And the change your lm.eiv.fit accordingly.
Or you could have a look at the .() function of the plyr package
I have a function for fitting a type of linear regression and have written
methods for it as shown below. I exclude the main lm.eiv.fit function as it is
large and I don't think necessary for the reproducible example. But, I can
certainly provide if that would be needed.
These methods allow for
Hello.
Relative importance of individual effects' variance: i.e., the standard
quas-demeaning parameter, as in any panel data textbook (e.g. Baltagi,
Wooldridge).
Cheers,
Giovanni
-Messaggio originale-
Da: alfonso.carf...@uniparthenope.it [mailto:alfonso.carf...@uniparthenope.it]
Invi
Hi all,
what indicates the parameter theta in the summary of a random effect
panel model estimated with the plm function?
example:
data("Produc", package = "plm")
zz <- plm(log(gsp) ~ log(pcap) + log(pc) + log(emp) + unemp,
model="random", data = Produc, index = c("state","year"))
summar
Thanks Prof. Ripley,
my.model$family <- binomial() worked
Do you have any suggestions on how to store the models? I used the
models mostly for projections, so being able to use "predict" would be
handy. Refitting the model would be an option but it would not ensure
AFAIK that the resulting mo
hi: my guess is that no one is answering because it's too hard to follow
your code
because it contains so many indices and variables and is without comments.
I don't know where you got that info about the distribution of the
coefficient when doing the ADF test but if you could write the code more
Hi Simon,
thank you for your help. You suggest to use 'alpha_zero' and 'beta_zero' in
the 'L' function. But that's not exactly what I'm trying to do. Maybe it
helps if I show where the 'Likelihood_cov' model developed from. The basic
likelihood look like this:
### Likelihood function ###
Likeliho
On 09/09/2013 10:02, Clara wrote:
Thanks Prof. Ripley,
my.model$family <- binomial() worked
Do you have any suggestions on how to store the models? I used the
models mostly for projections, so being able to use "predict" would be
handy. Refitting the model would be an option but it would not en
i have a data matrix with some x variables complete and some y variables
incomplete. i want to use the simplest regression imputation to fill in the
missing data. (form a regression line with all complete cases and predict
the missing values). is there any package that can do so? if not how should
On 09/09/2013 08:48, Clara wrote:
Hi all,
I have an older, probably about 2 years old glm object, I am not sure
with which version of glm it was produced. I have tried to
Before R 3.0.0, which was a change in major version number.
summary(my.model) but I get an error.
> summary(my.model)
E
I am following your suggestion and move this discussion to the R-devel list.
Best
Simon
On Sep 9, 2013, at 7:58 AM, Prof Brian Ripley wrote:
> On 09/09/2013 02:09, David Winsemius wrote:
>>
>> On Sep 8, 2013, at 8:00 AM, Simon Zehnder wrote:
>>
>>> Dear R-Users and R-Devels,
>>>
>>> I am w
Hi all,
I have an older, probably about 2 years old glm object, I am not sure
with which version of glm it was produced. I have tried to
summary(my.model) but I get an error.
> summary(my.model)
Error in .Call("binomial_dev_resids", y, mu, wt, PACKAGE = "stats") :
"binomial_dev_resids" not
Can anyone help me to see why maximum t-stat quantile has no Negative
values?
¿It is posible lm estimation is not fine?
-- Forwarded message --
From: Jose Narillos de Santos
Date: 2013/9/8
Subject: ADF test
To: r-help
Hi all,
I try to simulate an ADF test rolling on a windo
Hello,
Did you load it?
> require(Hmisc)
Regards,
Pascal
2013/9/9 Rhode Early CHARLES
> Hi,
> I was trying to install.packages("Hmisc"). I followed all the steps but
> nothing appears to have changed on my R.
> my labels still don't appear.
> Rhode
>
> --A Dieu soit la Gloire
>
The time has come to shake the cobwebs off of this analysis. I have
more data now and need to run the same tests, the same way as above.
My question is this--some of the pairs include NAs, and so are gumming
up the works. I'm not sure how to exclude them using the lhs ~ rhs
syntax. Any ideas her
Hi,
I was trying to install.packages("Hmisc"). I followed all the steps but
nothing appears to have changed on my R.
my labels still don't appear.
Rhode
--A Dieu soit la Gloire
[[alternative HTML version deleted]]
__
R-help@r-proj
60 matches
Mail list logo
