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
bas
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
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 wrote:
> The book that help
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 HTM
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 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 wrote:
> Hi Andrew: Not that I've gone through it all yet but the draft of hadley's
> book at https://g
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 H
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
2013
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
To: r-help@r-project.org
Cc:
Sent: Tuesday, March 5, 2013 9:33 PM
Subject: [R] How to
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 t
> 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
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 convol
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 <- runif(1)*trX[,1
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., inste
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
http://cran.r-project.org/doc/FAQ/R-FAQ.html#Why-doesn_0027t-R-think-these-numbers-ar
faq 7.31
---
Jeff NewmillerThe . . Go Live...
DCN:Basics: ##.#. ##.#. Live Go...
Live: OO#.. Dead: OO#.. Playing
Research Engine
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 -
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, decre
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 +
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
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.
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
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)
if(length(
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 wrote:
> Hi there,
>
> I am trying to write a to
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 wrote:
> Hi there,
>
> I am trying to write a tool which involves impl
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
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<- unlist(lapply(seq_len(nrow(dat1
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<- read.csv("rightest.csv",sep=",
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 wit
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
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:
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(ben<10, NA, -ben))
user system elapsed
1.310.241.56
system.time({BEN1 <- -ben; BEN1[BEN1> -10] <- NA})
u
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 (i.e
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 simulati
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 A&B and B&C.
--
View this message in
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
d
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 a
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 wind
On Tue, Mar 5, 2013 at 7:07 PM, Marbles 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 mean valu
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 g
On Tue, Mar 5, 2013 at 3:30 PM, Matthijs Daelman
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
> 2010-01-08, 270
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 <- matri
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 l
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
Thi
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 wrote:
> Hi,
>
> I'm having problems with the Zelig package - when using the below R
> displays the follwing message (I'm running R i386
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.
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]]
__
R-help@r-
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 proce
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
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 mo
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 weir
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
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.
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 [mailto
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, an
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
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 yo
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
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 <- 1
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
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 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 wrote:
>
> HI,
>
> Assume I have the d
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 ?
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)
c(c(tail(
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 such
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. 201
Dear Julien,
Check
citation('stats')
HTH,
Jorge.-
On Wed, Mar 6, 2013 at 12:05 AM, Julien Mvdb 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 Vettori
>
> [
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
pl
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 Wi
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
https:
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, Аскар Нысанов wrote:
>
>
> Hi everyo
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 o
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.
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 p
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
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
han
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 Bi
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 n
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
thierr
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
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
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 rea
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 bu
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 wo
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 d
On Mon, Mar 4, 2013 at 1:39 PM, Tammy Ma 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 61%
>
>
86 matches
Mail list logo
