I have installed R 2.7.2 in windows, now i am trying to install R oracle... I
am using Oracle 10g. i have downloaded the Roracle 0.5-9 sources, and i am
trying to compile it using Vc++, i found from readme files and forum that,
we need to use the makefile.win in the \src folder and nmake utility
Dear R users,
I am building a hierarchical model on a large data set. It can take quite
some time to finish one fit, I was just wondering whether it is possible to
store the fit object (the result) to a file for later (offline) analysis.
thanks
Julia
--
View this message in context:
Arthur Weiss wrote:
Hi everyone,
I am using the package spatstat for ploting kernel maps of my data.
It is a marked point pattern, the result of mosquito surveillance in a
area in a week.
For each trap, the number of individuals captured is the mark of the
point.
plot(density(X,
Dear R helpers,
Hi I am working on credit scoring model using logistic regression. I have main
sample of 42500 clentes and based on their status as regards to defaulted / non
- defaulted, I have genereted the probability of default.
I have a hold out sample of 5000 clients. I have calculated
liujb liujulia7 at yahoo.com writes:
I am building a hierarchical model on a large data set. It can take quite
some time to finish one fit, I was just wondering whether it is possible to
store the fit object (the result) to a file for later (offline) analysis.
Like with any other R object,
Hello,
If I have two correlation matrices (e.g. one for each of two treatments) and
then perform cor() on those two correlation matrices is this third
correlation matrix interpreted as the correlation between the two
treatments?
In my sample below I would interpret that the treatments are 0.28
Bing Shen wrote:
Hello,
I'm not sure this is doable but I'm having trouble running my R script
with multithreaded capability.
With 16x2.93Ghz CPUs available, only one is running with 100%. Any
suggestions?
Yes: Run 16 R jobs or use some other way to make your job parallel (e.g.
the snow
It worked for me, do you have the latest version of ggplot2 released a
few days ago (ggplot2_0.7) ?
Baptiste
On 9 Oct 2008, at 20:55, stephen sefick wrote:
Error in `[.data.frame`(df, , var) : undefined columns selected
I got this error in a fresh R session after rerunning all of the
I don't know where you got the idea that 'axes' was a valid argument.
From ?plot.Sample
...: further graphical parameters.
and 'axes' is *not* a graphical parameter. (It also suggests you call
plot, not plot.Sample.)
'xaxt' and 'yaxt' are graphics parameters which control whether
Hi everybody!
I am searching for an R implementation of multi-scale SSA (singular spectrum
analysis), which was introduced by Yiou etal (2000) (Data-adaptive wavelets
and *multi*-*scale *singular-spectrum
analysishttp://linkinghub.elsevier.com/retrieve/pii/S016727890452)
and further described
On 10-Oct-08 08:07:34, Michael Just wrote:
Hello,
If I have two correlation matrices (e.g. one for each of two
treatments) and then perform cor() on those two correlation
matrices is this third correlation matrix interpreted as the
correlation between the two treatments?
In my sample below
Richard.Cotton wrote:
The rectangles being drawn extend higher than the top of the panel.
(Your y
axis ranges from 0 to 50, but the bars go up to 100.)
Thankyou - I can also make the bars on the lower panels vanish by
tinkering with ylim.
In the top row of plots, depending upon the shape
I haven't heard Rpad mentioned yet, so I will mention it. Rpad allows you to
run R code inside a browser window. You can many any kind of GUI you want
using html forms and then call the R function from the html page. The
Apache/Rpad install works on any platform and anyone who uses it just
wmtsa has some wavelet processing things, but It sounds like you need
a very specific implementation. Sorry I would not be of more help-
good luck.
Stephen
On Fri, Oct 10, 2008 at 7:26 AM, Maik Renner [EMAIL PROTECTED] wrote:
Hi everybody!
I am searching for an R implementation of
Zitat von Duncan Murdoch [EMAIL PROTECTED]:
On 10/10/2008 8:13 AM, Oliver Bandel wrote:
Hello,
I tried to use rgl.snapshot and it failed.
The error message was not very verbose:
==
plot3d( motion[[idx+2]], motion[[idx+1]], motion[[idx]] )
Maithili Shiva wrote:
Dear R helpers,
Hi I am working on credit scoring model using logistic regression. I have main
sample of 42500 clentes and based on their status as regards to defaulted / non
- defaulted, I have genereted the probability of default.
I have a hold out sample of 5000
Hello,
I tried to use rgl.snapshot and it failed.
The error message was not very verbose:
==
plot3d( motion[[idx+2]], motion[[idx+1]], motion[[idx]] )
rgl.snapshot(filename=/tmp/shot_01.png, fmt=png)
[1] failed
==
There was a graphic created by
Hi,
I'm using the add-on package sound.I have the following
q-loadSample(a.wav)
plot.Sample(q,axes=FALSE)
this gives me the error.
Error in plot.default(sound(s)[1, ], type = l, col = red, ylim = c(-1,
:
formal argument axes matched by multiple actual arguments
I'm guessing the package has
[EMAIL PROTECTED] wrote:
r11 -- r16 are variables showing a reason for usage of a product in 6
different situations. Each variable is a factor with 4 levels imported
from a SPSS sav file with labels ranging from not important to very
important, and NA's for a sample of N = 276.
(1) I need a
Not to detract from John's very interesting work, but
there is a problem with citations as a measure of
penetration. When R totally succeeds, there will be
no citations. How many citations are there for Excel?
However, for the moment, citations are probably a
pretty good measure.
Using traffic
I am building a hierarchical model on a large data set. It can take
quite
some time to finish one fit, I was just wondering whether it is possible
to
store the fit object (the result) to a file for later (offline)
analysis.
See ?save.
Regards,
Richie.
Mathematical Sciences Unit
HSL
Hello, R users
I like to draw a cluster tree with a hclust and plot command.
But,instead of labels, it is giving number of variables in dendogram.
Could anyone tell me , how to print the labels of variables used for drawing.
Thank you.
--
Vijaykumar Yogesh Muley
Senior Research Fellow,
Dr.
--- On Fri, 10/10/08, Oliver Bandel [EMAIL PROTECTED] wrote:
From: Oliver Bandel [EMAIL PROTECTED]
Subject: Re: [R] plot-parameter pch without influence when plotting a
data-frame
To: [EMAIL PROTECTED]
what is the result of:
[...]
Try something like
mm - matrix(1:128,
Not an answer but I find that
DF - data.frame(let = letters[1:3], num = 1:3,
stringsAsFactors = FALSE)
is very handy. Damn it, if I want a factor I'll tell the machine I do :)
Or you can also set this option globally with
options(stringsAsFactors = TRUE)
Gabor
On Thu, Oct 9, 2008 at 8:19 AM, John Kane [EMAIL PROTECTED] wrote:
--- On Thu, 10/9/08, Oliver Bandel [EMAIL PROTECTED] wrote:
From: Oliver Bandel [EMAIL PROTECTED]
Subject: Re: [R] plot-parameter pch without influence when plotting a
data-frame
To: Gerhard Schön, UKE Hamburg [EMAIL
I would like to do the following in ggplot:
what am I missing?
River.Mile -c(202, 198, 190, 185, 179, 148, 119, 61)
TSS - c(1:8)
DOC - seq(2, by= 0.6, length.out=8)
z - data.frame(River.Mile, TSS, DOC)
xyplot(TSS+DOC~River.Mile, data=z, auto.key=TRUE)
thanks
--
Stephen Sefick
Research Scientist
Another solution would be
faclevels-c(Not,Little,Somewhat,Very, NA)
group - factor(1:6)
svdf-data.frame(value = sample(group, 400, TRUE), category =
sample(faclevels,400,TRUE))
library(ggplot2)
ggplot(svdf, aes(x = value, fill = category)) + geom_histogram(position
= dodge)
HTH,
Thierry
I want to make illustrations a la boxplot for a data material with
several groups. Some of the groups are very small, so a boxplot
doesn't make sense. I would like to use stripchart for these groups.
On the other hand, some of the groups are very large, so for them
stripchart isn't so good. So my
Given a Date object or simply a year is there an R function to tell me if the
it is a leap year or not? I was hoping for something like 'is.leapyear'. I
probably can build my own function (year divisible by 4 etc.) but I would
rather use an existing function if it is available.
Thank you.
Hi,
I would like to estimate a linear expendire with Systemfit package.
(method: SUR)
As someone could show me how to define the equations?
Thanks.
--
Think before you print !
**
Disclaimer: This e-mail may contain
Göran Broström wrote:
I want to make illustrations a la boxplot for a data material with
several groups. Some of the groups are very small, so a boxplot
doesn't make sense. I would like to use stripchart for these groups.
On the other hand, some of the groups are very large, so for them
Zitat von John Kane [EMAIL PROTECTED]:
--- On Thu, 10/9/08, Oliver Bandel [EMAIL PROTECTED] wrote:
From: Oliver Bandel [EMAIL PROTECTED]
Subject: Re: [R] plot-parameter pch without influence when plotting
a data-frame
To: Gerhard Schön, UKE Hamburg [EMAIL PROTECTED]
Cc: R-help
Dear Stephen,
You need the data in long format. Try melt()ing it first.
River.Mile -c(202, 198, 190, 185, 179, 148, 119, 61)
TSS - c(1:8)
DOC - seq(2, by= 0.6, length.out=8)
z - data.frame(River.Mile, TSS, DOC)
zMelt - melt.data.frame(z, id.vars = River.Mile)
library(ggplot2)
ggplot(zMelt, aes(x
Hi everybody!
I am searching for an R implementation of multi-scale singular
spectrum analysis (SSA), which was introduced by Yiou etal (2000)
(Data-adaptive wavelets and multi-scale singular-spectrum analysis)
and further described in Ghil etal (2002).
For SSA alone I found an package recently
Hello,
I am unable to start an R session by double-clicking on an existing
.RData file from the Windoze file explorer.
A dialogue appears with the message
Fatal error: unable to restore saved data in .RData,
and I notice in the R Console the last few lines are:
Loading required
Hi Marie!
On Friday 10 October 2008 12:40:23, Marie Vandresse wrote:
I would like to estimate a linear expendire with Systemfit package.
(method: SUR)
If I remember correctly, the linear expenditure system (LES) is linear in
income but non-linear in the parameters. Hence, you have to estimate
Hello,
I may have found the problem, although I don't understand it.
When I commented out the following line in my Rprofile.site file the
problem _seemed_ to disappear.
# ans as shorthand for .Last.value
# makeActiveBinding(ans, function() .Last.value, .GlobalEnv)
I decided to do
On Fri, Oct 10, 2008 at 9:18 AM, stephen sefick [EMAIL PROTECTED] wrote:
I would like to do the following in ggplot:
what am I missing?
River.Mile -c(202, 198, 190, 185, 179, 148, 119, 61)
TSS - c(1:8)
DOC - seq(2, by= 0.6, length.out=8)
z - data.frame(River.Mile, TSS, DOC)
it did help to create the separate vectors. The problem is, I would like
to continue working with each vector on its own and make some further
analysis.
For that i need to allocate each vector in the filePatterns a unique name,
which I just can't do.
I'm not so great in R loops and the help
Hello. I have an instrument, actually a step motor indexer that I want to drive
from R. It uses an RS-232 connection. I already have a few C subroutines that I
access from R but I would like to be able use R's connection related functions
- readLines, writeLines, cat, read.table, etc -
If you are able to create a .Rdata file for which this fails can you
record the last few commands before you create the .Rdata file?
luke
On Fri, 10 Oct 2008, Thompson, David (MNR) wrote:
Hello,
I am unable to start an R session by double-clicking on an existing
.RData file from the Windoze
I just got it and it works wonderfully.
thank you both for your help
Stephen
On Fri, Oct 10, 2008 at 4:01 AM, baptiste auguie [EMAIL PROTECTED] wrote:
It worked for me, do you have the latest version of ggplot2 released a few
days ago (ggplot2_0.7) ?
Baptiste
On 9 Oct 2008, at 20:55,
Thanks, Duncan, I was encountered with the same problem as Michael,
and just now I found the version you provided in your homepage could
solve this problem well.
Regards,
Yihui
--
Yihui Xie [EMAIL PROTECTED]
Phone: +86-(0)10-82509086 Fax: +86-(0)10-82509086
Mobile: +86-15810805877
Homepage:
On 10/10/2008 8:13 AM, Oliver Bandel wrote:
Hello,
I tried to use rgl.snapshot and it failed.
The error message was not very verbose:
==
plot3d( motion[[idx+2]], motion[[idx+1]], motion[[idx]] )
rgl.snapshot(filename=/tmp/shot_01.png, fmt=png)
[1] failed
rkevinburton at charter.net writes:
Given a Date object or simply a year is there an R function to tell me if the
it is a leap year or not? I was hoping
for something like 'is.leapyear'. I probably can build my own function (year
divisible by 4 etc.) but I
would rather use an existing
Dear Paulo,
Hello. I have an instrument, actually a step motor indexer that I want to
drive from R. It uses an RS-232 connection. I already have a few C
subroutines that I access from R but I would like to be able use R's
connection related functions - readLines, writeLines, cat, read.table,
Hi,
I thought about this but programming it seems rather difficult so I was
wondering if a function exists for this in R (as most of the times it turns out
that it does):
I have a map (shapefile) and for about 50 points on that map (GPS locations) I
have heights. Is there a function which can
On Fri, 10 Oct 2008, Thompson, David (MNR) wrote:
Hello,
I may have found the problem, although I don't understand it.
When I commented out the following line in my Rprofile.site file the
problem _seemed_ to disappear.
# ans as shorthand for .Last.value
# makeActiveBinding(ans,
Dear Ben Bolker,
Thanks a lot for your help.
I have two more questions:
1) My goal is maximize the function ( r = Y*log(comb)) but the
parameters found, minimized the function (r = Y*log(comb)).
2) What this function do?
model2 - function(p) {
do.call(Model,as.list(p))
}
Best
Dear R-Users,
I'm trying to estimate GARCH parameters implied by options prices (on
wednesdays only) by minimizing the average mean squared dollar error between
the market and the model price, but I always get the following error message:
Error in chol.default(Hessian) :
the leading minor
Hi,
I have three vectors and I would like to plot them in one window. So R
should overlay them in the
quartz. I achieved this by plotting the first vector: plot(v1) and add the
other two with
points(v2) and points(v3). Is there another way to achieve the same result?
Cheers
--
View this
Hi,
thanks again for your help.
I played around with your suggestions to change the appearance of my graph,
but there is still a problem I could not fix.
The vector, which I want to plot, contains 60 entries. After 3 bars I want
to have a large gap
between the next 3 bars. But what always
Dear all,
I am experiencing problems with glmmmPQL. I am trying to analyze
binomial data with some spatial autocorrelation. Here is my code and
some of the outputs
colnames(d.glmm)
[1] BV Longitude Latitudenb_pc_02 nb_expr_02
[6] pc_02 nb_pc_07nb_expr_07 pc_07
I played around with your suggestions to change the appearance of my
graph,
but there is still a problem I could not fix.
The vector, which I want to plot, contains 60 entries. After 3 bars I
want
to have a large gap
between the next 3 bars. But what always happens is, that R took at the
perhaps leapyear()
http://finzi.psych.upenn.edu/R/library/pheno/html/leapyear.html
--- On Fri, 10/10/08, [EMAIL PROTECTED] [EMAIL PROTECTED] wrote:
From: [EMAIL PROTECTED] [EMAIL PROTECTED]
Subject: [R] Leap year?
To: r-help r-help@r-project.org
Received: Friday, October 10, 2008, 10:29
List,
I have a question about why plots sent to certain devices can be of so much
worse quality than others. Whenever I plot to x11 or pdf the graphs are of
excellent quality. eps is a bit poorer, but of passable quality, and jpeg
or tiff are terrible. I tried all manner of parameters, on many
Vijaykumar Muley vijay.muley at gmail.com writes:
I like to draw a cluster tree with a hclust and plot command.
But,instead of labels, it is giving number of variables in dendogram.
Could anyone tell me , how to print the labels of variables used for drawing.
Thank you.
Please post str(of
?matplot
On Fri, Oct 10, 2008 at 9:50 AM, mentor_ [EMAIL PROTECTED] wrote:
Hi,
I have three vectors and I would like to plot them in one window. So R
should overlay them in the
quartz. I achieved this by plotting the first vector: plot(v1) and add the
other two with
points(v2) and
Hi, I would like to use the quantile function, but I have some doubts to
choose the value of type parameter. Does this depend on the
shape/structure of my distribution? I cannot unde rstand what mean the
expressions Discontinuous sample quantile and Continuous sample quantile,
respectively.
My
Dear R-users
When running a glm polynomial model with one explanatory variable (example
Y~X+X^2), with a poisson or binomial error distribution, the predicted values
obtained from using the predict() function and those obtained from using the
coefficients from the summary table as is in an
I am trying to lag a time series. My data is in a matrix, but I coerce it
into a ts object.
But when I lag it and then look at the result, nothing has changed. What am
I doing wrong?
residsq-resid^2
residsq-as.ts(residsq)
lag1residsq-lag(residsq,-1)
residsq[1:5]
1 2
Hello R-Users,
I have recently run into several problems using vgam() in the VGAM
package. I am hoping someone might have some solutions...
Briefly, I have been trying to fit GAM models for zero-altered negative
binomial models.
1. When fitting smoothed parameters (e.g. s(X, df=2))
lag1resid is a time series, try
str(residsq)
str(lag1residsq)
The problem is that when you subscript lag1resid you
don't get a time series out from that. See ?window.ts
or try the zoo or xts packages where subscripting of time
series works.
On Fri, Oct 10, 2008 at 1:27 PM, Scotty Nelson [EMAIL
I don't know what you mean by XCoef x X. But your problem is (as it works if
you specify normal in a glm) that the functional relationship between your
predictors, i.e. Intercept+X+X^2, and Y is not linear for a binomial or a
poisson distribution.
Generalized linear model implies that the model
Werner Wernersen pensterfuzzer at yahoo.de writes:
Hi,
I thought about this but programming it seems rather difficult so I was
wondering if a function exists for
this in R (as most of the times it turns out that it does):
I have a map (shapefile) and for about 50 points on that map (GPS
Thanks, really it is difficult to find solution
Jarek
Peter Dalgaard pisze:
Jarek Jasiewicz wrote:
Hi
It is probably simple but how to force in plot command number format:
2, 3, 4 etc. instead of 2e+04, 3e+04 etc.
Not all that simple! It's one of those cases where you need to
thanks. i'll try zoo or xts. ts is a p.i.t.a.
Gabor Grothendieck wrote:
lag1resid is a time series, try
str(residsq)
str(lag1residsq)
The problem is that when you subscript lag1resid you
don't get a time series out from that. See ?window.ts
or try the zoo or xts packages where
Eduardo Marinho eduardomarinho at gmail.com writes:
Dear all,
I would like to use a GWR model in order to spatially predict food
insecurity in Africa. I have a georeferenced village data-bases and I've run
a classic regression model (taking into account the spatial dependence of
the
Note that there is actually a good reason why this is so.
ts can only represent regularly spaced series but the
result of a subscripting operation might not be a regularly
spaced series so you would not be able to guarantee that
it could be a ts series as well. That's the reason for window.ts
Jean-Baptiste Ferdy jean-baptiste.ferdy at univ-montp2.fr writes:
Dear all,
I am experiencing problems with glmmmPQL. I am trying to analyze
binomial data with some spatial autocorrelation. Here is my code and
some of the outputs
m.1 -
LFRC feliperiehs at yahoo.com.br writes:
Dear Ben Bolker,
Thanks a lot for your help.
I have two more questions:
1) My goal is maximize the function ( r = Y*log(comb)) but the
parameters found, minimized the function (r = Y*log(comb)).
Oh. Oops. Just change the sign ( r
Yihui Xie wrote:
Thanks, Duncan, I was encountered with the same problem as Michael,
and just now I found the version you provided in your homepage could
solve this problem well.
Thanks very much, Yihui --
Using Duncan's homepage version also worked for me, and your test helped
me slay the
I have done what P. Dalgaard has suggested and I don't find a
descrepancy between the number of values and the number of labels: there
15 each...
Any hint on what might go wrong here ?
Here is the output
The SAS format from proc contents
VISITF
Hi Guys,
I am in desperate need with sampling. I am suppose to sample from a dataset.
I use the following code:
for(i in 1:5){temp[i]-sample(T2,40,replace=F)
+ show(temp)}
but all the samples are the same, but I want them to be different.any
suggestion?
Thanks
--
View this message in
Jean-Louis Abitbol wrote:
I have done what P. Dalgaard has suggested and I don't find a
descrepancy between the number of values and the number of labels: there
15 each...
Any hint on what might go wrong here ?
Actually, I think you got it:
factor(1,c(NA,1:4),c(1:5))
Error in factor(1,
You appear to have a label for an NA value, and that's the problem.
I think you need to send an example to the package maintainers (see the
posting guide) to help them debug what looks like a bug or undocumented
restriction in their code (and it seems this is the same codebase in the
two
Hi R community,
I'm running R 2.7.2 on Windows XP SP2.
I'm trying to (1) plot loess lines for each of my groupings using the same
color for each group; (2) plot loess predicted values.
The first part is easy:
data1 -
Hi all,
I am using the function plotCI with the following command:
plotCI(m.residuos.p.2 [1:41],li=m.residuos.p.3 [1:41],ui=m.residuos.p.4
[1:41],lty=1,ylab=)
This generates exactly what I want except for the fact that I wanna drawn a
line linking the points (m.residuos).
How could I do that?
On 11/10/2008, at 10:15 AM, Caio Azevedo wrote:
Hi all,
I am using the function plotCI with the following command:
plotCI(m.residuos.p.2 [1:41],li=m.residuos.p.3
[1:41],ui=m.residuos.p.4
[1:41],lty=1,ylab=)
This generates exactly what I want except for the fact that I wanna
drawn a
Thanks, Roger,
It seems simply great to me. I am far from being able to detect the
drawbacks of the system, if any. I don't understand why this kind of
programs are not of general use...
With my best regards,
Ricardo
Bos, Roger wrote:
I haven't heard Rpad mentioned yet, so I will mention
Hi,
I was wondering if there is a function in R which performs bivariate non
parametric smoothing
which allows for the possibility of including some weights in the smoothing
(for each data points
in my grid I have some predefined weights that I would like to include in the
smoothing).
On 10/10/2008 8:55 AM, Oliver Bandel wrote:
Zitat von Duncan Murdoch [EMAIL PROTECTED]:
On 10/10/2008 8:13 AM, Oliver Bandel wrote:
Hello,
I tried to use rgl.snapshot and it failed.
The error message was not very verbose:
==
plot3d( motion[[idx+2]], motion[[idx+1]],
Would it be possible to get data from uk.finance.yahoo.com instead
from finance.yahoo.com.
Thanks
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide
I'm trying to use R to determine the quality of a cubic Bezier curve
approximation of an elliptical arc. I know the four control points
and I want to compute (x,y) coordinates of many points on the curve.
I can't find anything in either the base distribution or CRAN that
does this; all the
Hi
I have to teach a course of R in two day , I need some an advice about what
must I teach?
Thanks in advance
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do
85 matches
Mail list logo
