This is really off-topic here (a programming question about a
non-recommended compiler system, and about a contributed package), but
reading README.packages will give you some clues (as rw-FAQ 8.3 says).
The current version is at
https://svn.r-project.org/R/trunk/src/gnuwin32/README.packages
P
Use ts() in R. See e.g. the chapter on times series in MASS4.
Just replacing rts() by ts() in your example gives what I guess you were
expecting in S-PLUS.
On Wed, 30 Apr 2008, Louise Hoffman wrote:
Dear readers,
S-Plus have a rts function (Regular Time Series), which is used like so:
fvek
I am reading some documentation about Cross Spectrum Analysis as a technique
to compare spectra.
My understanding is that it estimates the correlation strength between
quasi-periodic structures embedded in two signals. I believe it may be
useful for my signals analysis.
I was referred to the R
On Tue, Apr 29, 2008 at 11:03 PM, Louise Hoffman
<[EMAIL PROTECTED]> wrote:
> Dear readers,
>
> S-Plus have a rts function (Regular Time Series), which is used like so:
>
> fveks<-read.csv('http://louise.hoffman.googlepages.com/veks.csv',header=TRUE,sep=',')
> attach(fveks)
> acf(ts.intersect(rts(H
Georg,
A boxplot does not have a mean-bar. It has a median bar.
Bill Venables
CSIRO Laboratories
PO Box 120, Cleveland, 4163
AUSTRALIA
Office Phone (email preferred): +61 7 3826 7251
Fax (if absolutely necessary): +61 7 3826 7304
Mobile: +61 4 8819 4402
Home Phone:
?boxplot
The boxplot does not show the mean but the median:
median(c(0,15,0,60,0,0,0,0,0,60,60,0,60,0,30,0))
> b <- boxplot(c(0,15,0,60,0,0,0,0,0,60,60,0,60,0,30,0))
b$stats
[,1]
[1,]0
[2,]0
[3,]0
[4,] 45
[5,] 60
_
Professor Michael Kubovy
Univ
Dear R user group,
I am working with boxplots and cannot solve the following problem:
> data<-c(0,15,0,60,0,0,0,0,0,60,60,0,60,0,30,0)
> data
[1] 0 15 0 60 0 0 0 0 0 60 60 0 60 0 30 0
> boxplot(data)
The boxplot has the first hinge and the mean-bar both at 0... Why is this
the mean
Dear readers,
S-Plus have a rts function (Regular Time Series), which is used like so:
fveks<-read.csv('http://louise.hoffman.googlepages.com/veks.csv',header=TRUE,sep=',')
attach(fveks)
acf(ts.intersect(rts(HC.f),rts(Ta.f),rts(GR.f),rts(W.f)))
Warning the csv file is 750kB.
Can the same be don
Is this what you want:
> ES <- 1
> ifelse(length(ES %in% c(1,2,3)) == 0, TRUE, ES %in% c(1,2,3))
[1] TRUE
> ES <- 4
> ifelse(length(ES %in% c(1,2,3)) == 0, TRUE, ES %in% c(1,2,3))
[1] FALSE
> ES <- NULL
> ifelse(length(ES %in% c(1,2,3)) == 0, TRUE, ES %in% c(1,2,3))
[1] TRUE
>
On Tue, Apr 29, 20
brevity is refreshing
On Tue, Apr 29, 2008 at 3:40 PM, Mittal Meghna <[EMAIL PROTECTED]> wrote:
>
> __
> 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/posti
-BEGIN PGP SIGNED MESSAGE-
Hash: SHA256
Dear List:
I'm looking for a form of the empty set such that if ES is said
representation,
~ES %in% c(1,2,3)
evaluates to TRUE.
Thank you in advance for your assistance.
Sincerely,
Jason Q. McClintic
- --
Jason Q McClintic
UST MB 1945
2115
I'd just like to thank all you guys for stepping in so promptly with help. I
haven't yet had a chance to implement any of your code yet, but just by
looking over what you've suggested, I think I have enough to guide me. So
thanks once again!
--
View this message in context:
http://www.nabble.com
On 29/04/2008 6:52 PM, Sebastian Eiser wrote:
Hello R-world!
This is my first post to the R list, so I feel that I need to say thanks to
the community (I'm using R now for 5 months)!
I'd greatly appreciate if somebody could explain some odd behavior of
".Fortran" to me. Let's assume a primitive
Hi,
when I try installing new packages (in this case DBI) I run into the
following problem:
Paket 'DBI' erfolgreich ausgepackt und MD5 Summen abgeglichen
Warnung: kann temporäre Installation
'C:\Programme\R\R-2.7.0\library\file5f906952\DBI' nicht nach
'Ý xlDÌú [EMAIL PROTECTED]' verschieben
Hello R-world!
This is my first post to the R list, so I feel that I need to say thanks to
the community (I'm using R now for 5 months)!
I'd greatly appreciate if somebody could explain some odd behavior of
".Fortran" to me. Let's assume a primitive Fortran subroutine
barfoo.f
--
> results[i]<- nls(Tw ~ mu + ((alpha - mu)/(1 + exp(gamma*(B -
> Mean_Air,
> data = dem16,
> start = list(mu = 0.0001, alpha = 21.8, gamma = 0.22, B = 12.8))
> }
If you have a variable that codes "site" then you can try something
like this to get the parameters over sites.
Has anyone come up with some routines that will perform kernel matching
pursuit, or orthogonal matching pursuit? I was hoping for something similar
to what matlab has with wavelab:
http://www-stat.stanford.edu/~wavelab/
Best,
-Justin
--
Justin Donaldson
PhD Candidate, Informatics
Indiana Unive
Thanks Kingsford! I will cc r-help.
varPower() works for me. I want to use lm because predict.gls does not give
lwr and upr values, and also the bptest and ncv.test only work with lm
models...
On 4/29/08, Kingsford Jones <[EMAIL PROTECTED]> wrote:
>
> hi Tom,
>
> Basically, a simple way to model
Much thanks to Chuck Cleland for the following solution:
t(sapply(split(mydf, mydf$TIME), function(x){coefficients(summary(glm(X ~ A
+ B, data = x, family=binomial)))}))
To see what's in the 12 columns of the matrix returned by those lines, look
at the results this way:
lapply(split(mydf, mydf
Hi,
when I try installing new packages (in this case DBI) I run into the
following problem:
Paket 'DBI' erfolgreich ausgepackt und MD5 Summen abgeglichen
Warnung: kann temporäre Installation
'C:\Programme\R\R-2.7.0\library\file5f906952\DBI' nicht nach 'Ý xlDÌú
[EMAIL PROTECTED]' verschieben
If you can(one dimensional only), try using lowess() instead. Probably in a
for loop as Ray suggested.
loess() is more powerful and flexible, but you pay for it in extra
complexity and time. Maybe in this case, it's not worth it.
-- Bert Gunter
Genentech
-Original Message-
From: [EMAIL P
I'm familiar with mazes that have a designated starting and ending points
but don't see any here (or is this just a maze
Clint BowmanINTERNET: [EMAIL PROTECTED]
Air Dispersion Modeler INTERNET: [EMAIL PROTECTED]
Air Quality Program VOICE:
Hi,
I am new to R. I need to communicate with Oracle from R program in windows
xp. So I am planning to use ROracle. I downloaded the ROracle src from the
site below:
http://cran.r-project.org/web/packages/ROracle/index.html
I am using R version 2.6.2 and Oracle 10g. I followed the instruct
In addition to Tony's suggestion, have a look at the following sequence, which
I suspect is because the call to apply will duplicate your 1.5GB matrix,
whereas the for loop doesn't [I stand to be corrected here].
> x <- matrix(runif(21), 21)
> unix.time({res <- numeric(ncol(x)); for(i in 1:l
Thanks! I wasn't on a clean R session so I didn't catch this. I intended the
'i' to be the same from the last for loop, but as you point out on the first
run there is no 'i'!
Try this:
plotMaze <- function(z, text=FALSE, lwd=2, ...) {
N <- nrow(z)
y <- -matrix(rep(1:N, N), N, N, byrow=F
Hi,
>
> Any one has the code that I can borrow? Or any suggestions?
>
> Highly appreciate.
>
> Hongyan Du
__
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/post
It's quite possible that much of the time spent in loess() is setting up the
data (i.e., the formula, terms, model.frame, etc.), and that much of that is
repeated identically for each call to loess(). I would suggest looking at the
code of loess() and work out what arguments it is calling simp
McGehee, Robert wrote:
I had some fun this afternoon coding up a 'maze generator' in R. I
thought I'd pass along the fruits of my labor for everyone's amusement.
As written, every point is connected to every other point, so feel free
to 'start' and 'finish' anywhere you like.
Have fun!
--Rober
Sorry about this. cairoDevice 2.8 (just uploaded to CRAN) should fix this
problem.
On Tue, Apr 29, 2008 at 1:08 PM, Josh Gilbert <[EMAIL PROTECTED]> wrote:
> I sent this to R-Help and the listed maintainer of cairoDevice, I hope
> that was the right thing to do.
>
> For some reason, Cairo_png put
Hello,
I am trying to apply a least squares non-linear regression to my
dataset, dem16. I can apply it to a subset based on a single site
fine, but I want to apply it to each (of 197) sites. I am stumbling
over a "for" loop. Site and Mean_Air are variables in my dataset and
I would like to
Hi, I need to fit models and use AIC method to campare the best fitted model
manually.
When i extract AIC by using extractAIC, it gave me the df and AIC values.
Now the problem is, how can I compare the AIC values from two models?
is there anyway to extract AIC with no df so that I can compare dir
Jim
Thanks again
yes its 64 bit version of R.
So you suggest writing out 9 million columns and then loading
them individually and applying the function of each of these
columns!
But since we are using a for loop here too so won't it end up
taking almost same? Plus we aren't we increasing the i/o
ov
Hello all R-helpers,
I've performed an experiment to test for differential effects of
elevated temperatures on three different groups of corals. I'm
currently performing a cox proportional hazards regression with
censoring on the survivorship (days to mortality) of each individual
in the
I had some fun this afternoon coding up a 'maze generator' in R. I
thought I'd pass along the fruits of my labor for everyone's amusement.
As written, every point is connected to every other point, so feel free
to 'start' and 'finish' anywhere you like.
Have fun!
--Robert
PS. Feel free to pass
Beck, Kenneth (STP) <[EMAIL PROTECTED]> [Tue, Apr 29, 2008 at 07:21:40PM CEST]:
> mxx=max(cpx_list$nMV);
> mxy=max(trend_list$nMV);
> if (mxx>mxy)
> mxy=mxx
> else
> mxx=mxy
>
Can't this be replaced by
mxx <- max(c(cpx_list$nMV, trend_list$nMV))
mxy <- mxx
?
--
Johannes Hüsing
are the "chunks" on which you need to apply the function rolling windows?
do they overlap?
I have some c++ template utilities that I use for window functions (on
timeseries objects) which you are welcome to copy and modify to fit your
problem.
they are available here:
git://repo.or.cz/fts.git
gi
Try:
mxx = 2
mxy = 4
if (mxx > mxy) {
print ("mxx reigns supreme")
} else {
print ("mxy reigns supreme")
}
-Original Message-
From: [EMAIL PROTECTED] [mailto:[EMAIL PROTECTED]
On Behalf Of Beck, Kenneth (STP)
Sent: Tuesday, April 29, 2008 9:22 AM
To: r-help@r-project.org
Subject: [R
I apologize for sending this twice, I didn't realize that gmail
helpfully defaulted to rich formatting.
For some reason, Cairo_png puts a box around a figure when you call
plot.new. It looks like box was called with black and a transparent
background. Example:
> library(cairoDevice)
> Cairo_png('c
I sent this to R-Help and the listed maintainer of cairoDevice, I hope that
was the right thing to do.
For some reason, Cairo_png puts a box around a figure when you call
plot.new. It looks like box was called with black and a transparent
background. Example:
> library(cairoDevice)
> Cairo_png('ca
Hi Jim,
Thanks for your prompt response,
I am using a fairly powerful Mac with Leopard OS and 17GB RAM
and 2x3 GhZ intel zeon processor so I do not think the system
is paging. I also using the Rmpi and snow utilities to
parallelize it but even then it takes 3.5-4 hours to just
complete one chunk o
I'm afraid I haven't quite got the hang of this yet. Here's a Hello
World:
#include
#include
#include
#include
extern int R_running_as_main_program; /* in ../unix/system.c */
int main(int ac, char **av)
{
SEXP aleph, beth;
// R_running_as_main_program = 1; // A
//
Hi All,
I am calling some c code from R.
It successfully makes the dll and .so files.
When I run .so in Linux is works prefect but hangs in windows. Though the
dll loads but it never returns back from the c function
Can you please suggest the possible cause of this
thank you
vidhu
[[alte
__
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.
Thanks for the help Patrick
I ran
yum provides libg2c*
Looking in available packages for a providing package
Available package: gcc-g77.i386 0:3.4.6-9 from base matches with
/usr/lib/gcc/i386-redhat-linux/3.4.3/libg2c.a
Available package: gcc-g77.i386 0:3.4.6-9 from base matches with
/usr/lib/g
Thank you Jorge, what I want to have is
> z <- A1+A2
>z
4
The problem is that A1, A2,... still unknown, so that I have the following error
> A1
object "A1" not found
Thank you in advance!
- Original Message
From: Jorge Ivan Velez <[EMAIL PROTECTED]>
To: Diego Culattoni <[EMAIL PROTEC
Lisa wrote:
Hi, I have a question about how to extract paramters from a fitted model. I
can extract coefficients and std, but from some other statistics, I dont
know how to extract. Can anyone help?
Here it is an example:
coxout<-coxph(Surv(t,t.censor)~x)
coxout
Call:
coxph(formula = Su
I am not entirely sure what it is that you want to do.
functionname<-function(a,b,c){
+ a<-sum(a)
+ b<-sum(b)
+ c<-sum(c)
+ y<-(a+bx+cx^2)
+ x<-(x)
+ }
for each <- your are defining and object. it doesn't make much sense to
sum(a) unless you have made a something etc. Look at packages- type CRA
Respected R experts,
I am trying to apply a user function that basically calls and
applies the R loess function from stat package over each time
series. I have a large matrix of size 21 X 900 and I need
to apply the loess for each column and hence I have
implemented this separate user function
On an x86 machine operating under Gentoo Linux, I'm having trouble with
tcltk in R-2.6.1 (the latest version available through Gentoo's
emerge/portage system). The problem first came to my attention when
"install.packages("Rmetrics")" produced the following error message:
Loading required packa
Hi Diego,
Try this:
z<-As["A1"]+As["A2"]; names(z)<-"z"
z
4
HTH,
Jorge
On Tue, Apr 29, 2008 at 2:57 PM, Diego Culattoni <[EMAIL PROTECTED]>
wrote:
> Thank you Jorge, what I want to have is
> > z <- A1+A2
> >z
> 4
> The problem is that A1, A2,... still unknown, so that I have the following
See the source code of function:
getS3method("print", "coxph")
is this block:
tmp <- cbind(coef, exp(coef), se, coef/se, signif(1 -
pchisq((coef/se)^2, 1), digits - 1))
On 4/29/08, Lisa <[EMAIL PROTECTED]> wrote:
>
> Hi, I have a question about how to extract paramters from a fitted model.
> I
Hi, I have a question about how to extract paramters from a fitted model. I
can extract coefficients and std, but from some other statistics, I dont
know how to extract. Can anyone help?
Here it is an example:
> coxout<-coxph(Surv(t,t.censor)~x)
> coxout
Call:
coxph(formula = Surv(t, t.censor
Kenneth -
See ?if in the Details section. Specifically this part (at least in R
2.7),
"In particular, you should not have a newline between '}' and
'else' to avoid a syntax error in entering a 'if ... else'
construct at the keyboard or via 'source'."
The R interpreter can't 'see ahead' th
On 4/29/2008 1:21 PM, Beck, Kenneth (STP) wrote:
Why will this simple statement not work? I think I am following the
documentation for if(cond) statements, and I have tried wrapping the
cons.expr and alt.expr in {}, I get the same error. There is no example
in the help file, and this is not cover
FYI
This is a nice package that is not on CRAN for some reason. You can get it
at
http://www.csm.ornl.gov/esh/statoed/
Tom
--
View this message in context:
http://www.nabble.com/Statistical-Methods-and-Software-for-the-Analysis-of-Occupational-Exposure-Data-With-Non-Detectable-Values-tp169662
Hi Diego,
Is this what you want?
# Data set
set.seed(123)
x1<-paste("A", 1:6, sep = "")
x2<- round(rgamma(6,2,1))
x3<-paste("B", 1:6, sep = "")
x4<- round(rgamma(6,2,1))
data1 <- data.frame(x1,x2,x3,x4)
data1
x1 x2 x3 x4
1 A1 1 B1 0
2 A2 3 B2 0
3 A3 0 B3 3
4 A4 2 B4 2
5 A5 4 B5 2
6 A6
On 29/04/2008 12:19 PM, Maximillian Murphy wrote:
Dear All,
I've read the manual on "Writing R Extensions" and in particular the
part on calling R from C. (Most of the manual is about calling C
from R, not the other way around.)
The good news is that I can now call _some_ R from C, speci
I am try to create a plot using xyplot(). I created a function panel.fun() wich
generate segment, but I need use two columns from dataset2.
library(lattice)
panel.fun <- function(x, y, minx, maxx, miny, maxy)
{
panel.xyplot(x,y)
# since I don't know how to call dataset2 inside panel
> Hi,
>
> Any one has the code that I can borrow? Or any suggestions?
>
> Highly appreciate.
>
> Hongyan Du
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read
Why will this simple statement not work? I think I am following the
documentation for if(cond) statements, and I have tried wrapping the
cons.expr and alt.expr in {}, I get the same error. There is no example
in the help file, and this is not covered in the Introduction to R,
SimpleR or other tutor
On Tue, Apr 29, 2008 at 6:20 AM, tom soyer <[EMAIL PROTECTED]> wrote:
> Hi,
>
> I would like to use a weighted lm model to reduce heteroscendasticity. I am
> wondering if the only way to generate the weights in R is through the
> laborious process of trial and error by hand. Does anyone know if
Hey Anthony,
My previous function may not work in all cases. Say one of the
experiments yields these numbers:
1,2,3,6,7
Would you say that the proportion of consecutive numbers is 100%? If
so, this will work:
prop.diff=function(x){
d=diff(sort(x))
prop=sum((c(0,d==1)+c(d==1,0))>0)
prop=pr
On Tue, 29 Apr 2008, Anthony28 wrote:
I need to use R to model a large number of experiments (say, 1000). Each
experiment involves the random selection of 5 numbers (without replacement)
from a pool of numbers ranging between 1 and 30.
What I need to know is what *proportion* of those experime
Sorry I forgot these informations. I'm using R 2.70, under Linux (Debian
Etch). The png driver uses cairo, and the version of libcairo installed
is 1.2.4 (which is not the latest, but more recent versions are not
available is the debian stable repositories)
I tried bitmap(), and it does not
Hi Richard,
Thanks a lot! What I actually want to have
>A1
4
>A2
2
. and so on! Is this possible?
Thank you again :)!
- Original Message
From: "[EMAIL PROTECTED]" <[EMAIL PROTECTED]>
To: Diego Culattoni <[EMAIL PROTECTED]>
Cc: [EMAIL PROTECTED]; [EMAIL PROTECTED]
Sent: Monday
Dear All,
I've read the manual on "Writing R Extensions" and in particular the
part on calling R from C. (Most of the manual is about calling C
from R, not the other way around.)
The good news is that I can now call _some_ R from C, specifically
the R functions which have C header file
Hey Anthony,
There must be many ways to do this. This is one of them:
#First, define a function to calculate the proportion of consecutive
numbers in a vector.
prop.diff=function(x){
d=diff(sort(x))
prop=(sum(d==1)+1)/length(x)
return(prop)}
#Note that I am counting both numbers in a consecu
See help(warnings) and from there you'll also find help(warning). It
is described there.
/Henrik
On Tue, Apr 29, 2008 at 5:14 PM, Vidhu Choudhary
<[EMAIL PROTECTED]> wrote:
> Hi All,
> I want to hide all the system(R) generated warning messages and want only
> the messages from my script shoul
29/04/2008 17:04 trhermes a écrit
> I am outputting some graphs from SpatStat using a for loop. I want the png
> files to include the i and j from the nested loops in their names. For
> example, I want the file names to look like "i-j_plot.png". The code is
> below:
>
> ma <- levels(marks(X))
>
This will work:
my.list <- c(2, 28, 31, 4, 27)
sort(my.list)
diff(sort(my.list))
any(diff(sort(my.list)) == 1)
the middle two lines are only to illustrate what's going on.
Best wishes!
Charles Annis, P.E.
[EMAIL PROTECTED]
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEnginee
Richard Pearson wrote:
?rm
Richard
Ralf Goertz wrote:
How can I automatically exclude one variable from being saved in the
workspace when I quit an R session? The problem is I don't know how to
erase a variable once it has been created.
[...]
Ralf
More on the use of rm. If you want to
I am outputting some graphs from SpatStat using a for loop. I want the png
files to include the i and j from the nested loops in their names. For
example, I want the file names to look like "i-j_plot.png". The code is
below:
ma <- levels(marks(X))
n <- length(ma)
for(i in 1:n)
for(j in 1:
Two things:
1. the + Sex term is superfluous
2. the variable Subject needs to be a Factor, not a vector (as I
suspect it currently is). That is, add:
mydata.tab$Subject=as.factor(mydata.tab$Subject)
in the preamble before the aov call, and all should be fine.
On 29-Apr-08, at 4:00 AM, [EMA
On 29/04/2008 11:39 AM, Kenn Konstabel wrote:
Suppose X is a long vector of integers (typically about 3 elements). Is
there an efficient way to detect whether there are at least N consecutive
zeros in X, and if yes, where does this occur?
for example, suppose X is:
...
and N is 10. I wou
Suppose X is a long vector of integers (typically about 3 elements). Is
there an efficient way to detect whether there are at least N consecutive
zeros in X, and if yes, where does this occur?
for example, suppose X is:
1 2 3 4 5 0 0 0 0 0 0 0 0 0 0 0 0 0 0 3 4 2 0 1 2 0 0 0 2 2 2 2 2 ...
a
Readers,
I am having difficulty understanding how to enter commands into r.
I have data arranged as:
100, 200, 300
5.6, 6.7, 7.8
8.9, 9.0, 0.1
1.2, 2.3, 3.4
The data is saved in csv format and I use the command 'read.table' to
import into r.
The values 5.6...3.4 are a function of values 100,..
?rm
Richard
Ralf Goertz wrote:
How can I automatically exclude one variable from being saved in the
workspace when I quit an R session? The problem is I don't know how to
erase a variable once it has been created.
Background: I open a connection called "con" to a database server in my
~/.Rprof
Hi All,
I want to hide all the system(R) generated warning messages and want only
the messages from my script should be shown when i run myscript.R
Can you please suggest how can I hide the warnings
Thank you
Vidhu
[[alternative HTML version deleted]]
How can I automatically exclude one variable from being saved in the
workspace when I quit an R session? The problem is I don't know how to
erase a variable once it has been created.
Background: I open a connection called "con" to a database server in my
~/.Rprofile. Obviously, the connection expi
It appears that these are time series (i.e. there is only one value of
River Mile in any one data frame) so if that's right you would be better
off representing them as time series. Using zoo:
> library(zoo)
> feb06.z <- zoo(feb06[,2], feb06[,1])
> may06.z <- zoo(may06[,2], may06[,1])
> jun06.z <
Hi,
I'm plotting a bar chart like this:
ggplot() +
geom_bar(data=res,aes(fill=f1,x=f2,y=y),stat="identity",position="dodge")
f1 contains quite a few levels and the plot is really quite difficult to
read when the order of bars on the graph and on the legend does not match.
This problem has been
On Tue, Apr 29, 2008 at 03:53:02PM +0200, Julien Roux wrote:
> Hi list,
> I created a function to plot my data:
> plot_function(vector)
> I want to write the name of the argument vector in the legend/title of
> the plot.
> For example if I call:
> > plot_function(my_vector)
> I want "my_vector" to
Hello Andy,
thanks for your answer and sorry that I did not check the rfNews.
You are right, presently I am only looking at variable importance and
therefore I am very happy to hear that I won`t have problems using
ordered factors.
Since I just started to fiddle around with randomForest it mi
merge can only merge two objects at a time- I would like to merge more than
two objects at a time.
s.d <- structure(list(RiverMile = c(202L, 198L, 190L, 185L, 179L, 148L,
119L, 61L)), .Names = "RiverMile", row.names = c(NA, -8L), class =
"data.frame")
#s.d is all of the river miles that can occur
On 29/04/2008 9:53 AM, Julien Roux wrote:
Hi list,
I created a function to plot my data:
plot_function(vector)
I want to write the name of the argument vector in the legend/title of
the plot.
For example if I call:
> plot_function(my_vector)
I want "my_vector" to be written in the legend or ti
If you are using the latest version (4.5-25), you will see in rfNews() that
that's the problem I need to fix. The package was able to handle ordered
factors, but some more stringent checks for factor levels consistency
introduced in 4.5-23 broke the support for ordered factors in prediction.
On Tue, Apr 29, 2008 at 10:20 AM, Alberto Monteiro
<[EMAIL PROTECTED]> wrote:
>
> Gabor Grothendieck wrote:
> >
> > By write, do you mean print?
> >
> No, I mean "save to a file".
>
> I solved the problem with this:
>
> library(XML)
> x <- xmlTreeParse("Irpf2008/aplicacao/dados/12345678901/12345678
Gabor Grothendieck wrote:
>
> By write, do you mean print?
>
No, I mean "save to a file".
I solved the problem with this:
library(XML)
x <- xmlTreeParse("Irpf2008/aplicacao/dados/12345678901/12345678901.xml")
sink("ihatetheirs.xls")
print(x)
sink()
and then I can edit the saved file to cut some
s.d <- structure(list(RiverMile = c(202L, 198L, 190L, 185L, 179L, 148L,
119L, 61L)), .Names = "RiverMile", row.names = c(NA, -8L), class =
"data.frame")
#s.d is all of the river miles that can occur in all of the data frames that
I want to put together
feb06 <- structure(list(RiverMile = c(202L, 1
Alberto Monteiro wrote on 04/29/2008 08:37 AM:
Is there any function to write a XML structure, after it was
read using xmlTreeParse?
Ex:
library(XML)
x <- xmlTreeParse("Irpf2008/aplicacao/dados/12345678901/12345678901.xml")
# write it...
Calling:
library(help=XML) # or
help(package=XML)
li
By write, do you mean print?
In that case using the basic.xml file that comes with the XML package:
library(XML)
basic.xml <- system.file("exampleData", "basic.xml", package = "XML")
# try this
con <- xmlTreeParse(basic.xml)
root <- xmlRoot(con)
root
# or if you are using internal nodes:
con
You forgot to mention your version of R and your platform (and R has at
least 3 separate png() devices).
In 2.7.0, specify the size in inches and increase the resolution.
In earlier versions, try bitmap(). From NEWS for 2.7.0
o Considerable efforts have been made to make the default outp
Hi list,
I created a function to plot my data:
plot_function(vector)
I want to write the name of the argument vector in the legend/title of
the plot.
For example if I call:
> plot_function(my_vector)
I want "my_vector" to be written in the legend or title and so retrieve
the name of this object
Hi all,
I prepared a few charts with ggplot2, and was happy with the results as
displayed on screen.
I tried to draw them with a PNG driver instead. But I ran into several
problems while trying to increase the resolution of the picture.
Mainly, when I increased the picture size (e.g png(width
Is there any function to write a XML structure, after it was
read using xmlTreeParse?
Ex:
library(XML)
x <- xmlTreeParse("Irpf2008/aplicacao/dados/12345678901/12345678901.xml")
# write it...
Alberto Monteiro
PS: please, brazilians, don't be offended by my foul language!
___
Johannes Graumann wrote:
> Vioplots have great appeal to me, as they manage to squeeze so much
> information into so little space ...
> Now some evaluation has made me suspicious about the implementation in the
> package vioplot and I would like to hear what you say about the appended
> results
Dear R-users,
I have problems specifying a model with nested random terms. I want to
reproduce a variance component model made in Mintab, the model
specification in Minitab is:
Response = Prod P Prod*P Sample(Prod P) Extract(Prod P Sample)
Random Prod Sample Extract
The dataset is bala
Hi,
I would like to use a weighted lm model to reduce heteroscendasticity. I am
wondering if the only way to generate the weights in R is through the
laborious process of trial and error by hand. Does anyone know if R has a
function that would automatically generate the weights need for lm?
Thank
How about this
result <- numeric(10)
for(i in 1:10){
x <- sample(1:30, 5, replace = FALSE)
x <- sort(x)
result[i] <- any(diff(x) == 1)
}
> -Original Message-
> From: [EMAIL PROTECTED]
> [mailto:[EMAIL PROTECTED] On Behalf Of Anthony28
> Sent: Tuesday, April 29, 2008 8:52 AM
> To:
On Tue, Apr 29, 2008 at 01:49:01PM +0200, Philipp Pagel wrote:
>
> > When exporting to PDF a graph with a legend, in the final PDF, the
> > text is going beyond the legend box.
> > > dev2bitmap("test.pdf", type="pdfwrite", h=6, w=6)
> > The legend looks OK on the screen. I noticed that the size
1 - 100 of 149 matches
Mail list logo
