That's the only one I pull up with a quick search, but it seems pretty
good -- what do you need from it?
Best,
Michael
On Fri, May 18, 2012 at 7:16 PM, JWDougherty wrote:
> Can someone point me to an implementation of the "runs up and runs
> down" test, or does such beast exist in R? From web se
Unfortunately the task is slightly harder here -- ugarchsim calls a
network of S4 code that is a little hard to trace
You can get started by typing
getMethod("ugarchsim", "uGARCHfit")
to see the first step which calls a variety of functions with names like
.sgarchsim
You'll have to read a
I had such good luck with my previous question to r-help, (a few minutes
ago) that I thought I would try again with the following query:
Suppose I have
xNm <- "gamma"
I would like to be able to do
plot(1:10,xlab = )
and get the x axis label to be the Greek letter gamma
(rather than
On 19/05/12 16:31, baptiste auguie wrote:
Try this,
bquote(.(L1) * .(L2))
HTH,
b.
Yep! That works! Yew Bewdy!!! Thanks very much.
cheers,
Rolf
On 19 May 2012 16:17, Rolf Turner wrote:
In the context in which I am actually working it is necessary for
me to build parts o
The smooth.sspline function in sspline v0.1-5 is giving different results
depending whether I run it in 32-bit or 64-bit platform. Correct results
(TRUE result below) occur on the 32-bit platform. I have run it on
Windows 7 x32 +R2.15.0 (works) and R2.14.2 (works), Windows 7 x64 +R2.15.0
(d
Try this,
bquote(.(L1) * .(L2))
HTH,
b.
On 19 May 2012 16:17, Rolf Turner wrote:
>
> In the context in which I am actually working it is necessary for
> me to build parts of a text string to place on a plot in two separate
> stages.
>
> The following toy example illustrates what I am trying to
In the context in which I am actually working it is necessary for
me to build parts of a text string to place on a plot in two separate
stages.
The following toy example illustrates what I am trying to do:
L1 <- substitute(list(paste("Bias in ", gamma,", ")))
L2 <- substitute(list(beta==b0),lis
Thanks Petr and David.
Sorry David if I was not clear enough. The last comment line highlights the
end objective (hopefully more clearly this time).
Petr: I kind of see your line of thought, but still cannot see how it works
on a specific example like this one.
set.seed(1)
dframe <- matrix(
Thanks Rich,
I was looking for something I could run *inside* of R, as part of my R script.
How are you using ncurses?
--
Noah Silverman
UCLA Department of Statistics
8117 Math Sciences Building
Los Angeles, CA 90095
On May 18, 2012, at 5:14 PM, Rich Shepard wrote:
> On Fri, 18 May 2012, Noah
On Fri, 18 May 2012, Noah Silverman wrote:
I need to design a fairly simple front-end for someone to use an R script
system that I've built. My thought was to just use the text based menus
available in the base R package, perhaps in some kind of loop.
Noah,
ncurses will do all you want in
Hello,
I need to design a fairly simple front-end for someone to use an R script
system that I've built. My thought was to just use the text based menus
available in the base R package, perhaps in some kind of loop.
How have other people done this? Any "best practices" that you can recommend?
On May 18, 2012, at 1:44 PM, Robbie Edwards wrote:
Hi all,
I'm trying to model some data where the y is defined by
y = summation[1 to 50] B1 * x + B2 * x^2 + B3 * x^3
Hopefully that reads clearly for email.
cumsum( rowSums( cbind(B1 * x, B2 * x^2, B3 * x^3)))
Anyway, if it wasn't for
Can someone point me to an implementation of the "runs up and runs
down" test, or does such beast exist in R? From web searches the
"runs up and runs down" test is commonly used for testing pseudo-random
number generators and in simulations. John C. Davis describes its use
in geology in his book o
Following up on Rolf's post:
1) cumulative summation (cumsum) maybe?
2) In fact, you should probably **not** fit the non-summation version
as you have stated. See ?poly.
I would guess that context is important here. Based on (my
interpretation) of the rather strange nature of your request, I
sus
I would like to nominate the following quote from a posting to R-sig-Geo
by Roger Bivand as a fortune:
Some OSX users know that OSX is really Unix [...]; others think that
OSX is cooler
Windows, and they have ontological problems with non-Apple phenomena
and constructs.
cheers,
Thanks for assisting here Paul.
It turns out there is a bug in this release of Rattle. It is fixed and
2.6.19 will include the fix.
In the mean time you can download an updated 2.6.18 from
http://www.togaware.com.au/repository/
Regards,
Graham
On 19 May 2012 01:14, Paul Miller wrote:
> Hi Avid
On 19/05/12 05:44, Robbie Edwards wrote:
Hi all,
I'm trying to model some data where the y is defined by
y = summation[1 to 50] B1 * x + B2 * x^2 + B3 * x^3
Hopefully that reads clearly for email.
Anyway, if it wasn't for the summation, I know I would do it like this
lm(y ~ x + x2 + x3)
Whe
Your issue is that nls returns "singular gradient", but I believe the real
trouble is that
if or ifelse are not in the derivatives table. My own (experimental, but in
R-forge)
package nlmrt has nlxb to mirror nls but be much more aggressive in finding a
solution,
and it gives an error msg that i
On May 18, 2012, at 2:26 PM, Giggles wrote:
> I am a newbie and can't figure out how to recode a numeric value. In
> my data (pharm311), I have a column called "explain" and I need to
> find all the 6's and change it to NA (blank). Could someone help?
>
> I'm sorry if this is too basic, I starte
On May 18, 2012, at 3:26 PM, Giggles wrote:
I am a newbie and can't figure out how to recode a numeric value. In
my data (pharm311), I have a column called "explain" and I need to
find all the 6's and change it to NA (blank). Could someone help?
is.na(pharm311$explain) <- pharm311$explain==6
Hi all,
First:
I have a small line of code I'm applying to a variable which will be
placed in a matrix table for latex output of accuracy measures:
acc.aarima <- signif(accuracy(forecast(auto.arima(tix_ts,
stepwise=FALSE), h=365)), digits=3).
The time series referred to is univariate (daily c
I am a newbie and can't figure out how to recode a numeric value. In
my data (pharm311), I have a column called "explain" and I need to
find all the 6's and change it to NA (blank). Could someone help?
I'm sorry if this is too basic, I started messing with R this week and
got stuck with this probl
I am looking into using the Matching package for some multi-level analysis
with multiply imputed datasets, and so far I have been unable to get a
handle on how to appropriately deal with these issues when using genmatch()
in the Matching library. Here are my basic questions:
-Is the 'weights' opt
Hello,
There are two simple ways,
1. Type the function name without parenthesis at an R prompt;
2. Download the source from CRAN or R-Forge, R is open source.
Hope this helps,
Rui Barradas
jaimie villanueva wrote
>
> hi
>
> someone can show me how can i get the source code of a function. Is
You could use parse(text=a)[[1]]:
> numericDeriv(parse(text=a)[[1]], c(t(l),recursive=TRUE))
[1] 3
attr(,"gradient")
[,1] [,2] [,3] [,4] [,5] [,6]
[1,]103000
or construct your expression directly:
> lnames <- array(lapply(l, as.name), dim=dim(l))
> aa <-
I missed a couple line of codes in the previous e-mail. Here is whole code
again:
load <- matrix(c(3,0,1,4,1,3),nrow=3,ncol=2,byrow=TRUE)
l <- matrix(nrow=nrow(load),ncol=ncol(load))
for(i in 1:nrow(load)) {
for(j in 1:ncol(load)) { l[i,j]=paste("l",i,j,sep="")}}
for(i in 1:nrow(load)){
for(j in
Hi,
I am stuck on something for a couple days, I am almost about to give up.
This looks simple, but I can't figure out. I hope I can get some help here.
I am trying to do some symbolic and numerical derivations. Let me explain
the problem. Let's say, I have a matrix as follows:
> load <- matrix(
On May 18, 2012, at 10:52 AM, Duncan Murdoch wrote:
On 12-05-18 5:15 AM, gholi bahrami wrote:
Dear Sir/Madam
I am a BSc student studying statistics in Ferdowsi University of
Mashhad,Iran.
Translating R messages to Farsi is my thesis and I have already
downloaded *poEdit* to translate thos
Hi all,
I am running 64-bit R 2.15.0 on windows 7. I am trying to use read.delim
to read from a file that has 2-byte unicode (CJK) characters.
Here is an example of the data (it is tab-delimited if that gets messed up):
HITId HITTypeId Title
2Q69Z6KW4ZMAGKKFRT6Q4ONO6MJF68 2LVJ1LY58B72OP36GNBHH16
On May 18, 2012, at 12:56 PM, Matthew Ouellette wrote:
Dear R help list,
I am very new to R and I apologize in advance if this has been
answered
before. I have done my best to google/R search what I need but no
luck.
Here is what I am attempting:
I have hundreds of .csv files that I nee
Hello,
Try the following.
# Make some data
alldata <- list(matrix(rnorm(12), ncol=3), matrix(sample(100), ncol=10))
(alldata <- lapply(alldata, function(x){colnames(x) <- c("Name",
LETTERS[2:ncol(x)]); x}))
# This does the trick
all.order <- lapply(alldata, function(x) order(x[, "Name"]))
lappl
Have you read "An Introduction to R" ? This is a pretty basic
exercise that you should be able to do by yourself if you have. If you
haven't, you need to (or some other basic R tutorial -- there are many
around).
-- Bert
On Fri, May 18, 2012 at 10:41 AM, Adrian Johnson
wrote:
> Hi:
>
> I am loo
You appear to have a good start.
If you type
alldata[[1]]
do you get what you expect for the first file?
This is not tested, but I would start with something like this:
sorteddata <- lapply(alldata, function(df) df[order(df$Name),] )
## then this will overwrite
for (id in seq(filenames)) {
Hi all,
I'm trying to model some data where the y is defined by
y = summation[1 to 50] B1 * x + B2 * x^2 + B3 * x^3
Hopefully that reads clearly for email.
Anyway, if it wasn't for the summation, I know I would do it like this
lm(y ~ x + x2 + x3)
Where x2 and x3 are x^2 and x^3.
However, sin
Hi:
I am looking for some help in making two boxplots next to each other.
I have a data like this:
N1 T1 N2 T2 N3 T3 N4 T4 ... Nn Tn
7 8.2 4 5 8104 5 . 10 11
I want to have box plot for all Normal samples (N1,N2,N3,N4Nn)
and another box plot
that´s a success. "Das ist ein Erfolg". :) Maybe you had that in your german
lessons, too ;)
So if you are still interested in learning our language, feel free to ask
me.
I really learned a lot in this forum, so i am always happy to be able to
give something back
--
View this message in context
hi
someone can show me how can i get the source code of a function. Is a S4
class or Method. (I'm not an expert in R environment)
Exactly, Function "ugarchsim" from library (rugarch).
I need to know (in detailed ) how the variance and mean ecuation of a
arma/garch process are calculated.
With oth
Dear R help list,
I am very new to R and I apologize in advance if this has been answered
before. I have done my best to google/R search what I need but no luck.
Here is what I am attempting:
I have hundreds of .csv files that I need to sort based on a single column
of alphanumeric data. All o
Dear David,
Thank you for your time and consideration,
My data set is Sea Surface temperature (sst) of 50 year monthly global
data (1x1 degree resolution ). My goal is to calculate the EOF1 , EOF2 &
EOF3 at each Lat-Lon location. How to use 'princomp' command usefully to
calculate these
Hi Francesca,
> for(i in 1:length(x1<-c(100,1000,1))){
j<-x1[i]
x1[i]<-mean(j)
}
> x1
[1] 100 1000 1
A.K.
- Original Message -
From: Francesca
To: r-help@r-project.org
Cc:
Sent: Friday, May 18, 2012 10:59 AM
Subject: [R] How to fix indeces in a loop
Dear Contributor
Thanks Milan, it is running now. It seems part of the problem, as you
suggested were the packages. It seems that although I just installed Rweka,
Snowball and the like they were out of date. So updataing fixed
stemDocument. As for removeWords, that began working once I cut my data in
half. Appa
Hi Uwe:
Thanks a lot for your reply.
I contemplated the use of the predict() function but realized that it could
not return what I needed (i.e., the terminal nodes and not the classes). I
went back to the post of Professor Ripley
(http://tolstoy.newcastle.edu.au/R/e4/help/08/07/17702.html) and
Thanks a lot!
Francesca
On 18 May 2012 18:21, arun wrote:
> Hi Francesca,
>
>> for(i in 1:length(x1<-c(100,1000,1))){
> j<-x1[i]
> x1[i]<-mean(j)
> }
>
>> x1
> [1] 100 1000 1
>
>
>
> A.K.
>
>
>
> - Original Message -
> From: Francesca
> To: r-help@r-project.org
> Cc:
> Sent
Check the posting guidelines and give us a small reproducible example using
dput(). It is here
http://www.R-project.org/posting-guide.html
You say you want "PCs of a spatial data set (single variable)", but you must
mean something else. It sounds like your variables are highly correlated
with one
Hi Francesca,
Here's one approach:
loopvalues <- c(100,1000,1)
results <- numeric(length(loopvalues))
for(loopindex in 1:length(loopvalues)) {
i <- loopvalues[loopindex]
results[loopindex] <- mean(i)
# assuming your intent is actually something other
# than ta
Hi Avideh,
Sorry to hear you're still having a problem. Tried usinging this two options on
my 64-bit and 32-bit versions of R 2.15.0. Everything works fine. So I've at
least been able to determine that there's nothing wrong with Rattle or with the
combination of Ratttle and either the 64-bit an
If you are unfamiliar with Windows command-line, be warned that it is only
vaguely similar to sh. You might find that reading the documentation on
developing R extensions has hints to what you need to do.
In general, you have to modify environment variables manually often to make
command line t
On May 18, 2012, at 5:41 AM, Nick Riches wrote:
Hi
I used the package Natural Language Toolkit about a year ago to parse
strings into parts of speech. However, I cannot find it on CRAN.
Are you perhaps confused about which 'unnatural language' it was
written in?
http://nltk.org/
Has i
Dear Contributors,
I have an easy question for you which is puzzling me instead.
I am running loops similar to the following:
for (i in c(100,1000,1)){
print((mean(i)))
#var<-var(rnorm(i,0,1))
}
This is what I obtain:
[1] 100
[1] 1000
[1] 1
In this case I ask the software to print out
On 12-05-18 5:15 AM, gholi bahrami wrote:
Dear Sir/Madam
I am a BSc student studying statistics in Ferdowsi University of Mashhad,Iran.
Translating R messages to Farsi is my thesis and I have already downloaded
*poEdit* to translate those messages.However, the created .mo files for R.pot
and R
Hi everybody!
Does anyone know how to obtain a inter-item-correlation-table
(with p-values or significance-levels)? (as SPSS does, either spearman or
pearson)
Repeatedly using cor.test() is pretty exhausting as the table size
increases...
Thank you!
David
[[alternative HTML version del
Sorry about that.
Basically I have two fixed factors (experimental group and sex) and a random
factor
(brood). I need the random factor to control for pseudoreplication because I
have siblings,
which are pseudoreplicates because share genes/environment.
however, it is not a nested design. Then I
Hi, after re-reading I think that I misunderstood your question. You
don't provide many details, but I suppose that the "brood" effect is
nested within the fixed effects, so you don't mean a multivariate
approach for a split-plot or a repeated-measures design, but the
analysis of a multivariate mix
Thank you Robert,Petr and Bert.
I know my question is generic one, which does specify what is range to
values which comes as "unusual" value.
I start with R a week ago and trying to learn it.
Can somebody tell me how to set the range (say >700) for the outlier
value?
Thanks & Regards,
Prakash
Hello everyine,
I am trying to run R in Batch mode in windows, but can't find any help
in google. In particular, the only information I get is that I need to
have an R.exe file, than I can't find in the bin folder. Can you
please provide me with a link or refence on how to run R in batch mode
in w
You're looking in the wrong area: capacity planning is fitting observations
to a hyperbola, not a straight line (:-))
See the R work by Neil Gunther, instead: http://www.perfdynamics.com/
--dave
--
View this message in context:
http://r.789695.n4.nabble.com/Extrapolation-for-capacity-planning-t
On May 18, 2012, at 6:37 AM, Axel Urbiz wrote:
Would I be able to accomplish the same if x.sample was created from x
instead of x.sorted. The problem is that in my real problem, I have
to sort
with respect to many variables and thus keep the sample indexes
consistent
across variables. So I
# Hi Manish, try this
# locate the y axis central positions
barObj <- plot(factor(mtcars$gear, labels = paste(3:5, "Gears")), horiz = TRUE)
# create the plot
plot(factor(mtcars$gear, labels = paste(3:5, "Gears")), horiz = TRUE)
# add axes
axis(side = 2, at = c(barObj), labels = FALSE)
Chris C
Hi Avideh,
I have Rattle installed on the 64-bit version of R 2.15.0 and this seems to
work OK. I had a dataset of my own loaded into Rattle when I saw your email.
Tried creating the plot using these data but got a warning saying I should
limit the number of variables to 40. So I changed to the
-Original Message-
From: Manish Gupta
Sent: Friday, May 18, 2012 2:52 AM
To: r-help@r-project.org
Subject: [R] How to create axis y axis for horizontal bar plot?
Hi,
i am working on bar plot but i need to generate y axis for horizontal bar
plot. In the attached diagram x-axis is there
Dear all,
I am trying to find the PCs of a spatial data set (single
variable). I want to calculate the PCs at each Lat-Lon location.
The* 'princomp'* command gives the approximate standardized data,
(i.e* pca$scores*), stranded deviation ..etc. I tried*
'pca$loadings'*also, but
No it's not a bug in R.
Indexing can only occur by integers so there's an automatic conversion
to integer inside of `[`. Your "insert" variable is actually slightly
less than 8 (note as.integer(insert) == 7) by some amount that's too
small to print, but the computer still picks up on.
More concre
Dear R users
I used glmnet generating a regression model, now I need to convert it to pmml
format, but I noticed pmml r package doesn't support glmnet object, has anyone
find a way solving this problem? I was thinking convert glmnet object to glm
object, has anyone tried it?
Many thanks
Yan
Well, it's not quite the dput(head(X, 25)) I asked for, but it
worksthank you :-)
Do take a look at the examples for aggregate -- then the right syntax
is pretty clear
z <- aggregate(x, timeSequence(from = start(x), to = end(x), by = "day"), sum)
but there seems to be some little issue that
Dear Sir/Madam
I am a BSc student studying statistics in Ferdowsi University of Mashhad,Iran.
Translating R messages to Farsi is my thesis and I have already downloaded
*poEdit* to translate those messages.However, the created .mo files for R.pot
and RGui.pot are not useful at all.
I would like
Hi,
i am working on bar plot but i need to generate y axis for horizontal bar
plot. In the attached diagram x-axis is there with scale 0 to 12 but i need
y axis. How can i implement it?
http://r.789695.n4.nabble.com/file/n4630478/barplot2.jpg
Thanks
--
View this message in context:
http://
Hi
I have one similarity matrix which is not 0-1 scaled.
Higher values mean higher similarity in the matrix.
I would like bottom-up plot. I know how to plot root up plot.
Does anyone know?
Br,
Luffy Liu
[[alternative HTML version deleted]]
Dear all,
I am trying to find the PCs of a spatial data set (single
variable). I want to calculate the PCs at each Lat-Lon location.
The* 'princomp'* command gives the approximate standardized data
(i.e* pca$scores*), stranded deviation ..etc. I tried*
'pca$loadings'*also, but it
Dear community,
I have a little problem filling in a vector with loop output.
I think the problem is the following: I am caluclating the index which
indicates where the loop output should be placed at by a formula. But when I
want to fill in or read out some specific index places this doesn't work
Dear R users
I generated a model using glmnet, I need to convert it to pmml, but R pmml
package doesn't support glmnet, has anyone come across similar problem? any
idea to solve it?
Many thanks
YAn
--
View this message in context:
http://r.789695.n4.nabble.com/glmnet-object-to-pmml-tp4630493.h
On Fri, May 18, 2012 at 06:37:03AM -0400, Axel Urbiz wrote:
> Would I be able to accomplish the same if x.sample was created from x
> instead of x.sorted. The problem is that in my real problem, I have to sort
> with respect to many variables and thus keep the sample indexes consistent
> across var
I suppose you mean this:
http://www.r-project.org/doc/Rnews/Rnews_2007-2.pdf
Another approach, less general but in my opinion easier to use (and
equivalent in many situations), is provided by Anova() in package car.
See:
http://socserv.mcmaster.ca/jfox/Books/Companion/appendix/Appendix-Multivariat
Would I be able to accomplish the same if x.sample was created from x
instead of x.sorted. The problem is that in my real problem, I have to sort
with respect to many variables and thus keep the sample indexes consistent
across variables. So I need to first take the sample and then sort it
with res
You're leaving out a step -- take a look at the examples for viewFin
and note that you have to download the financial data first with
getFin()
Michael
On Thu, May 17, 2012 at 9:53 PM, Brian Edmundson
wrote:
> Dear All,
>
> I'm new at R, but I really just need a couple of things. The first thing
Well, the easiest way would be to use SAS ;-)
In R take a look at the ar() function. You can view documentation by
typing ?ar at the console and see some worked examples by running
example(ar)
Hope this helps,
Michael
On Fri, May 18, 2012 at 12:11 AM, anil wrote:
> Hi,
>
> I want to find out, h
Hi
I used the package Natural Language Toolkit about a year ago to parse
strings into parts of speech. However, I cannot find it on CRAN.
Has it been withdrawn?
Thanks
Nick
--
Lecturer in Speech and Language Pathology
Room 1.9
King George VI Building
Queen Victoria Road
University of Newcastl
On Thu, May 17, 2012 at 06:14:37PM -0400, Nathan Stephens wrote:
> I have a very simple maximization problem where I'm solving for the vector
> x:
>
> objective function:
> w'x = value to maximize
>
> box constraints (for all elements of w):
> low < x < high
>
> equality constraint:
> sum(x) = 1
Marc Girondot yahoo.fr> writes:
>
> Le 18/05/12 00:14, Nathan Stephens a écrit :
> > I have a very simple maximization problem where I'm solving for the vector
> > But I get inconsistent results depending on what starting values I. I've
> > tried various packages but none seem to bee the very sol
On May 18, 2012, at 09:10 , Hans W Borchers wrote:
> peter dalgaard gmail.com> writes:
>>
>> On May 18, 2012, at 00:14 , Nathan Stephens wrote:
>>
>>> I have a very simple maximization problem where I'm solving for the vector
>>> x:
>>>
>>> objective function:
>>> w'x = value to maximize
>>>
Is there any way to analyse high frequency data in R like
cleaning,manipulation and volatility etc.I know there are packages like
RTAQ and Realized for analysing high frequency data but they are only valid
for NYSE stocks and have well defined data format.
Please help.
[[alternative HTML
Le 18/05/12 00:14, Nathan Stephens a écrit :
I have a very simple maximization problem where I'm solving for the vector
x:
objective function:
w'x = value to maximize
box constraints (for all elements of w):
low< x< high
equality constraint:
sum(x) = 1
But I get inconsistent results dependi
On 17/05/2012 20:35, Sophie Baillargeon wrote:
Hi,
When I run the following code :
Y<- c(rep(0,35),1,2,0,6,8,16,43)
cst<- log(choose(42, 42:1))
beta<- 42:1
tau<- (beta^2)/2
fit<- glm(formula = Y ~ offset(cst) + beta + tau, family = poisson)
fit
fit$converged
glm prints a warning saying that th
peter dalgaard gmail.com> writes:
>
> On May 18, 2012, at 00:14 , Nathan Stephens wrote:
>
> > I have a very simple maximization problem where I'm solving for the vector
> > x:
> >
> > objective function:
> > w'x = value to maximize
> >
> > box constraints (for all elements of w):
> > low < x
83 matches
Mail list logo