On 06/02/2013 22:47, Roy Mendelssohn - NOAA Federal wrote:
Hi All:
We are developing (actually updating) a piece of code that runs in
R.
Due to some limitations of some the libraries we use, at the present
time we would be limited to the 32-bit version of R on Windows. Does
anyone have a feel (
In R I have 3 variables x,y and z that contain:
| x y z
[1,] 1 12 122
[2,] 1 13 113
[3,] 1 11 111
[4,] 3 13 123
[5,] 4 14 134
[6,] 5 15 155|
I am creating a coplot with the following syntax:
|given.depth<- co.intervals(z, number=3, overlap=.1)
coplot(x~y|z,pch=19,col
x<-matrix(c(1,11,11,2,3,4),2,3)
y<-matrix(x(7,8,9,10,11,12),2,3)
if you want to merge them by row, then use:
rbind(x,y)
if you want to merge them by column, then use
cbind(x,y)
hope this can help you out*
Alice Xiong
On Wed, Feb 6, 2013 at 8:55 AM, Mat wrote:
> Hello together,
>
> i have
Dear sir,
Which mailing list is appropriate to ask for kstest in spatstat?
Regards,
Hiroshi Saito
*
Hiroshi Saito:
E-MAIL: saito.hiro...@lab.ntt.co.jp
http://www9.plala.or.jp/hslab/
PHONE: +81 422 59 4300
FAX: +81 422 59 6364
__
to start with, you should now what kinds of independent variables in your
data=Rail.
it looks like you depend variable is "travel".
is "Rail" the only IV in your model? it seems has been classifed to L, Q, C.
from the output, it seems the final model should be:
E(travel)=66.5+54*Rail.L-4.7 Rail.Q-
Hi,
Anyway, just using some random combinations:
dnew<- expand.grid(4:10,5:10,6:10,3:7,4:5,6:8)
names(dnew)<-c("m","n","x1","y1","x","y")
resF<- cbind(dnew,d2[rep(1:nrow(d2),nrow(dnew)),])
row.names(resF)<- 1:nrow(resF)
head(resF)
# m n x1 y1 x y m1 n1 cterm1_P1L cterm1_P0H
#1 4 5 6 3 4 6
After applying the NLS for a model like y=exp(a*x), and I get a result
showing the summary as:
Estimate Std. Error t value Pr(>|t|)
2.6720 1.4758 1.811 0.3212
My question is what this t-statistics tests? And what's the meaning of Pr?
New to R. Thanks.
[[alternative HTML version d
Blast! -- dumb typo:
On Wed, Feb 6, 2013 at 8:14 PM, Bert Gunter wrote:
> Just wanted to point out...
>
> On Wed, Feb 6, 2013 at 7:20 PM, Rich Shepard wrote:
>> On Wed, 6 Feb 2013, David Winsemius wrote:
>>
>>> Wouldn't this just be:
>>> save(itis.table, file="/usr/lib/R/library/bio.infer/data/i
See abort() in the R.oo package, with some discussion here
https://stat.ethz.ch/pipermail/r-devel/2012-September/064838.html
/Henrik
On Wed, Feb 6, 2013 at 7:05 PM, Jeff Newmiller wrote:
> stop() does not have momentum... it is designed to play nice.
>
> Perhaps you want ?options(error=).
> Alte
On Wed, 6 Feb 2013, David Winsemius wrote:
Wouldn't this just be:
save(itis.table, file="/usr/lib/R/library/bio.infer/data/itis.ttable.rda")
David,
Almost certainly. I did not find the save function in the indices of my R
library or the R data doc. Now I know it exists I'll read about it th
On Wed, 6 Feb 2013, MacQueen, Don wrote:
My sense is that even if this were possible (which I doubt),
Don,
It is possible.
... it would be outside the scope of what would be considered appropriate
in R (assuming I understand the question). After all, it's someone else's
package, theirs to
Dear All,
Here is a hypothetical sample (sorry for the clumsy code):
A1 <- matrix(1:5, nrow=5, ncol=1)
A2 <- matrix(6:10, nrow=5, ncol=1)
A3 <- matrix(11:15, nrow=5, ncol=1)
A4 <- matrix(16:20, nrow=5, ncol=1)
A5 <- matrix(21:25, nrow=5, ncol=1)
A6 <- matrix(26:30, nrow=5, ncol=1)
B1 <- matrix(c
stop() does not have momentum... it is designed to play nice.
Perhaps you want ?options(error=).
Alternately, you want to trigger a core dump?
---
Jeff NewmillerThe . . Go Live...
DCN
is it possible to throw a stop() that is so hard that it will escape
even tryCatch?
/iaw
Ivo Welch (ivo.we...@gmail.com)
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R
Hello,
I drew a plot of weight and height of people and fitted it with a
polynominal regression x^2.
(using curve())
Now I would like to draw the confidence interval line for the fitted curve.
Please kindly advise the code for the purpose.
Thank you.
Elaine
[[alternative HTML version de
On Feb 6, 2013, at 4:53 PM, Rolf Turner wrote:
>
> Perhaps you want to look at
>
>?plotmath
>
>cheers,
>
>Rolf Turner
Definitely ... he _will_ want to look at ?plotmath, but there is a lot on that
page and teh parts that explain the parsing of expressions were not at all
c
Hello,
If you know what a lowpass filter is, you should be able to use "filter"
from "stats" package.
And you can have a look at:
?triang
HTH,
Pascal
Le 07/02/2013 09:50, Janesh Devkota a écrit :
Hello R users,
I am trying to use R to do the low pass filter analysis for the tidal data.
I
Perhaps you want to look at
?plotmath
cheers,
Rolf Turner
On 02/07/2013 11:41 AM, David Arnold wrote:
Hi,
I'd like to add "alpha" in latex code to my image. Any suggestions?
library(UsingR)
BagA = c(rep(10,6),rep(20,5), rep(30,4),rep(40,3),rep(50,2),60,70)
BagA
BagB = c(10,
Hello R users,
I am trying to use R to do the low pass filter analysis for the tidal data.
I am a novice in R and so far been doing only simple stuffs on R. I found a
package called signal but couldn't find the proper tutorial for the low
pass filter.
Could anyone point me to the proper tutorial
For what it is worth Spencer, I can start rJava in both 32-bit R and
64-bit R for Windows 7. [And could even before Simon fixed the error
message).
And yes, I have both 32-bit Java and 64-bit Java 1.7.0_13 installed.
They should be separate entries under your control panel.
Rob
--
Robert
Since no-one else seems to have responded to your question, I'll
chip in with my two bob's worth:
(1) You can probably get what you might want to get from predict()
by looking at X$fitted.values --- where X is the object returned by
segmented().
(2) But predict() is probably the wrong way to
Most Windows computers in my neighborhood have been upgraded to 64 bit.
Clint BowmanINTERNET: cl...@ecy.wa.gov
Air Quality Modeler INTERNET: cl...@math.utah.edu
Department of Ecology VOICE: (360) 407-6815
PO Box 47600
On Feb 6, 2013, at 9:23 PM, David Winsemius wrote:
>
> On Feb 6, 2013, at 9:01 AM, Nicole Ford wrote:
>
>> hello, all.
>>
>> I am having his problem. i looked in archives and could not find a solution.
>
> Did you look in the rhelp archives or in the package archives?
>
>>
>> --- Please
Hi All:
We are developing (actually updating) a piece of code that runs in R. Due to
some limitations of some the libraries we use, at the present time we would be
limited to the 32-bit version of R on Windows. Does anyone have a feel (or
even real knowledge) if most people these days run R
Hi,
I'd like to add "alpha" in latex code to my image. Any suggestions?
library(UsingR)
BagA = c(rep(10,6),rep(20,5), rep(30,4),rep(40,3),rep(50,2),60,70)
BagA
BagB = c(10, 20,rep(30,2),rep(40,3),rep(50,4),rep(60,5),rep(70,6))
BagB
par(mfrow = c(2, 1))
DOTplot(BagA)
abline(v=60,lwd=2,lty=2,col="
On Feb 6, 2013, at 9:01 AM, Nicole Ford wrote:
> hello, all.
>
> I am having his problem. i looked in archives and could not find a solution.
Did you look in the rhelp archives or in the package archives?
>
> --- Please select a CRAN mirror for use in this session ---
> Warning message:
> pa
You can use system.file to fine the package root, then
data file path relative to this.
For example
system.file(package = "stats") # The root of package 'stats'
See ?system.file
On 6 February 2013 22:02, David Winsemius wrote:
>
> On Feb 4, 2013, at 11:37 AM, Rich Shepard wrote:
>
>> The bio.in
On Feb 4, 2013, at 11:37 AM, Rich Shepard wrote:
> The bio.infer package contains a data frame
> /usr/lib/R/library/bio.infer/data/itis.ttable.rda that needs to be modified.
> After loading the bio.infer package and attaching the data frame with the
> data() function, I wrote the data frame to a
On Feb 6, 2013, at 3:16 AM, Laura de Mergelina wrote:
> Where can I mark the option to not receive emails from users?
I didn't mean to hit reply, but rather to copy the link that is at the bottom
of every posting. ( So I guess my "empty reply actually could be construed as
being relevant.)
>
On Feb 6, 2013, at 3:16 AM, Laura de Mergelina wrote:
> Where can I mark the option to not receive emails from users?
>
> --
> Laura de Mergelina
>
> __
> R-help@r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do re
My sense is that even if this were possible (which I doubt), it would be
outside the scope of what would be considered appropriate in R (assuming I
understand the question). After all, it's someone else's package, theirs
to control and maintain, etc. What if it's an installation of R on a
multi-use
Hi everyone,
I'm trying to vectorize an application of lda to each 2D slice of a 3D
array, but am running into trouble: It seems there are quite a few 2D
slices that trigger either the "variables are collinear" warning, or worse,
trigger a "variable appears to be constant within groups" error and
Hi Mat,
just have a look at ?rbind
cheers.
Am 06.02.2013 15:55, schrieb Mat:
> Hello together,
>
> i have probably a easy question, but how can i sum 2 data.frames among each
> other.
> I have 2 data.frames which look like this one:
>
> Cu.No. place level
> 1 123
On Feb 6, 2013, at 11:24 AM, David Winsemius wrote:
>
> On Feb 6, 2013, at 8:46 AM, Christofer Bogaso wrote:
>
>> Hello again,
>>
>> I was looking for some way on How to delete repeated appearance in a
>> String. Let say I have following string:
>>
>> Text <- "ahsgdvasgAbcabcsdahj"
>>
>> Her
Hi Christopher,
what is the rule to omit "ah" which is also repeated in Text?
The following might be a start:
Text <- "ahsgdvasgAbcabcsdahj"
#finds first repetion of substrings of length 2 or more, here "ah"
gsub("(?i)([a-z]{2,})(.*)\\1","\\1\\2",Text,perl=T)
#finds all repetions of substrings of
On Feb 6, 2013, at 8:46 AM, Christofer Bogaso wrote:
> Hello again,
>
> I was looking for some way on How to delete repeated appearance in a
> String. Let say I have following string:
>
> Text <- "ahsgdvasgAbcabcsdahj"
>
> Here you see "Abc" appears twice. But I want to keep only 1
> occurrenc
Have you read An Introduction to R or other R tutorial? DO you know
about ordered factors and contrasts in linear models?
-- Bert
On Wed, Feb 6, 2013 at 5:00 AM, jing tang wrote:
> Hi,
> I am reading the book "Mixed Effects Models in S and S-Plus" and come
> across an example with the Rail data.
Hi,
maxN<- 9
res3new<- aggregate(.~m1+n1,data=res2[,c(1:2,12:13)],max)
d2<-res3new[res3new[,3]<0.01 & res3new[,4]<0.01,]
d2
# m1 n1 cterm1_P1L cterm1_P0H
#2 3 2 0.00032 0.0025
(d2$m1+2)
#[1] 5
maxN-(d2$n1+2)
#[1] 5
dnew<-expand.grid((m1+2):(maxN-(n1+2)),(n1+2):(maxN-m),0:m1,0:
hello, all.
I am having his problem. i looked in archives and could not find a solution.
--- Please select a CRAN mirror for use in this session ---
Warning message:
package ‘manipulate’ is not available (for R version 2.15.2)
thank you in advance.
_
Dear List,
It's the first time that I try to request your support.
I intent do use the BRB (Biased Random Bridge) function with a boundary.
I believe I have made sure to fulfill all the rules for boundaries (segment
length, angle of junction)
When calculating the UD I get the following Error:
Hi John,
Do you want to minimize correlation or |correlation|?
Here is a simple code, with example, that shows how to optimize the magnitude
(that is corr^2) of correlation between y and x^p, where p is a positive real
number. You can modify it to suit your needs.
corxy.opt <- function(x, y,
Hi,
I got an error message with:
vlist <- apply(mm, list)
Error in match.fun(FUN) : argument "FUN" is missing, with no default
#assuming that
vlist <- apply(mm,2,list)
mapply("%*%",mlist,vlist[1:2],SIMPLIFY=FALSE)
#[[1]]
# [,1]
#[1,] 19
#[2,] 22
#[3,] 25
#[4,] 28
#
#[[2]]
# [,1]
Hello together,
i have probably a easy question, but how can i sum 2 data.frames among each
other.
I have 2 data.frames which look like this one:
Cu.No. place level
1 123London A
2111 Paris B
Cu.No. place level
1 333
Where can I mark the option to not receive emails from users?
--
Laura de Mergelina
__
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 provi
To detect heteroscedasticity for a multiple linear OLS regression (no time
dependencies):
What if the residuals vs. fitted values plot shows well behaved residuals
(cloud) - but the some of the x versus residuals plots are a megaphone?
Also, it seems that textbooks and internet tutorials in R
Hi Martin,
thanks for the answer.
I do agree ; in my example the initalize methods are not very releavant. It
was actually a toy example. I came across the issue in the following case
- I developped an object (say DB - for database) that was able to
handle in a nice way a connection and que
Hi,
I am reading the book "Mixed Effects Models in S and S-Plus" and come
across an example with the Rail data.
I tried to use lm(travel~Rail,data=Rail) and got the following result:
Call:
lm(formula = travel ~ Rail, data = Rail)
Residuals:
Min 1Q Median 3Q Max
-6.6667 -1.
>> Duncan Murdoch
>> on Mon, 04 Feb 2013 13:23:57 -0500 wrote:
DM> On 04/02/2013 12:59 PM, Ross Boylan wrote:
>> I am tweaking an R package for which I have the source; the relevant code
>> is in R not C. I'm making changes to the package code.
>>
>> What is the best workflow for
Grrr ... new trial with code here: http://pastebin.com/RjHNNG9J
Maybe the amount of inline-code prevented posting?
Hello,
I am writing a simple peak detector and it works quite well ... however
there's one special case below, that I can't get my head wrapped around ...
the problem is in the "De
1. Please do not double post.
2. This is a statistics, not an R question. Post on a statistics list,
like stats.stackexchange.com instead.
3. Is this a homework question? -- we don't do homework here.
-- Bert
On Wed, Feb 6, 2013 at 8:59 AM, Aminreza Aamini wrote:
> Dear all,
>
> How can i obta
Dear all,
How can i obtain odds ratio in logistic regression when the model
contains interaction terms in R?
how can i obtain OR for a special case in R?
Thanks in advance for any help.
Amin
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mail
Hi Eik,
thank you so much - it works perfectly!
Thank you and best wishes
Alain
Â
Eik Vettorazzi hat am 6. Februar 2013 um 17:01
geschrieben:
> Hi Alain,
> here is a one-liner for a df without the rowSum column
>
> df<-data.frame(id=c("x01","x02","x03","x04","x05","x06"),a=c(1,2,NA,4,5,6
Dear all,
How can i obtain odds ratio in logistic regression when the model
contains interaction terms in R?
how can i obtain OR for a special case?
Thanks in advance for any help.
Amin
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/l
Hello again,
I was looking for some way on How to delete repeated appearance in a
String. Let say I have following string:
Text <- "ahsgdvasgAbcabcsdahj"
Here you see "Abc" appears twice. But I want to keep only 1
occurrence. Therefore I need that:
Text_result <- "ahsgdvasgAbcsdahj" (i.e. the f
Thanks Rolf and Arun! -David
On Wed, Feb 6, 2013 at 6:13 AM, arun wrote:
> Hi,
>
> I got an error message with:
> vlist <- apply(mm, list)
>
> Error in match.fun(FUN) : argument "FUN" is missing, with no default
> #assuming that
> vlist <- apply(mm,2,list)
>
> mapply("%*%",mlist,vlist[1:2],SIM
Hi Alain,
here is a one-liner for a df without the rowSum column
df<-data.frame(id=c("x01","x02","x03","x04","x05","x06"),a=c(1,2,NA,4,5,6),b=c(2,4,6,8,10,NA),c=c(NA,3,9,12,NA,NA))
(df$wSum<-apply(sweep(df[,-1],1,rowSums(df[,-1],na.rm=T),"/"),1,function(x)sum(x*w,na.rm=T)))
hth.
Am 06.02.2013 1
a nice aspect about roxygen (and perl pod) is that it compiles into
the standard .Rd files. it's not a substitute, but a complement tool.
it may also make a good start for docs when one starts writing code,
and be eventually yanked out in favor of direct changes in the .Rd
files that it has creat
On 06/02/2013 9:49 AM, Hadley Wickham wrote:
>> I have heard of people using noweb to do this, but I can't point to any
>> examples. I'd actually recommend against it. Good documentation files
>> don't
>> make good source files.
>
> the compiler package in base R is, apparently, developed using
On 06/02/2013 9:44 AM, Martin Morgan wrote:
On 02/06/2013 03:31 AM, Duncan Murdoch wrote:
> On 13-02-05 7:43 PM, ivo welch wrote:
>> Dear R experts---
>>
>> after many years, I am planning to give in and write my first R
>> package. I want to combine my collection of collected useful utility
>>
>> I have heard of people using noweb to do this, but I can't point to any
>> examples. I'd actually recommend against it. Good documentation files
>> don't
>> make good source files.
>
> the compiler package in base R is, apparently, developed using noweb
> https://svn.r-project.org/R/trunk/src/
On 02/06/2013 03:31 AM, Duncan Murdoch wrote:
On 13-02-05 7:43 PM, ivo welch wrote:
Dear R experts---
after many years, I am planning to give in and write my first R
package. I want to combine my collection of collected useful utility
routines.
as my guide, I am planning to use Friedrich Leis
> Other people have recommended Roxygen, but honestly I haven't seen a package
> documented with Roxygen where the documentation was adequate.
> It looks as though it's great to get initial documentation created, but does
> not appear to encourage followup.
I don't think that's a problem with roxy
Dear R-List,
I am sure there must be a very simple way to do this - I just do not know how...
This is what I want to do:
#my dataframe
df<-data.frame(id=c("x01","x02","x03","x04","x05","x06"),a=c(1,2,NA,4,5,6),b=c(2,4,6,8,10,NA),c=c(NA,3,9,12,NA,NA),sum=c(3,9,15,24,15,6))
id a b
Thank you, David. I believe you have made the thing more clear to me.
Best Regards,
Xie YM
在 2013-02-06 16:14:18，"David Winsemius" 写道：
>
>On Feb 5, 2013, at 10:48 PM, 谢一鸣 wrote:
>
>> Thanks David, your solution perfectly solve my problem. But Do you
>> know what exactly the mechanism of "Compu
To answer my own question: I thought for a while that centering is the
solution to the problem, but it is not:
> integrate(dnorm, 0, Inf, mean = 25)
3.187474e-05 with absolute error < 5.9e-05
> integrate(dnorm, -25, Inf, mean = 0)
3.187474e-05 with absolute error < 5.9e-05
Here "Don't do this"
On 13-02-05 7:43 PM, ivo welch wrote:
Dear R experts---
after many years, I am planning to give in and write my first R
package. I want to combine my collection of collected useful utility
routines.
as my guide, I am planning to use Friedrich Leisch's "Creating R
Packages: A Tutorial" from Sep
> I've checked my iteration loop with other dataset, i've not
> problems with integer numbers or other dataset with decimals
> but not all decimals dataset.
Any R data set that is 'not all decimals' will be represented internally as
floating point* (eg c(1, 2, 3, 1.5) consists of four floating
I've checked my iteration loop with other dataset, i've not problems with
integer numbers or other dataset with decimals but not all decimals dataset.
I would like it works with all type of data.
Thanks
--
View this message in context:
http://r.789695.n4.nabble.com/How-to-menage-decimals-in-or
Hi there,
>
> I am using the function pgls in the CAPER package, everything seems to run
> fine except for one of my variables.
> When using this variable as continuous the anova works without problem,
> but if using the same variable as factor (in fact, this is what it is), the
> program returns
I have a problem with the LaplacesDemon.hpc which terminates with an
error. I am not sure if this is in the parallel or the LaplacesDemon
package. And if it is a unsuffiency of my code or of the packages.
The code runs WELL as NON paralellized code.
> Fit <- LaplacesDemon.hpc(Model, Data=M
Hi
This is an error I get with your loop code due to this line
capture.output(print(lll), file = "C:/Chi_Square_fix/temp.txt", append = TRUE)
Error in file(file, if (append) "a" else "w") :
cannot open the connection
In addition: Warning messages:
1: In chisq.test(lll[["mat4"]], correct = FALSE
On 02/06/2013 07:51 PM, mary wrote:
I'm sorry,
I have a data set contains numbers with more decimal places, the data are
available in R "data (milk)".
I made a loop to get a robust correlation matrix but in calculations,
decimals cause me problems ( especially I do not know how R rounds numbe
Hi Chunlin,
You have to write your own myCmatrix.
Cheers
Andrew
On Wednesday, February 6, 2013, chunlin liu wrote:
> Thanks, Andrew.
>
>
>
> I studied the function "estimable" in package "gmodels"
>
> and realized that it could combine the model coefficients together
>
> if a correct cm matri
I'm sorry,
I have a data set contains numbers with more decimal places, the data are
available in R "data (milk)".
I made a loop to get a robust correlation matrix but in calculations,
decimals cause me problems ( especially I do not know how R rounds numbers
by default) so the diagonal should
Dear Rafael,
Thank You very much for Your reply.
I need to use non-overlapping blocks and such that each block ends with an
element having an index multiple of 3 in the original series.
For example if I have the numbers: 1 2 3 4 5 6 7 8 9, the samples I could take
are:
123 456 789, 456 789
Hello,
RsiteSearch('garch')
require(sos)
findFn('garch')
HTH,
Pascal
Le 06/02/2013 16:50, jpm miao a écrit :
Hi everyone,
Do we have GARCH-in-mean VAR in R?
Is there any GARCH package in R?
Thanks,
Miao
[[alternative HTML version deleted]]
__
Thanks, Andrew.
I studied the function "estimable" in package "gmodels"
and realized that it could combine the model coefficients together
if a correct cm matrix is provided. For example, if
mylm <- lm(y ~ x1 + factor(x2) + x3*factor(x4), mydata)
Assume that factors x2 and x4 have 6 and
When I wrote the message earlier, I didn't actually try running example for
function 'stop' from R CMD BATCH.
Now, I have run a script with this content from R CMD BATCH.
example(stop)
It turns out that execution of the example still didn't continue after first
error. But, there was no "Execut
On Feb 5, 2013, at 10:48 PM, 谢一鸣 wrote:
Thanks David, your solution perfectly solve my problem. But Do you
know what exactly the mechanism of "Computing on the language" in
pdf file: "R language definition", Can that "THING" do the same job?
I wouldn't consider that to be "computing on t
79 matches
Mail list logo
