On Mon, Mar 4, 2013 at 1:39 PM, Tammy Ma metal_lical...@live.com wrote:
HI,
Assume I have the data frame generated from R as the following:
Product Price market_share
A 10010%
B 1109%
C 12020%
D 90
I appreciate the help and suggestions. I was afraid that this question would
be considered off topic, but thought I would give it a try to see if anyone
else gets these results from R output files. As I do not know what a hex
editor or hexbin are I guess I will not be able to continue this
I would try write.csv or write.table with row.names = FALSE - just a guess,
and do follow up on R-Sig-Geo if need be.
Cheers, Mike
On Tuesday, March 5, 2013, Kerry wrote:
I appreciate the help and suggestions. I was afraid that this question
would be considered off topic, but thought I would
Hi,
I have two different imaging modalities (for the identification of areas of
infarcted myocardium) that I need to compare regarding agreement and
consistency.
However, I don't think that methods like Cohen's Kappa, PCC, Bland-Altmann
or ICC are sufficient here as there is not just a pairwise
You might have a look at:
http://www.portfolioprobe.com/2012/12/17/a-look-at-historical-value-at-risk/
which points to a function for historical VaR.
As Nello said, we really need to know what it is
that you think doesn't work, before we can help
you with what you have.
It probably doesn't
Le lundi 04 mars 2013 à 20:28 +, Lopez, Dan a écrit :
Hi,
We have comment questions from a survey that we need to categorize.
What package and functions can I use in R to help do this?
If you are asking for some kind of clustering based on vocabulary, then
have a look at the tm package
Hi everybody!
Does anyone know a good way to color my images so that
when I print them out on a non-color-printer the colors used
can be distinguished well? As I have many categories I would
not want to assign the colors c(black, grey, white) by
hand.
Thank you!
[[alternative HTML
Have a look at bpy.colors()
Best regards,
ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek / Research Institute for Nature and
Forest
team Biometrie Kwaliteitszorg / team Biometrics Quality Assurance
Kliniekstraat 25
1070 Anderlecht
Belgium
+ 32 2 525 02 51
+ 32 54 43 61 85
On 03/05/2013 02:42 PM, Jyoti Sharma wrote:
hello sir
myself Jyoti Sharma, and i am working as preoject fellow in IGIB Delhi.
I need your help to know how to create a package as well as how to
post that package to the CRAN mirror, for public use. i have used
package.skeleton() command but its
That is bpy.colors() in the contributed sp package on CRAN, not in R itself.
On Tuesday, March 5, 2013, ONKELINX, Thierry wrote:
Have a look at bpy.colors()
Best regards,
ir. Thierry Onkelinx
Instituut voor natuur- en bosonderzoek / Research Institute for Nature and
Forest
team Biometrie
On Tue, 5 Mar 2013, David Studer wrote:
Hi everybody!
Does anyone know a good way to color my images so that
when I print them out on a non-color-printer the colors used
can be distinguished well? As I have many categories I would
not want to assign the colors c(black, grey, white) by
hand.
Hello,
Another document that could help is Friedrich Leisch, Creating R
Packages: A Tutorial, that can be found at
http://cran.r-project.org/doc/contrib/Leisch-CreatingPackages.pdf
Hope this helps,
Rui Barradas
Em 05-03-2013 11:15, Jim Lemon escreveu:
On 03/05/2013 02:42 PM, Jyoti Sharma
Dear all,
is it possible to print angle brackets (LaTeX notation: \langle, \rangle)? I
found that lceil and lfloor
are available, see demo(plotmath). But langle and rangle are not.
I tried to print utf8 characters directly as well without success.
I have also read something about a tikzDevice
1. A censored observation
2. It does not relate to either
3. See ?print.survfit . Recall also that the mean of a positive random
variable is the integral from 0 to infinity of the survival function. The
truncated mean is the integral from 0 to tau (819 in your case) of the survival
function.
On 3/5/2013 12:51 AM, Ingo Reinhold wrote:
Hi,
I'm trying to use the data which I generate within R to make images in .bmp
format to be lateron printed by a printer.
My first thought was the RImageJ package, but this seems to be discontinued. What I am
currently doing is generating a matrix
Cross-posted, verbatim, on stackoverflow:
http://stackoverflow.com/q/15203347/271616
--
Joshua Ulrich | about.me/joshuaulrich
FOSS Trading | www.fosstrading.com
R/Finance 2013: Applied Finance with R | www.RinFinance.com
On Mon, Mar 4, 2013 at 7:07 AM, Аскар Нысанов nysanas...@mail.ru
The question is in the title.
Then, I would like to know how I should refer to the documentation
regarding the use of each functions.
Thanks,
Julien Mehl Vettori
[[alternative HTML version deleted]]
__
R-help@r-project.org mailing list
On 11/3/2011 3:30 PM, Brian Diggs wrote:
Well, I figured it out. Or at least got it working. I had to run
initexmf --mkmaps
because apparently there was something wrong with my font mappings.
I
don't know why; I don't know how. But it works now. I think
installing
the font into the
Dear All,
I have a table as following
position type count
1 2 100
1 3 51
1 5 64
1 8 81
1 6 32
2 2 41
2 3 85
and so on
Normally if would have a vector of 2,3,4,5... by position position and
Dear Julien,
Check
citation('stats')
HTH,
Jorge.-
On Wed, Mar 6, 2013 at 12:05 AM, Julien Mvdb julien.m...@gmail.com wrote:
The question is in the title.
Then, I would like to know how I should refer to the documentation
regarding the use of each functions.
Thanks,
Julien Mehl
Julien,
I would just try your best given the journal's style guide and wait for
them to change it. For what it is worth, my last paper was corrected as
follows.
(In text)...one-way analysis of variance (ANOVA) in R [version 2.14.2
(R Core Team, 2012)].
(reference section)
R Core Team.
Does anyone know of any literature on the kappa statistic with plsda? I
have been trying to find papers that used plsda for classification and have
yet to come across this kappa value. All the papers I come across
typically have R2 as an indicator of model fit. I want to make sure I
conduct
If you wanted to do a t.test
res1-do.call(cbind,lapply(seq_len(nrow(m)),function(i)
do.call(rbind,lapply(split(rbind(m[i,-1],n),1:nrow(rbind(m[i,-1],n))),
function(x) {x1- rbind(x,m[i,-1]); t.test(x1[1,],x1[2,])$p.value}
res2-do.call(cbind,lapply(seq_len(ncol(res1)),function(i)
On 03/05/2013 05:00 AM, r-help-requ...@r-project.org wrote:
Hello, I create a plot from a coxph object called fit.ads4: plot(survfit(fit.ads4))
Questions: 1. What is the cross mark in the plot ? 2. How does the cross mark in the
plot relate to either the rmean or the median from survfit ?
This request is very completely satisfied by RExcel. Please look at
rcom.univie.ac.at for full
information including download information. Follow up should be on the
rcom email list.
You can sign up at the website.
On Mon, Mar 4, 2013 at 8:39 AM, Tammy Ma metal_lical...@live.com wrote:
HI,
I'm working on a survival analysis project and, based on previous work of
my colleagues on similar projects, am trying to estimate an AFT model with
a generalized gamma distribution. When I execute the following code:
flexsurvreg(timesurv ~ cycles + rcycle + cyctime, dist=gengamma)
I get the
This is my first time to use R.
For clarification:
I made an oa design for 4 factors with the following levels: one with 4 levels,
2 with 3 levels, 1 with 6 levels.
Using DoE package, I have generated 72 runs (setting the columns=min3).
However, the numbers of generalized words of lengths 3
Dear List,
I have the code below adapted from the lattice-package examples to draw two
spheres. I would now like to give
both different surface colors, e.g. one red and one blue.
## 3-D surface parametrized on a 2-D grid
n -
Hi,
I'm having problems with the Zelig package - when using
the below R displays the follwing message (I'm running R
i386 2.15.3 for Windows and have updated all the Zelig
packages):
z.out-zelig(Surv(psurv2, pcens2) ~ ren_sup3 + age,
data=data_urgent, model=coxph)
** The
To who(m) it may concern:
I'm so beginner with R , I have a data set and want to compare the Generalized
gamma distribution with different
K (gamma shape parameter) for my data, in other hand wanna to fit the gengamma
distribution when K has different value for example k=1, 2, 5, 10 ,…
would
I've updated to the latest version of R, but still the problem persists.
Another thing I noticed (but failed to mention before) is that when I
initially open the change working directory dialog box, the little line
labeled Folder: (under the window showing the folder tree) does show the
current
On 05/03/2013 10:22 AM, Glenn Stauffer wrote:
I've updated to the latest version of R,
That's a little ambiguous. Do you mean 2.15.3? I didn't put the new
change into that, because it was too late. You can try a build of the
development version, which will be released as 3.0.0 in April,
I'm sorry - yes I meant 2.15.3. But I think it probably does not matter.
Others are not seeing the behavior that I am (with any version), so I don't
think R is the problem. But I will post if I ever get it figured out.
Thanks,
Glenn Stauffer
-Original Message-
From: Duncan Murdoch
On 05/03/2013 12:05, Donatella Quagli wrote:
Dear all,
is it possible to print angle brackets (LaTeX notation: \langle, \rangle)? I
found that lceil and lfloor
are available, see demo(plotmath). But langle and rangle are not.
I tried to print utf8 characters directly as well without success.
Hi there:
I'm having a weird problem with my startup procedure. R.app is reading an
unknown .Rprofile file.
First, I'm on a Mac Os 10.6.8 running R.app 2.15.0
On startup
getwd()
[1] /Users/simon
But: the contents of my .Rprofile file in my home directory when viewed with a
text editor are:
Preamble: There is a dedicated R-Mac-SIG list where this would have been more
appropriately directed. and followups should be directed there (and the R-help
address removed so we do not have two simulatnaous threads.)
On Mar 5, 2013, at 8:38 AM, Simon Kiss wrote:
Hi there:
I'm having a weird
I'm working with capscale and permutest for the first time, and having
trouble getting statistical analyses for more than one constraining
variable. I've read the documentation, but setting first=FALSE or using
by=axis doesn't seem to be helping. capscale seems to be fine, I receive
output for
On R 2.15.2 and ArcGIS 9.3.1, it works for me in ArcCatalog but you have to
follow the particulars here:
http://webhelp.esri.com/arcgisdesktop/9.3/index.cfm?TopicName=Accessing_delimited_text_file_data
For example:
write.table(test, '***.tab', sep = '\t', row.names = F)
The extension .tab and
Dear R experts,
I am trying to simulate correlated binary data and have stumbled upon the
following problem:
With the help of binarySimCLF or mvpBinaryEp I have been able to
simulate correlating binary vectors given certain mean values and a desired
correlation. My problem is that these
Dears, a simple question here: is there any AR1 and AR2 test for FGLS-FE
fitted with the pggls function of plm package?
(one example would be the Baltagi-Wu LBI test)
Thanks for your attention!
[[alternative HTML version deleted]]
__
We are happy to inform you that registration for useR! 2013 is now open,
see
http://www.R-project.org/useR-2013
This meeting of the R user community will take place at the University
of Castilla-La Mancha, Albacete, Spain, July 10-12, 2013. Pre-conference
tutorials will be offered on July 9.
I think Zelig uses model=cox.ph (
https://github.com/zeligdev/ZeligMisc/blob/master/tests/coxph.R)
-thomas
On Wed, Mar 6, 2013 at 2:54 AM, Stephen Knight stephenkni...@doctors.org.uk
wrote:
Hi,
I'm having problems with the Zelig package - when using the below R
displays the follwing
For anyone looking for an intermediate to advanced level short course on R
this summer, I will be presenting at the BYU Summer Institute of Applied
Statistics June 19-21, 2013.
The official page is here:
http://statistics.byu.edu/r-beyond-basics-38th-annual-summer-institute-applied-statistics
Hi
Using the ggplot2 package, I would like to obtain a plot that contains two
time series that have data points on different dates.
For instance, one data frame looks like:
date1, value1
2010-01-05, 2921.74
2010-01-08, 2703.89
2010-01-14, 3594.21
2010-01-20, 3659.22
The other data frame looks
Dear all,
I got an error message when running the following code.
Can anyone give any suggestions on fixing this type of error?
Thank you very much in advance.
Hanna
integrand - function(x, rho, a, b, z){
+ x1 - x[1]
+ x2 - x[2]
+ Sigma -
On Tue, Mar 5, 2013 at 3:30 PM, Matthijs Daelman
matthijs.dael...@gmail.com wrote:
Hi
Using the ggplot2 package, I would like to obtain a plot that contains two
time series that have data points on different dates.
For instance, one data frame looks like:
date1, value1
2010-01-05, 2921.74
On Mar 5, 2013, at 3:28 AM, David Schellenberger Costa wrote:
Dear List,
I have the code below adapted from the lattice-package examples to draw two
spheres. I would now like to give
both different surface colors, e.g. one red and one blue.
## 3-D surface parametrized on a 2-D grid
On Tue, Mar 5, 2013 at 7:07 PM, Marbles max-ihm...@gmx.de wrote:
Dear R experts,
I am trying to simulate correlated binary data and have stumbled upon the
following problem:
With the help of binarySimCLF or mvpBinaryEp I have been able to
simulate correlating binary vectors given certain
Dear r-help,
Somewhere in my innocuous function to rotate an object in Cartesian space
I've created a monster that completely locks up my computer (requires a
hard reset every time). I don't know if this is useful description to
anyone - the mouse still responds, but not the keyboard and not
Hi,
I need to do some analysis on historic daily radiosonde data I download
from the Wyoming Univ. web page (
http://weather.uwyo.edu/upperair/sounding.html).
I am trying to use the RadioSonde package (V 1.3), but the format of the
files from Wyoming don't match what RadioSonde is expecting.
Has
Hi,
I’m attempting to export data (split into multiple files from one large
dataset) from R to excel using the excel.link package. The code for export
is as follows:
for(i in practicesNN){
#Create relevant data for input
#Separate out all parts of data – PracticeName is removed from example
Thank you for your response!
I know that it is not enough to uniquely specify the correlation. That is
why i would like to simulate it so that i can see how the resulting
correlations between A and C are distributed in dependence of the
correlations between AB and BC.
--
View this message in
Btw: I have done this for non-binary variables already. If anyone is
interested, it looks like this.
I am pretty new to R, so excuse the potentially unelegant code.
# Implement library
library(ecodist)
# Prepare ACOutcome vector
ACOutcomes = c()
# Set desired variablelength Number of
Hi there,
I am trying to write a tool which involves implementing logistic
regression. With the batch gradient descent method, the convergence is
guaranteed as it is a convex problem. However, I find that with the
stochastic gradient decent method, it typically converges to some random
points
Tena koe Benjamin
I haven't looked at you code in detail, but in general ifelse is slow and can
generally be avoided. For example,
ben - 1:10^7
system.time(BEN - ifelse(ben10, NA, -ben))
user system elapsed
1.310.241.56
system.time({BEN1 - -ben; BEN1[BEN1 -10] - NA})
user
Sure looks as if that second --... line is causing
another attempt at parsing for varkey and unitkey which for 70.0 and 59.6
just won't make sense to getsonde.
It's a pain but I'd experiment by removing that second ... line.
Clint BowmanINTERNET:
I see you have profiling calls in there. Have you used them?
It is often fruitful to see how the time for a function grows as size of the
input or output grows. Have you tried that?
A concrete suggestion is to change
for(i in 1:n){
if(is.na(or.deg[i])==TRUE) {or.deg[i] - 0}
}
to
Hi there!
Today I tried to estimate models using both plm and pgmm functions, with an
interaction between X1 and lag(X2, 1). And I notice two issues.
Let Y=b_1 * X_1 + b_2 * X_2 + b_3 * X_1 * x_2 + e be our model.
1) When using plm, I got different results when I coded the interaction
term with
HI Elisa,
Just noticed the order of elements in vec1:
You have to replace `vec1`
dat2- as.dist(dat1[,-1],upper=F,diag=F)
vec1- as.vector(dat2)
head(vec1)
#[1] 5.796656 43.523023 38.193750 44.730182 6.511703 2.904954 #the order is
based on columns
#with
dat1-
Hi,
I guess this is what you wanted.
Attaching a plot from a subset (lstSub)
dat1- read.csv(rightest.csv,sep=,,header=TRUE,check.names=FALSE)
label1=c(0-25,25-50,50-75)
Name1-unlist(lapply(0:123,function(i) rep(i+1,i)))
dat1New- dat1[,-1]
vec1-
A potentially ridiculous question, but why does R use wilcox (e.g.,
pwilcox or wilcox.test) instead of the full name Wilcoxon? I've
browsed (but not scoured) the help files and Peter Dalgaard's book,
but I'm coming up empty.
Purely for brevity or have I missed something massive?
## Reproducible
I may be missing something, but what does this have to do specifically
with R? I believe this is OT here and you need to post elsewhere, e.g.
perhaps on stats.stackexchange.com.
-- Bert
On Tue, Mar 5, 2013 at 1:36 PM, Ivan Li machinelearning2...@gmail.com wrote:
Hi there,
I am trying to write
Perhaps I should have added and FWIW, R, like essentially all
statistical software, has logistic regression already built in, if I
understand what you mean by the term (which I may not), via glm's.
-- Bert
On Tue, Mar 5, 2013 at 1:36 PM, Ivan Li machinelearning2...@gmail.com wrote:
Hi there,
Hi all,
Thanks for the suggestions. Updating the function as below to break the
problem into chunks seemed to do the trick - perhaps there is a relatively
small limit to the size of a vector that R can work with?
Best
rotate - function(x,y,tilt,threshold){
df.main-data.frame(x,y)
Actually, the problem number 2 is easy to solve: instead of using I(X1 *
lag(X2,1)), one should use X1:lag(X2,1). It works.
The issue number 1 remains, though. And also affects this solution for
number 2. It means: results are different with one uses X1:lag(X2,1)
whithin the formula or uses a new
Indeed! Thank you!
__
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.
On 2013-03-05 14:51, R. Michael Weylandt wrote:
A potentially ridiculous question, but why does R use wilcox (e.g.,
pwilcox or wilcox.test) instead of the full name Wilcoxon? I've
browsed (but not scoured) the help files and Peter Dalgaard's book,
but I'm coming up empty.
Purely for brevity or
On Mar 5, 2013, at 5:54 AM, Stephen Knight wrote:
Hi,
I'm having problems with the Zelig package - when using the below R displays
the follwing message (I'm running R i386 2.15.3 for Windows and have updated
all the Zelig packages):
z.out-zelig(Surv(psurv2, pcens2) ~ ren_sup3 + age,
Hello everone,
Anyone who knows how to force a cubic smoothing spline to pass
through a particular point?
I found on website someone said that we can use cobs package to
force the spline pass through certain points or impose shape
constraints (increasing,
Hello everyone,
I have a basic question regarding logical operators.
x-seq(-1,1,by=0.02)
x
[1] -1.00 -0.98 -0.96 -0.94 -0.92 -0.90 -0.88 -0.86 -0.84 -0.82 -0.80 -0.78
[13] -0.76 -0.74 -0.72 -0.70 -0.68 -0.66 -0.64 -0.62 -0.60 -0.58 -0.56 -0.54
[25] -0.52 -0.50 -0.48 -0.46 -0.44
not sure if you meant to use both 0.2 and 0.02, but i believe your
unexpected results are a floating point issue..
start here
http://r.789695.n4.nabble.com/That-dreaded-floating-point-trap-td3418142.html
and here
faq 7.31
---
Jeff NewmillerThe . . Go Live...
DCN:jdnew...@dcn.davis.ca.usBasics: ##.#. ##.#. Live Go...
Live: OO#.. Dead: OO#..
Dear useRs,
I've just encountered a serious problem involving the F-test being carried
out in aov() and anova(). In the provided example, aov() is not making the
correct F-test for an hypothesis involving the expected mean square (EMS) of
a factor divided by the EMS of another factor (i.e.,
The following code fails to train a nnet model in a random dataset using
caret:
nR - 700
nCol - 2000
myCtrl - trainControl(method=cv, number=3, preProcOptions=NULL,
classProbs = TRUE, summaryFunction = twoClassSummary)
trX - data.frame(replicate(nR, rnorm(nCol)))
trY -
Your subject line is patent nonsense. The aov() and anova() functions
have been around for decades. If they were doing something wrong
it would have been noticed long since.
You should realize that the fault is in your understanding, not in these
functions.
I cannot really follow your
I would guess that it's the brevity argument;
witness pnorm, pchisq, punif, ppois, pbinom, phyper, ...
But maybe Kurt knows better since he's down as the author
of pwilcox.
I can only guess it was the brevity argument also. The [...]wilcox functions
were added to S+ in June 1990, with no
Dear R user
I have data created using code below
b-matrix(2:21,nrow=4)
b[,1:3]=NA
b[4,2]=5
b[3,1]=6
Now the data is
b
[,1] [,2] [,3] [,4] [,5]
[1,] NA NA NA 14 18
[2,] NA NA NA 15 19
[3,] 6 NA NA 16 20
[4,] NA5 NA17 21
I want to
Hi,
b[b[,4]15 (b[,1]4|is.na(b[,1])) (b[,2]4|is.na(b[,2])),]
# [,1] [,2] [,3] [,4] [,5]
#[1,] 6 NA NA 16 20
#[2,] NA 5 NA 17 21
A.K.
- Original Message -
From: HJ YAN yhj...@googlemail.com
To: r-help@r-project.org
Cc:
Sent: Tuesday, March 5, 2013 9:33 PM
thanks for your answer.
The problem in my case was that Mac was loading R (32-bit) version rather
than R64 (64-bit).
Since many libraries in this Mac are compiled for 64 architecture. So
naturally the installation did not work at all.
Please make use of my case for future reference.
soichi
Dears, I am specifying a panel model with GLS function and corCAR1 and cor
AR1 parameters to correct for serial autocorrelation. But I have one
seemingly trivial question: those models fitted that way are random effects
models or pooling models?
Thanks again,
Tomas Notes
[[alternative
Dear Mark--
I've just spent an hour and a half reading chapters from Hadley's book. It
is phenomenal. Thanks for pointing it out to me
--andrewH
On Mon, Mar 4, 2013 at 9:04 PM, Mark Leeds marklee...@gmail.com wrote:
Hi Andrew: Not that I've gone through it all yet but the draft of hadley's
Thanks, David! That bookfinder.com search is awesome! I checked four sites
and the best price i found for the white Bokk used was $99 + $4 shipping.
This was a quarter of that. So i just bought it.
The Venables and Ripley book was actually part of my previous
budget-busting splurge. I agree with
Dear all,
I was reading last night the lm and the Formula manual page, and 'I have to
admit that I had tough time to understand their syntax. Is there a simpler
guide for the dummies like me to start with?
I would like to thank you in advance for your help
Regards
Alex
[[alternative
Thanks, Andrew! I'll put it on my list.
I have not been through much of it yet, but the exercises on count data are
excelent and at least one of them is immediately helpful to a current
project.
With appreciation, andrewH
On Mon, Mar 4, 2013 at 7:28 PM, Andrew Koeser arborkoe...@yahoo.com wrote:
Dear Patrick--
After the official Core Team's R manuals and the individual function help
pages, I have found The R Inferno to be the single most useful piece of
documentation when I have gotten stuck with a R problems. It is the only
introduction that seems to be aware of the ambiguities present
Victor,
It's a bit clunky, but you can use the 'pcls' function in package 'mgcv'
for this, by adapting the examples in the help file. The examples
themselves deal with inequality constraints imposing monotonicity, but
'pcls' also allows you to impose equality constraints. The examples are
86 matches
Mail list logo