On 28 Oct 2013, at 13:07 , kmmoon100 wrote:
> Hello everyone,
>
> This is Kangmin.
>
> I am trying to produce shape and scale of my wind data. My data is based on
> wind speed frequency with 1km/hr increment. data is described below.
>
> Windspeed (km/h)Frequency
> 1 351
> 2 147
>
On 29-10-2013, at 00:35, kmmoon100 wrote:
> Hi Berend,
>
> Thank you for your reply.
> How can I use dput function for this type of data?
> I looked up the description of the function but I still can't understand how
> to use it for solving my error.
>
You don't use dput() to solve your error
Hi Berend,
Thank you for your reply.
How can I use dput function for this type of data?
I looked up the description of the function but I still can't understand how
to use it for solving my error.
Regards,
Kangmin.
--
View this message in context:
http://r.789695.n4.nabble.com/Optimization-f
I have a class Details that contains information needed by FirstSet to do some
calculations then a super class that returns Details and FirstSet. The problem
seems to be in FirstSet where I use the function getAnumber(id).
setClass("Details",
representation(
ID = "charact
Hello,
Thanks a lot for the replies.
I tried to use the "optim" function in R, and chose the "L-BFGS-B" method
of optimization. Now the result is very sensitive to the starting values. I
use the same function, but I give the line of code I use for "optim".
optim(0.25, f, method= "L-BFGS-B", low
> One can do better yet by using the logspace_add(logx, logy) function
> available in the C API to R, but there is currently not an R-language
> interface
> to that (or logspace_subtract(logx,logy)). They calculate log(exp(logx) +-
> exp(logy))
> with minimal roundoff error.
>
> You could proba
On Mon, 28 Oct 2013 15:03:36 -0700
jwd wrote:
> On Mon, 28 Oct 2013 13:40:51 +
> "Ortiz, John" wrote:
...
> >
> It is not clear what you want. Why is there an 'x = "log"' term in
> plot()? If you want to plot log(x) try using:
Got that backward typing, should be log = "x"
__
On Mon, 28 Oct 2013 13:40:51 +
"Ortiz, John" wrote:
>
> Dear list users,
>
>
> I'm doing a plot integrating Grid output with Base Graphics output
> (gridBase, Murrell 2012).
>
> My goal is to produce a xy plot where each point is represented by a
> pie. I could get it using the attach co
Note that f(x) returns exactly zero for all x above 1, hence its estimated
derivative is 0 everywhere in that region. You are asking optimize to find
the non-zero part of a function which is 0 in more than 80% of its domain.
You can put a trace on f() to see where optimize() looks:
> trace(f,
Your question still makes no sense at all to me. You provide no example
code
and no data. You do not specify, in any comprehensible way, what model you
are trying to fit.
And as Jeff Newmiller said, if you are concerned with "Revo R" why are
you posting
to "r-help"?
Please get your act to
This could be described as a bug, perhaps. Or it could be described as
an indication that numerical optimization is inevitably tricky. Notice that
if you narrow down your search interval from [0,5] to [0,0.5] you get the
right answer:
> optimize(f, c(0, 0.5), maximum= TRUE,tol=1e-10)
$maximum
[
Hello,
This seems like a bug, removing 'maximum = TRUE' has no effect except
that the returned value says it's a minimum, not a maximum.
Rui Barradas
Em 28-10-2013 17:00, Shantanu MULLICK escreveu:
Hello Everyone,
I want to perform a 1-D optimization by using the optimize() function. I
want
Hello,
Instead of string manipulation you should think of POSIXt objects. And
maybe ?seq.POSIXt will help.
first <- strptime("2011-10-12 10:59:00","%Y-%m-%d %H:%M:%S")
last <- strptime("2011-10-12 11:13:00","%Y-%m-%d %H:%M:%S")
seq(first, last, by = "5 min")
Hope this helps,
Rui Barradas
Please avail yourself of search tools before posting here, and reference why
they did not answer your question.
RSiteSearch( "chess" )
---
Jeff NewmillerThe . . Go Live...
DCN:
Hello again, first off thank you for your suggestion Mr. Rigby, I'll take
a look at the GAMLSS package.
I have a (slightly) related followup question regarding the 'fitdistr'
function. I was examining my data, a sample of which is attached, using
this function and I am confused about the interpre
On 10/29/2013 02:30 AM, Ahmed Attia wrote:
Hi
I have a question about drawing a linear line in x, y plot. I usually
use the following code, but for this time the x values are to small
(-0.08 to -0.02)
I wrote the following code, but r does not draw the line. However, it
does not give an error w
On 10/29/2013 01:51 AM, Johannes Radinger wrote:
Hi,
I'd like to follow up an older discussion on adding significance stars to
boxplots.
(
http://r.789695.n4.nabble.com/indicating-significant-differences-in-boxplots-td862335.html
)
As suggested by Jim Lemon, there is following solution for a sin
Thanks a lot Simon, that's useful.
I will take a look at this process-pinning things.
Arnaud
2013/10/28 Simon Zehnder
> First,
>
> use only the number of cores as a number of thread - i.e. I would not use
> hyper threading, etc.. Each core has its own caches and it is always
> fortunate if a p
I'm trying to run bGLMM on a large dataset and it is failing after trying to
create huge vectors with the error:
Error: cannot allocate vector of size 553.7 Gb
The dataset that I'm using is fairly large, around 1.5 million observations.
The command is as follows:
boost1 <- bGLMM(A~ B + C, rnd = l
Hi,
I'm trying to figure out a way to make the output of a call to system()
invisible on the R mac GUI. Is there a mac equivalent for the calls
invisible=TRUE, show.output.on.console=FALSE, intern=TRUE for the mac
environment? Thanks.
~John
__
R-help
Hello Everyone,
I want to perform a 1-D optimization by using the optimize() function. I
want to find the maximum value of a "logistic" function. The optimize()
function gives the wrong result.
My code:
f= function (k) {
T_s = 20
result = (2- 2/(1+ exp(-2*T_s*k)))
return(result)
}
optimize(f, c(0
Have there been attempts to build a chess engine in R? If so, is it
available anywhere?
I recognize that R really isnt the right language for a chess engine, but
I was more curious if it had been attempted.
Thank you.
Rob Bernard
--
--
Robert N. Bernard
minus...@alumni.princeton.
Hi All,
This is my first time to seek help from the R-forum (and also conduct a formal
data-mining analysis). I searched the archive a bit but didn't get responses
that could totally address my question. Any comments would be highly
appreciated.
I am using the rpart function to analyze factors
Dear all,
in my code I have written the following list
TimeFramesShort <-list(c(strptime("2011-10-12 10:59:00","%Y-%m-%d
%H:%M:%S"),strptime("2011-10-13 11:02:00","%Y-%m-%d %H:%M:%S")) #
rest items were truncated
I was wondering if could somehow take the first element of the lis
Hi Arun,
Sorry for late reply.
I would like to thank you for your pointer. This is what I wanted. Thanks
for your help.
Thanks and regards,
On Sun, Oct 27, 2013 at 10:33 PM, arun wrote:
> Hi,
> DF$Col2
> # [1] a e b b a b c e b a c c e e a c d c c e
> #Levels: a b c d e
>
>
> "b" is not foun
Hi,
You may try:
x <- 1:5
set.seed(49)
mat1 <- do.call(rbind,lapply(1:1000,function(y) sample(x,3)))
#or
mat2 <- matrix(0,ncol=3,nrow=1000)
set.seed(49)
for(i in seq_len(nrow(mat2))) mat2[i,] <- sample(x,3)
all.equal(mat1,mat2)
#[1] TRUE
#or
set.seed(49)
mat3 <- t(replicate(1000,sample(x,3))
On 28/10/2013 16:19, Arnaud Mosnier wrote:
Hi all,
I am quite new in the world of parallelization and I wonder if there is a
way to increase the speed of creation of a parallel socket cluster. The
time spend to include threads increase exponentially with the number of
It increases linearly in
First,
use only the number of cores as a number of thread - i.e. I would not use hyper
threading, etc.. Each core has its own caches and it is always fortunate if a
process has enough memory; hyper threads use all the same cache on the core
they are running on. detectCores() gives me for exampl
Thanks Simon,
I already read the parallel vignette but I did not found what I wanted.
May be you can be more specific on a part of the document that can provide
me hints !
Arnaud
2013/10/28 Simon Zehnder
> See library(help = "parallel)
>
>
> On 28 Oct 2013, at 17:19, Arnaud Mosnier wrote:
>
On Oct 28, 2013, at 6:40 AM, Ortiz, John wrote:
Dear list users,
I'm doing a plot integrating Grid output with Base Graphics output
(gridBase, Murrell 2012).
My goal is to produce a xy plot where each point is represented by a
pie. I could get it using the attach code,
but now I wan
See library(help = "parallel”)
On 28 Oct 2013, at 17:19, Arnaud Mosnier wrote:
> Hi all,
>
> I am quite new in the world of parallelization and I wonder if there is a
> way to increase the speed of creation of a parallel socket cluster. The
> time spend to include threads increase exponentiall
On Oct 28, 2013, at 8:30 AM, Ahmed Attia wrote:
Hi
I have a question about drawing a linear line in x, y plot. I usually
use the following code, but for this time the x values are to small
(-0.08 to -0.02)
That is not the problem.
I wrote the following code, but r does not draw the line.
Dear all,
I am trying to implement a function which removes aliased terms from a model.
The challenge I am facing is that with "alias" I get the aliased coefficients
of the model, which I have to translate into the terms from the model formula.
What I have tried so far:
--8<---
Hi all,
I am quite new in the world of parallelization and I wonder if there is a
way to increase the speed of creation of a parallel socket cluster. The
time spend to include threads increase exponentially with the number of
thread considered and I use of computer with two 8 cores CPU and thus
sh
On 28-10-2013, at 17:16, Baro wrote:
> thank you so much. I will write my questions as you mentioned
>
>
Final advice.
Please reply to the list only. Everybody can then see your reply and respond
if necessary.
Forwarded to the list.
Berend
__
On 28-10-2013, at 17:01, Baro wrote:
> Hi experts,
>
> I want to user haar wavelet transform. If I am using this simple command,
> it works nice:
>
> k<-c(1,2,3,4,5,6,7,8)
> ywd<-wd(k,filter.number=1,family="DaubExPhase")
>
> but if the K is a list like this:
>
> *[[1]]*
> *[1] 401*
> *
> *
A few problems ...
This statement doesn't make sense.
seq(-0.08, -0.02, len = -0.02)
Perhaps you meant
seq(-0.08, -0.02, by = 0.02)
The xlim= and ylim= are arguments to higher level plot functions, like
plot(), and won't work for functions lines() or abline().
Are you trying to limit the range o
Hi experts,
I want to user haar wavelet transform. If I am using this simple command,
it works nice:
k<-c(1,2,3,4,5,6,7,8)
ywd<-wd(k,filter.number=1,family="DaubExPhase")
but if the K is a list like this:
*[[1]]*
*[1] 401*
*
*
*[[2]]*
*[1] 481*
*
*
*[[3]]*
*[1] 480*
*
*
*[[4]]*
*[1] 482*
*
*
*[
Not certain about .por but this works with ordinary SPSS files:
require(Hmisc)
dat <- spss.get(...) # gets variable labels, etc.
contents(dat)
html(contents(dat), ...)
The last command produces a hyperlinked data dictionary, e.g., for each
variable the number of levels is given and you click o
-- Forwarded message --
From:
Date: Mon, Oct 28, 2013 at 10:53 AM
Subject: Undeliverable mail
To: hill0...@umn.edu
Cc: postmas...@vs-m.tc.umn.edu
Your message was not delivered to the following recipients:
r-h...@r-proj.org: No such host
Original-Recipient: rfc82
You still don't seem to get the hint that support that is SPECIFIC to
Revolution R is off- topic on this mailing list. Since you have identified your
subject so clearly, we can only safely respond that we cannot answer you.
If you change your mind and decide to ask what can be done in R (which s
Ben,
It helps us help you if you provide a simple example with code that anyone
can run. When I created a simplified version of you situation, the NAs in
the response variable yielded NAs in the residuals just as the help for
na.exclude indicates.
?na.exclude
Can you replicate the problem (and
Hi experts,
I want to user haar wavelet transform. If I am using this simple command,
it works nice:
k<-c(1,2,3,4,5,6,7,8)
ywd<-wd(k,filter.number=1,family="DaubExPhase")
but if the K is a list like this:
*[[1]]*
*[1] 401*
*
*
*[[2]]*
*[1] 481*
*
*
*[[3]]*
*[1] 480*
*
*
*[[4]]*
*[1] 482*
*
*
*[
Hi
I have a question about drawing a linear line in x, y plot. I usually
use the following code, but for this time the x values are to small
(-0.08 to -0.02)
I wrote the following code, but r does not draw the line. However, it
does not give an error when it takes the code.
reg1<- lm(CWSI~NWI,
Can you trim down your example to a size where you can show us
the data (using dump() or dput()) and the commands you used so
one can just copy it from your mail and paste it into R to reproduce
the problem?
I don't see your problem when I made up data similar to what you described:
> timeseries
On 28-10-2013, at 16:07, Rui Barradas wrote:
> Hello,
>
> I can't reproduce your error:
>
> windfreq <-
> c(1351L, 2147L, 3317L, 4378L, 5527L, 6667L, 7865L, 8970L, 9987L,
> 10907L, 11905L, 12642L, 131000L, 14983L, 15847L, 16842L, 17757L,
> 18698L, 19632L, 20626L, 21599L, 22529L, 23325L, 24391L
Hello,
I can't reproduce your error:
windfreq <-
c(1351L, 2147L, 3317L, 4378L, 5527L, 6667L, 7865L, 8970L, 9987L,
10907L, 11905L, 12642L, 131000L, 14983L, 15847L, 16842L, 17757L,
18698L, 19632L, 20626L, 21599L, 22529L, 23325L, 24391L, 25356L,
26267L, 27230L, 28223L, 29190L, 30142L, 31124L, 32104
Hi,
A model gives me as output a netCDF file which I read with ncdf package. The
output consists of three similar matrices: one contains the variable value,
and the other two, the longitude and latitude coordinates. I need to do a
contour graph with such information, however, the R contour funct
Hi,
I'd like to follow up an older discussion on adding significance stars to
boxplots.
(
http://r.789695.n4.nabble.com/indicating-significant-differences-in-boxplots-td862335.html
)
As suggested by Jim Lemon, there is following solution for a single boxplot:
boxplot(count ~ spray, data = InsectS
Hello everyone,
This is Kangmin.
I am trying to produce shape and scale of my wind data. My data is based on
wind speed frequency with 1km/hr increment. data is described below.
Windspeed (km/h)Frequency
1 351
2 147
3 317
4 378
5 527
6 667
7 865
8 970
9 987
10 907
11 905
12 64
Dear Ben,
Thank you, I agree. I have forwarded this to the r-sig-mixed-models list.
Best,
Andreas
--
View this message in context:
http://r.789695.n4.nabble.com/Problems-with-lme-random-slope-intercept-model-tp4679104p4679168.html
Sent from the R help mailing list archive at Nabble.com.
Hello everyone,
This is Kangmin.
I am trying to produce shape and scale of my wind data. My data is based on
wind speed frequency with 1km/hr increment. data is described below.
Windspeed (km/h)Frequency
1 351
2 147
3 317
4 378
5 527
6 667
7 865
8
The output is as follows. My question is how to include only the interaction
terms where x is equal to "Maintained".
(x has two possible values "Maintained" and "Nonmaintained".)
--
View this message in context:
http://r.789695.n4.nabble.com/coxph-how-to-define-interaction-terms-tp4679162p46
Pavlos,
There are several ways to evaluate how well new data fit an old
regression.Part of the answer depends on what you are concerned about. For
example, if you are concerned about bias, you can test whether the mean of
the new data is within the expected range of the mean of that many new
valu
Dear all
I've got the following problem, I want to extract the residuals from
regression loops. The problem here is that some columns include NA's at the
beginning and end (i.e. each time series of stocks starts at different
points in time and ends at different points in time).
Dear list users,
I'm doing a plot integrating Grid output with Base Graphics output (gridBase,
Murrell 2012).
My goal is to produce a xy plot where each point is represented by a pie. I
could get it using the attach code,
but now I want to change the X axis to log scale.
when I introduce t
Hi Rolf,
Thanks for the response. I have re-phrased the problem. Hope this will help.
We're working on a project where our model tries to predict the value of
bookings as a time series of past bookings as well as some external flags.
The data was found to be stationary using the Dicky-Fuller tes
> Rolf Turner
> on Sat, 26 Oct 2013 09:34:14 +1300 writes:
> On 10/25/13 19:24, Alexandra Kuznetsova wrote:
>> Thank you for the reply! I need to use it in the package
>> for constructing the contrast matrices.
> Really?
>> Renaming a variable like y2 <- y^2
Rick,
I see a couple errors in my code ... you will need to attach the MASS
package for the eqscplot() function and that function's first arguments
should be any coordinates, I use (1, 1)
library(jpeg)
library(MASS)
# download image from internet for use in example
# http://www.tootsie.com/wallpa
Hi Hansol,
If you know the lattice system
library(lattice)
library(latticeExtra)
useOuterStrips(
xyplot(y ~ x|class, data = dat, ...)
)
See also ? combineLimits
Otherwise
mfrow(10,10)
and loop your data to plot
This is about as far as I can go - I have the same comment as Jim
Regards
Du
Hi experts,
Is there any packag for PAV and MPAV algorithms (Pattern Anomaly Value)
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://w
I’m trying to set up Cox Proptional Hazard model with interactions between
time and the covariates (which are categorical). The problem that I face is
that how to define the interactions, i.e. “x+cutStart:x”, properly.
The code below illustrates the problem. R gives the error message ” X matrix
de
Peter,
There is a function called codebook as part of the memisc package that probably
does what you want.
See also the package vignette:
http://cran.r-project.org/web/packages/memisc/vignettes/anes48.pdf
library(memisc)
fn <- "NILT12w2.por"
dat <- spss.portable.file(fn)
codebook(dat)
names(dat
Hi Jim and thanks for your answer... I might be too tired with my new born or
just exhausted.
I am sharing for everyone a small data snipset that you can load
https://www.dropbox.com/s/fh8jhwujgunmtrb/DataToPlotAsImage.Rdata
load("DataToPlotAsImage.Rdata")
require(plotrix)
browser()
test<-da
64 matches
Mail list logo
