Neotropical bat risk assessments wrote:
>
>
>How do people deal with R and memory issues?
>I have tried using gc() to see how much memory is used at each step.
>Scanned Crawley R-Book and all other R books I have available and the
> FAQ
>on-line but no help really found.
>R
On Apr 26, 2009, at 8:40 PM, Ronnen Levinson wrote:
Hi.
Is there an analog to abline() that can be used to plot a polynomial
fit?
For example, I can draw the straight-line fit
fit <- lm(y ~ x)
via
abline(coef=fit$coef)
but I'm not sure how to draw the polynomial fit
fit <- lm
This is easy...
qqnorm(res, datax = TRUE)
... unfortunately.
Bill Venables
http://www.cmis.csiro.au/bill.venables/
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Chris_d
Sent: Monday, 27 April 2009 9:38 AM
To: r-help@r-projec
Nigel Birney wrote:
>
> Hello all,
>
> I have to import numeric data from file but found it contains Infinite
> values which need to be eliminated. I tried to replace them in this way:
>
> data[which(data=="-Inf")] <- -0.3
> data[which(data=="+Inf")] <- 0.3
>
> But, somehow, the
?qqnorm
qqnorm(y, ylim, main = "Normal Q-Q Plot",
xlab = "Theoretical Quantiles", ylab = "Sample Quantiles",
plot.it = TRUE, datax = TRUE, ...)
Charles Annis, P.E.
charles.an...@statisticalengineering.com
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.
Use ?is.infinite
inf <- is.infinite(data)
data[inf] <- 0.3 * sign(data[inf])
On Sun, Apr 26, 2009 at 5:44 PM, Nigel Birney wrote:
>
> Hello all,
>
> I have to import numeric data from file but found it contains Infinite
> values which need to be eliminated. I tried to replace them in this way:
>
?is.infinite
Charles Annis, P.E.
charles.an...@statisticalengineering.com
phone: 561-352-9699
eFax: 614-455-3265
http://www.StatisticalEngineering.com
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Nigel Birney
Sent: Sunday, Apr
Try (re)reading ?qqnorm. Use datax = TRUE.
--sundar
On Sun, Apr 26, 2009 at 4:37 PM, Chris_d wrote:
>
> Hi all,
> I have just started using R to produce qqnorm plots. I am trying to
> switch the x and y axes so that the theoretical values are plotted on the y
> axis and my data on the x ax
Hello all,
I have to import numeric data from file but found it contains Infinite
values which need to be eliminated. I tried to replace them in this way:
data[which(data=="-Inf")] <- -0.3
data[which(data=="+Inf")] <- 0.3
But, somehow, the Infinite values stayed there. Any sugg
Hi all,
I have just started using R to produce qqnorm plots. I am trying to
switch the x and y axes so that the theoretical values are plotted on the y
axis and my data on the x axis. Can anyone help me with this?
Thanks
--
View this message in context:
http://www.nabble.com/Flipping-axe
Hi all,
I am saving a program's output to a file to be read by another algorithm.
But somehow such a simple operation (the reading) fails. I get:
Error in read.table("a_corr_data.txt", sep = ",", col.names = T, row.names =
F) :
more columns than column names
Here is the write statement:
w
Am 22.Apr.2009 um 21:45 schrieb Duncan Temple Lang:
> Hi Benno et al.
>
> I have had some code for reading RDA files via R functions
> and binary connections. It is available from
>
> http://www.omegahat.org/RDA
>
> or
> install.packages("RDA", repos = "http://www.omegahat.org/R";,
>
For base graphics, you can do it with par - specifically
mfrow or mfcol
par(mfrow=c(3,3)).
Depending on exactly what you want, you might also
wish to look into axis().
For fancier set-ups, the lattice package provides many
options.
Sarah
On Sun, Apr 26, 2009 at 8:32 PM, Christian Bustamante wr
Try using RDCOMClient or rcom:
library(RDCOMClient)
xl <- COMCreate("Excel.Application")
# next line optional
xl[["Visible"]] <- TRUE
wb <- xl[["Workbooks"]]$Open("/mydir/sample.xlsx")
sheet <- wb$Sheets("Sheet2")
xl[["DisplayAlerts
Hi.
Is there an analog to abline() that can be used to plot a polynomial fit?
For example, I can draw the straight-line fit
fit <- lm(y ~ x)
via
abline(coef=fit$coef)
but I'm not sure how to draw the polynomial fit
fit <- lm(y ~ poly(x,2))
I do see the function curve(), but not
I all,
I want to create a graph layout in a 3x3 matrix like this:
ylab |__| |__| |__|
___ ___ ___
ylab |__| |__| |__|
___ ___ ___
ylab |__| |__| |__|
xl xl xl
With this layout, then I'll insert the 9 plots. How ca I create it?
--
CdeB
__
Here's a ggplot2 based solution:
#load the ggplot2 library
library(ggplot2)
#here's the data provided by Andreas
countries <- c("U-lande", "Afrika syd for sahara", "Europa og
Centralasien", "Lantinamerika og Caribien","Mellemøstenog
Nordafrika","Sydasien","ØStasien og stillehaveet", "Kina",
"Bras
The code below shows what I'm trying to get rid of.
If there is no way to get rid of the loop, I will try to use package( inline
).
I'm just curious as to whether there is a "vector way" of doing this
algorithm.
#
-
Hi!
I'm manipulating XLSX data using RODBC, however a limitation which appears
to be driver based is that you can't clear or drop sheets from the XLSX
files, as per the following example:
> library(RODBC)
> xlsx<-odbcDriverConnect("DRIVER=Microsoft Excel Driver (*.xls, *.xlsx,
*.xlsm, *.xlsb);DBQ
Have a look at plotweb in the bipartite package.
On Sun, Apr 26, 2009 at 6:45 PM, Andreas Christoffersen
wrote:
> Hi there,
>
> I would like to make a 'bumps chart' like the ones described e.g.
> here: http://junkcharts.typepad.com/junk_charts/bumps_chart/
>
> Purpose: I'd like to plot the propor
In statistics, a bumps chart is more commonly called a parallel
coordinates plot.
Hadley
On Sun, Apr 26, 2009 at 5:45 PM, Andreas Christoffersen
wrote:
> Hi there,
>
> I would like to make a 'bumps chart' like the ones described e.g.
> here: http://junkcharts.typepad.com/junk_charts/bumps_chart/
Have a look at the plyr package and associated documentation -
http://had.co.nz/plyr
Hadley
On Sun, Apr 26, 2009 at 12:42 PM, wrote:
> After a year my R programming style is still very "C like".
> I am still writing a lot of "for loops" and finding it difficult to recognize
> where, in place o
Hi there,
I would like to make a 'bumps chart' like the ones described e.g.
here: http://junkcharts.typepad.com/junk_charts/bumps_chart/
Purpose: I'd like to plot the proportion of people in select countries
living for less then one USD pr day in 1994 and 2004 respectively. I
have already constru
Thanks for your replies. I ended up using the following:
df = data.frame(year
=c(1991,1991,1992,1992,1993,1993,1992,1991),x=rnorm(8),y=rnorm(8))
df
year x y
1 1991 0.5565083 -1.31364232
2 1991 0.1686598 -0.20344656
3 1992 -0.1010090 -0.65681852
4 1992 0.6130324 -0.107
Agreed on the redirect to SIG. mpich2-devel and -libs are installed.
V.
On Sun, Apr 26, 2009 at 4:51 PM, Dirk Eddelbuettel wrote:
>
> On 26 April 2009 at 14:50, Vince Fulco wrote:
> | On FC10 with openmpi and mpich2 installed, the command "R CMD INSTALL
> | Rmpi_0.5-7.tar.gz --configure-args=--
On 26 April 2009 at 14:50, Vince Fulco wrote:
| On FC10 with openmpi and mpich2 installed, the command "R CMD INSTALL
| Rmpi_0.5-7.tar.gz --configure-args=--with-mpi=/usr/lib64/openmpi" (or
| /usr/lib64/mpich2) fails with the error ''cannot find mpi.h". Doing a
| (s)locate indicates no header fil
On 26 April 2009 at 13:25, Tom La Bone wrote:
| I just installed Ubuntu 9.04 but there does not seem to be repository for
| binaries for this version. Are there going to be such repositories set up in
| the near future?
As I understand it Michael and Vincent are working on it right now.
This wou
Tena koe Ken
Would something along the following lines do what you require:
set.seed(1)
x <- runif(100)
y <- rep(NA, length(x))
y[x<0.25] <- -1
y[x>0.75] <- 1
y[-1][y[-length(y)]%in%1 & (x[-1]>=0.25 & x[-1]<0.5)] <- 0
# etc
HTH
Peter Alspach
> -Original Message-
> From: r-help-bou
Hello,
my first attempt at installing version 2.9.0 failed
because I got an error
"Error in library(pspline) : there is no package called 'pspline'"
Later I realised that this comes from HOME/.RProfil, and removing
that files "solves" that problem.
However, I'm actually glad that this error happ
I just installed Ubuntu 9.04 but there does not seem to be repository for
binaries for this version. Are there going to be such repositories set up in
the near future?
Tom
--
View this message in context:
http://www.nabble.com/R-64-bit-for-Ubuntu-9.04-64-bit-tp23246229p23246229.html
Sent fro
I was trying to resist responding to this question since the original
questioner
had already been admonished twice last october about asking questions
on R-help about posted code that was not only not a part of R-base,
but not even a part of an R package. But the quoted comment about
Stata is
Hi there,
I have one problem calculating an linear mixed model.
I have a repeated measurement project for several gases. These were measured in
3 different altitudes (sites) and three different positions within these sites
(plot). each I now wanted to find out if the gasfluxes rather depend on
On 26/04/2009, at 3:56 PM, Rebecca1117 wrote:
I am new in R.
I can simulate Arma, using Arima.sim
However, I want to simulate an Arima Model. Say (1-B)Zt=5+(1-B)at.
I do not
know how to deal with 5 in this model.
Can any one could help me?
Thank you very much!
If this is a homework pro
Hi,
I'm trying to use Information Gain for feature selection.
There is a InfoGain implementation in Weka:
*weka.attributeSelection.InfoGainAttributeEval*
Is it possible to use this function with RWeka? If yes how?
list_Weka_interfaces doesn't show it and there is no make function for
AttributeE
On FC10 with openmpi and mpich2 installed, the command "R CMD INSTALL
Rmpi_0.5-7.tar.gz --configure-args=--with-mpi=/usr/lib64/openmpi" (or
/usr/lib64/mpich2) fails with the error ''cannot find mpi.h". Doing a
(s)locate indicates no header file labeled as such. Would appreciate
any trailheads.
TI
Hello,
I would like to generate ramified structures like plant root systems, river
networks or trees and save the generated structure as an image. Does anyone
knows if there is a way to do that with R?
Thank you in advance,
Talita
[[alternative HTML version deleted]]
__
On 4/26/09, Qifei Zhu wrote:
> Hi David,
>
> Thanks! It looks much better now. but is there any way to add (x,y)
> coordinates as labels to all the points in the graph? Best case if I can
> enforce some conditions saying if (y>10,000) label, else no label. Any
> advice is appreciated.
Sure, w
Also you don't need the abs since you are
squaring it anyways and seq with length
argument is a bit cleaner:
ix <- seq(Nstart, length = NW)
sum(Y[ix]^ 2)
Also please read the last line on every
message to r-help. The code in your
post lacks all 4 of the asked for
ingredients: 1. there are no co
drmh wrote
>
>Hello again,
>In my situation, I have three variables: pretest, posttest, and cohesion.
>
>I want to work out the correlation between postest and cohesion.
>
cor(cohesion, posttest) gives you this.
>I looked at multiple sets of data and created ANOVA tables of them. However,
>as
I think you can replace your 'for' loop with vectorized operations:
if(NW > 0){
Rnorm <- Rnorm + sum(abs(Y[Nstart:(Nstart + NW - 1)]) ^ 2)
}
On Sun, Apr 26, 2009 at 1:42 PM, wrote:
> After a year my R programming style is still very "C like".
> I am still writing a lot of "for loop
After a year my R programming style is still very "C like".
I am still writing a lot of "for loops" and finding it difficult to recognize
where, in place of loops, I could just do the
same with one line of code, using "sapply", "lapply", or the like.
On-line examples for such high level function d
hesicaia a écrit :
Hello,
I have a specific question regarding figure layout. I am tryng to make a 2
by 1 figure but I would like to make the bottom figure slightly larger than
the top figure. I have read through the help posts and have tried to use
"fig=c(),new=T" and have also tried to u
ONKELINX, Thierry wrote:
Have a look at all.equal
matA <- matrix(1:4, ncol = 2)
matB <- matA
all.equal(matA, matB)
matB[1,1] <- -10
all.equal(matA, matB)
Hi Thierry,
Thanks, all.equal does indicate if it's all equal so that
works great!
Much nicer than my hack - thanks,
Esmail
David Winsemius wrote:
> identical( matrix((1:4), ncol=2), matrix((1:4), nrow=2))
[1] TRUE
> identical( matrix((1:4), ncol=2), matrix((2:5), nrow=2))
[1] FALSE
Thanks David, much appreciated.
Esmail
__
R-help@r-project.org mailing list
https://s
baptiste auguie wrote:
I'm not sure I'm following you but have you tried,
identical(matrix(c(1,1,1,1),ncol=2), matrix(c(1,1,1,1),ncol=2))
?all.equal
?isTRUE
?identical
and possibly the compare package,
compare(matrix(c(1,1,1,1),ncol=2),matrix(c(1,1,1,1),ncol=2))
HTH,
baptiste
Hi Babtiste
On Apr 26, 2009, at 1:02 PM, Esmail wrote:
I'm trying to compare two matrices made up of bits.
doing a simple comparison of
matA == matB
> identical( matrix((1:4), ncol=2), matrix((1:4), nrow=2))
[1] TRUE
> identical( matrix((1:4), ncol=2), matrix((2:5), nrow=2))
[1] FALSE
yields
unfortunately this problem is difficult do reproduce. If no, it would be
mentioned and removed earlier.
It probably heppens in some specific cirumstances. I I look for that
Jarek
Uwe Ligges pisze:
Jarek Jasiewicz wrote:
I tried several mirrors
But, what may be more important.
This error w
I'm not sure I'm following you but have you tried,
identical(matrix(c(1,1,1,1),ncol=2), matrix(c(1,1,1,1),ncol=2))
?all.equal
?isTRUE
?identical
and possibly the compare package,
compare(matrix(c(1,1,1,1),ncol=2),matrix(c(1,1,1,1),ncol=2))
HTH,
baptiste
On 26 Apr 2009, at 18:02, Esmail wro
Peter Dalgaard wrote:
> Nigel Birney wrote:
>> Hello all,
>>
>> Apologize for the newbie question. What's the easiest way to do a SQL
>> inner
>> table join in R?
>> Say I have a table containing column names A, B, C and another which has
>> columns named C, D, E. I would like to do an inner table
Have a look at all.equal
matA <- matrix(1:4, ncol = 2)
matB <- matA
all.equal(matA, matB)
matB[1,1] <- -10
all.equal(matA, matB)
HTH,
Thierry
ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek / Research In
Then post the material that would make sense for Windows.
What _does_ memory.limits() return? This _was_ asked and you did not
answer.
How many other objects do you have in your workspace?
How big are they?
Jim Holtman offered this function that displays memory occupation by
object and total
Hello again,
In my situation, I have three variables: pretest, posttest, and cohesion.
I want to work out the correlation between postest and cohesion.
I looked at multiple sets of data and created ANOVA tables of them. However,
as pretest and postest are sometimes correlated (with a statistic
I have a small database (file csv):
'District';'HouseType';'Income';'PreviousCustomer';'Outcome'
'Suburban';'Detached';'High';'No';'Nothing'
'Suburban';'Detached';'High';'Yes';'Nothing'
'Rural';'Detached';'High';'No';'Responded' ... itd.
After instruction str() as a result I get:
> str(ziel
i am using vglm for multiple logistic regression.
i have 1 response variable (total 4 category)
and 5 predictor.
Call:
vglm(formula = class ~ PC1 + PC2 + PC3 + PC4 + PC5, family = multinomial(),
na.action = na.pass)
Coefficients:
(Intercept):1 (Intercept):2 PC1:1 PC1:2
Hello,
I have a specific question regarding figure layout. I am tryng to make a 2
by 1 figure but I would like to make the bottom figure slightly larger than
the top figure. I have read through the help posts and have tried to use
"fig=c(),new=T" and have also tried to use "split.screen" and
I'm trying to compare two matrices made up of bits.
doing a simple comparison of
matA == matB
yields this sort of output.
[,1] [,2] [,3] [,4] [,5] [,6]
[1,] FALSE TRUE FALSE TRUE TRUE FALSE
[2,] TRUE TRUE TRUE TRUE TRUE TRUE
[3,] FALSE TRUE FALSE FALSE FALSE TRUE
[4,]
Neotropical bat risk assessments wrote:
Seems that the latest version R version 2.9.0 (2009-04-17)
has a glitch and will not install packages. Issue with unzipping?
Works fine with R version 2.8.1 (2008-12-22)
> install.packages("ff")
--- Please select a CRAN mirror for use in this session
Seems that the latest version R version 2.9.0 (2009-04-17)
has a glitch and will not install packages. Issue with unzipping?
Works fine with R version 2.8.1 (2008-12-22)
> install.packages("ff")
--- Please select a CRAN mirror for use in this session ---
trying URL 'http://cran.fhcrc.org/bin/w
On Sun, 26 Apr 2009, Uwe Ligges wrote:
Works for me with
library(MASS)
plot(lda(Species~., data=iris))
hence you may want to profide the data to enable us to reproduce your
problem...
He is trying to plot the results from a cross-validation. As the help
page clearly states, that is a list
Dear R users,
I am trying to do exact matching on a large dataset (500.000 obs), about equal
size of treatment and controll group, with replacement: As for the moment I use
the "Match" function of the "Matching" library. I match on 2 covariates and all
observations in the treatment group have a
On Sun, 26 Apr 2009 09:20:12 -0600 Neotropical bat risk assessments
wrote:
NBRA>
NBRA>How do people deal with R and memory issues?
NBRA>I have tried using gc() to see how much memory is used at each
NBRA> step. Scanned Crawley R-Book and all other R books I have
NBRA> available and the F
On Apr 26, 2009, at 11:20 AM, Neotropical bat risk assessments wrote:
How do people deal with R and memory issues?
They should read the R-FAQ and the Windows FAQ as you say you have.
http://cran.r-project.org/bin/windows/base/rw-FAQ.html#There-seems-to-be-a-limit-on-the-memory-it-uses_002
Jarek Jasiewicz wrote:
I tried several mirrors
But, what may be more important.
This error was mentioned earlier on VISTA and WindowsXP, I use Ubuntu 8.04
Unfortunately, I cannot reproduce this for "sp" on Windows XP.
Uwe
Uwe Ligges pisze:
Jarek Jasiewicz wrote:
since 2.9.0 version
Hi Tim,
There are two main problems with your implementation of Stochastic gradient
algorithm:
1. You are only implementing one cycle of the algorithm, i.e. it cycles over
each data point only once. You need to do this several time, until convergence
of parameters is obtained.
2. Stochastic
David Winsemius
wrote:
My understanding of the OP's request was for some sort of copy which did
change when entries in the original were changed; the sort of behavior that
might be seen in a spreadsheet that had a copy "by reference".
You misunderstood (my phrasing wasn't probably the best)
In that case, you would want a shallow copy, and you'd need to jump
through a lot of hoops to do that in R.
Hadley
On Sun, Apr 26, 2009 at 10:35 AM, David Winsemius
wrote:
> My understanding of the OP's request was for some sort of copy which did
> change when entries in the original were change
hadley wickham wrote:
I want to (1) create a deep copy of pop,
I have already said *I* do not know how to create a "deep copy" in R.
Creating a deep copy is easy, because all copies are "deep" copies.
You need to try very hard to create a reference in R.
Hi Hadley
Right you are .. I discove
David,
Good news! It seems that R has deep copy by default. I ran this simplified
test and it seems I can change 'pop' without changing the saved version.
POP_SIZE = 4
LEN = 8
pop=create_pop_2(POP_SIZE, LEN)
cat('printing original pop\n')
print(pop)
keep_pop = pop
pop[1,1] = 99
cat('printing c
My understanding of the OP's request was for some sort of copy which
did change when entries in the original were changed; the sort of
behavior that might be seen in a spreadsheet that had a copy "by
reference".
On Apr 26, 2009, at 11:28 AM, hadley wickham wrote:
I want to (1) create a d
>> I want to (1) create a deep copy of pop,
>
> I have already said *I* do not know how to create a "deep copy" in R.
Creating a deep copy is easy, because all copies are "deep" copies.
You need to try very hard to create a reference in R.
Hadley
--
http://had.co.nz/
_
How do people deal with R and memory issues?
I have tried using gc() to see how much memory is used at each step.
Scanned Crawley R-Book and all other R books I have available and the FAQ
on-line but no help really found.
Running WinXP Pro (32 bit) with 4 GB RAM.
One SATA drive p
I tried several mirrors
But, what may be more important.
This error was mentioned earlier on VISTA and WindowsXP, I use Ubuntu 8.04
Uwe Ligges pisze:
Jarek Jasiewicz wrote:
since 2.9.0 version I have a problem with installing packages:
install.packages("sp")
--- Please select a CRAN mirror
On Apr 26, 2009, at 9:43 AM, Esmail wrote:
David Winsemius wrote:
Yes. As I said before "I am going to refrain from posting
speculation until you provide valid R code
that will create an object that can be the subject of operations."
The code I have provided works, here is a run that may
I tried several mirrors
Uwe Ligges pisze:
Jarek Jasiewicz wrote:
since 2.9.0 version I have a problem with installing packages:
install.packages("sp")
--- Please select a CRAN mirror for use in this session ---
Loading Tcl/Tk interface ... done
Warning: unable to access index for repository
Ken-JP wrote:
set.seed(1)
x <- runif(100)
# I want to calculate y such that:
#
# 1. if x>0.75, y <- 1
# 2. else if x<0.25, y <- -1
# 3. else if y_prev==1 && x<0.5, y <- 0
# 4. else if y_prev==-1 && x>0.5, y <- 0
# 5. else y <- y_prev
#
# 1. and 2. are directly doable without looping.
#
# How d
Works for me with
library(MASS)
plot(lda(Species~., data=iris))
hence you may want to profide the data to enable us to reproduce your
problem...
Uwe Ligges
pgseye wrote:
Hi,
I've performed an lda and obtained a classification table for some of my
data:
efa.dfa<-lda(groups~.,efa.scores.
WilDsc0p wrote:
Dear List,
I have a question regarding "ipred" package. Under 10-fold cv, for different
knn ( = 1,3,...25), I am getting same misclassification errors:
#
library(ipred)
data(iris)
cv.k = 10 ## 10-fold cross-validation
bwpredict.kn
Caroline Wallis wrote:
Hi,
I'm using e1071 package to do fuzzy cluster analysis. My dataset (ra) has
5237 observations and 2 variables - depth and velocity. I used fuzzy cmeans
to create 6 fuzzy classes.
ra.flcust6<-cmeans(ra,6,iter.max=100,verbose=F,dist="euclidean",method="cmea
ns",m=
Jarek Jasiewicz wrote:
since 2.9.0 version I have a problem with installing packages:
install.packages("sp")
--- Please select a CRAN mirror for use in this session ---
Loading Tcl/Tk interface ... done
Warning: unable to access index for repository
http://piotrkosoft.net/pub/mirrors/CRAN/src
Dear Karin,
If I understand correctly what you want, the scatterplot function in the car
package isn't designed to produce it, but there are many ways to draw
side-by-side scatterplots. Here is one, using basic R graphics:
par(mfrow=c(1,2))
by(Data, Data$group,
function(x) {
plot(Pul
Thank you.
Unluckily what makes the problem only apparenttly simple (for me) is that we
have not differentiable functions and the parameter space is not continuous ...
which reduces dramatically the number of choices.
I would be grateful to chat with anyone who has tackled a similar problem.
Ma
On Apr 26, 2009, at 8:02 AM, Alessandro wrote:
Dear User,
thank for the attention. I have a data.frame with 5 columns (ex:ID,
a1,a2,a3,a4) and 1000 rows. I wish to find the absolute max value
for all
data.frame and save a new data.frame with the row where is that
value. Ex:
ID: 1,2,3,4,5,
David Winsemius wrote:
Yes. As I said before "I am going to refrain from posting speculation
until you provide valid R code
that will create an object that can be the subject of operations."
The code I have provided works, here is a run that may prove helpful:
POP_SIZE = 6
LEN = 8
pop=cre
On Sun, 26 Apr 2009 09:29:39 + (GMT) "nonu...@yahoo.de"
wrote:
ND> I'm realy new to R, so I hope you can help me, as I didn't find any
ND> solution in the common books.
Have a look at the online resources at:
http://cran.r-project.org/other-docs.html
There is also stuff on graphics.
Furth
Daryl Morris wrote:
Thanks for the various responses.
It was easier than I thought to get all the tools together and setup
Windows paths for the build process. I've been successful now on Windows.
BUT... how do I make a package which DOES NOT require R CMD INSTALL to
install? Obviously,
since 2.9.0 version I have a problem with installing packages:
install.packages("sp")
--- Please select a CRAN mirror for use in this session ---
Loading Tcl/Tk interface ... done
Warning: unable to access index for repository
http://piotrkosoft.net/pub/mirrors/CRAN/src/contrib
Warning messages
On Apr 26, 2009, at 7:48 AM, Esmail wrote:
Hello David,
Let me try again, I don't think this was the best post ever I've
made :-)
Hopefully this is clearer, or otherwise I may break this up into
three separate simple queries as this may be too long.
> "==" is not an assignment operator in
David Winsemius wrote:
On Apr 25, 2009, at 6:57 PM, reneepark wrote:
Hello,
I am using Dr. Harrell's design package to make a nomogram. I was able to
make a beautiful one without stratifying, however, I will need to
stratify
to meet PH assumptions. This is where I go wrong, but I'm not sure
Dear Alessandro,
Here is one way:
DF <- data.frame(ID,a1,a2,a3,a4)
Row <- which( DF == max(DF[,-1]), arr.ind = TRUE)[1]
DF[Row,]
# ID a1 a2 a3 a4
# 10 10 10 20 30 40
See ?which and ?max for more details.
HTH,
Jorge
On Sun, Apr 26, 2009 at 8:02 AM, Alessandro wrote:
> Dear User,
>
>
>
> t
Hi David,
Thanks! It looks much better now. but is there any way to add (x,y)
coordinates as labels to all the points in the graph? Best case if I can
enforce some conditions saying if (y>10,000) label, else no label. Any
advice is appreciated.
Best,
Tony
-Original Message-
From: David W
I want to generate some orthogonal tables in my experiment design.
I searched http://search.r-project.org/ (RSiteSearch) whith key words such as
"orthogonal table", "orthogonal design", "latin square" etc and get no useful
result. the same result get by searching via google's insite search in r-
I have various objects defined but I am trying to remove a set of elements in
one vector from another and using the loops at the end of this post but I am
getting the error at the very end of this post.
> str(x)
num [1:923, 1:923] 1 -0.00371 -0.00102 -0.00204 -0.00102 ...
- attr(*, "dimnames")=
Dear User,
thank for the attention. I have a data.frame with 5 columns (ex:ID,
a1,a2,a3,a4) and 1000 rows. I wish to find the absolute max value for all
data.frame and save a new data.frame with the row where is that value. Ex:
ID: 1,2,3,4,5,6,7,8,9,10
a1:1,2,3,4,5,6,7,8,9,10
a2:11,12,13,
Hello David,
Let me try again, I don't think this was the best post ever I've made :-)
Hopefully this is clearer, or otherwise I may break this up into
three separate simple queries as this may be too long.
> "==" is not an assignment operator in R, so the answer is that it
> would do neither.
Hi,
You could do this very easily using ggplot2,
#install.packages("ggplot2", dep=TRUE)
library(ggplot2)
c <- ggplot(mtcars, aes(y=wt, x=mpg)) + facet_grid(. ~ cyl)
c + stat_smooth(method=lm) + geom_point()
See more examples on Hadley's website: http://had.co.nz/ggplot2/
Hope this helps
Hi,
Have you considered using high-level plotting functions provided by
the ggplot2 or lattice package? Here's a dummy example,
x <- seq(0, 10, length=100)
y1 <- sin(x)
y2 <- cos(x)
y3 <- x^2/100
y4 <- 1/x
d <- data.frame(x, y1, y2, y3, y4)
library(reshape)
dm <- melt(d, id="x")
dm$type1
This is a nontrivial problem. This comes up often on the Statalist
(-qreg- is for cross-section quantile regression):
"
You want to fit a plane through the origin using the L-1 norm.
This is not as easy as with L-2 norm (LS), as it is more
than a matter of dropping a constant predictor yet otherwi
I'm sorry - I meant a "median survival" estimate, not a median "risk."
I see - I didn't realize that by stratifying it would pool the levels of the
stratified variable. Hm, that is unfortunate, considering the stratified
variable is one that I would like to keep in the nomogram.
Thank you for yo
Nigel Birney wrote:
Hello all,
Apologize for the newbie question. What's the easiest way to do a SQL inner
table join in R?
Say I have a table containing column names A, B, C and another which has
columns named C, D, E. I would like to do an inner table join on C and
produce a table A, B, C,
Hi. guys,
I am trying to write my own Stochastic Gradient Ascent for logistic
regression in R. But it seems that I am having convergence problem.
Am I doing anything wrong, or just the data is off?
Here is my code in R -
lbw <-
read.table("http://www.biostat.jhsph.edu/~ririzarr/Teaching/7
1 - 100 of 114 matches
Mail list logo
