On Thu, May 17, 2012 at 09:43:13PM -0400, li li wrote:
Dear all,
I have a function f - function(x,y,c){as.numeric(x*y c) }
I need to solve the value of c so that when I take the integral of
the function f from 0.2 to 0.8 with respect to x and from 0 to 1
with respect to y, the integral
peter dalgaard pdalgd at 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 high
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 the
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
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
On May 18, 2012, at 09:10 , Hans W Borchers wrote:
peter dalgaard pdalgd at 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
Marc Girondot marc_grt at 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
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
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
Hi.
As
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 anilshivku...@gmail.com wrote:
Hi,
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
brianedmunds...@gmail.com wrote:
Dear All,
I'm new at R, but I really just need a couple of
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
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:
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
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:
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
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
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]]
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:
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
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 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
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
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
-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
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
# 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
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
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:
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
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
On
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
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 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
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
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 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
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
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
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
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
Thanks a lot!
Francesca
On 18 May 2012 18:21, arun smartpink...@yahoo.com 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 francesca.panco...@gmail.com
To:
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 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.
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 francesca.panco...@gmail.com
To: r-help@r-project.org
Cc:
Sent: Friday, May 18, 2012 10:59 AM
Subject: [R] How to fix indeces in
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
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
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
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:
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
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,
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)) {
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
oriolebaltim...@gmail.com
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]))
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
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
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
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 -
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
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 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'
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 a
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
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
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
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 started
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
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)
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
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,
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
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 on
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
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 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
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
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 -
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 -
Try this,
bquote(.(L1) * .(L2))
HTH,
b.
On 19 May 2012 16:17, Rolf Turner rolf.tur...@xtra.co.nz 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
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
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 Turnerrolf.tur...@xtra.co.nz wrote:
In the context in which I am actually working it is necessary
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 = something involving xNm)
and get the x axis label to be the Greek
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
83 matches
Mail list logo
