On Tue, 17 Jun 2014 12:48:54 +
"Hiyoshi, Ayako" wrote:
> Dear Martyn and Professor Ripley,
>
> Thank you so much for your help. I used Window's large file search
> (it was useful! thank you), but there is no big files detected in C:
> drive . Perhaps I will have to reinstall Windows..
>
The
On Tue, 17 Jun 2014 08:58:12 PM Cleber Chaves wrote:
> Dear all,
> has a long time that I unsuccessfully try to obtain letters of each
equal
> means in a Tukey test of a GLM object, in a descending order. Could
> someone, please, help me?
>
> I try this:
>
> m1 <- glm(comp~pop)
> anova(m1,test="
Dear all,
has a long time that I unsuccessfully try to obtain letters of each equal
means in a Tukey test of a GLM object, in a descending order. Could
someone, please, help me?
I try this:
m1 <- glm(comp~pop)
anova(m1,test="F")
lsmeans(m1, cld~pop, adjust='tukey')
But the letters I obtain are n
On Tue, 17 Jun 2014 10:30:17 AM suganthiny cumarasamy wrote:
> Dear R mailing listers,
>
> I am try to find the different power calculation for
> p1=0.2 and p2=0.4 , with significant level=0.05 (one sided test)
> I would like to have a graph y -axis as a power and
> x-axis as a sample size .
>
If you are not determined to apply switch(), you can apply apply():
Data <- data.frame(V1=c(1, 5), V2=c(3, 4))
Data
## V1 V2
## 1 1 3
## 2 5 4
label <- data.frame(freq=c(1:5), lbl=c("One", "Two", "Three", "Four", "Five"))
label
## freq lbl
## 11 One
## 22 Two
## 33 Three
#
On 18/06/14 03:49, Christopher W Ryan wrote:
I upgraded to version 3.1.0 of R, still on Win XP, but that did not
solve the problem:
library(spatstat)
data(redwood)
str(redwood)
plot(redwood) # works OK
plot(Kest(redwood)) # works OK
plot(density(redwood)) # produces an "empty" plot
You were using 'rownames' incorrectly. Here is the correct way:
library(XLConnect)
wb <- loadWorkbook("writeWorksheet.xlsx", create = TRUE)
createSheet(wb, name = "CO2")
writeWorksheet(wb, CO2, sheet = "CO2",rownames="Row Names")
saveWorkbook(wb)
Jim Holtman
Data Munger Guru
What is the pro
That would depend on the actual type of object that your Data variable is. Use
the dput function so we can get your example days in a reproducible way [1].
[1]
http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example
-
Hi,
Is anyone travelling north after the useR! conference in Los Angeles and could
give me a ride?
I need to get to Medford, which is located along I5 in southern Oregon.
Thanks ahead,
Berry
PS: sorry to bother the help list, but I could not find anything specific for
useR.
On Jun 17, 2014, at 5:00 AM, Thomas Barningham wrote:
> Dear R users,
>
> I have a .txt file of time series data from the middle of 1999 through
> to the end of 2012. There are two columns in the file, the first is a
> decimal date, the second is an atmospheric tracer value and looks like
> this
Hello,
I have a simple Data table like below:
> Data
V1 V2
1 1 3
2 5 4
And I like to apply a switch to all the values of the table.
Exp for one value:
switch(as.character(Data[2,1]),
"1"="One",
"2"="Two",
"3"="Three",
"4"="Four",
"5"="Five"
)
What would b
Dear R mailing listers,
I am try to find the different power calculation for
p1=0.2 and p2=0.4 , with significant level=0.05 (one sided test)
I would like to have a graph y -axis as a power and
x-axis as a sample size .
I run this command for different value of power. and get the n and power
Dear R users,
I have a .txt file of time series data from the middle of 1999 through
to the end of 2012. There are two columns in the file, the first is a
decimal date, the second is an atmospheric tracer value and looks like
this:
Date APO
1999.47945560 -168.48
1999.47948410 -158.08
1999.4795
On Jun 16, 2014, at 11:51 PM, dila radi wrote:
> Dear R users,
>
> I used this codes to find the 1-day annual rainfall maximum
>
> bmFort <- blockmaxxer(Fort, blocks = Fort$year, which="Prec");bmFort
>
> and the output is given below:
>
> structure(list(obs = c(119, 507, 994, 1398, 1583, 1917
On Jun 17, 2014, at 9:20 AM, Knut Krueger wrote:
> Any Idea to prevent the warning message?
>
> library(XLConnect)
> wb <- loadWorkbook("writeWorksheet.xlsx", create = TRUE)
> createSheet(wb, name = "CO2")
> writeWorksheet(wb, CO2, sheet = "CO2",rownames=c("",1:84))
> saveWorkbook(wb)
>
> Warni
R users,
Abstracts are now being accepted for the 2015 ASA Conference on Statistical
Practice, February 19-21, New Orleans, Louisiana, USA. Each presentation
will have a 45-minute time slot.
If you are interested, you may submit your abstract on the website,
http://www.amstat.org/meetings/c
I'm running some spatial GAMs and am using the bam call (negbin family) and
am getting conflicting information on whether the model is converging or
not. When the model completes its run, I get a warning message that the
model did not converge. When I look at the object itself, I'm told that it
d
Suga,
Perhaps something like this ...
library(pwr)
n <- 50:100
power <- pwr.2p.test(h=ES.h(0.4,0.2), n=n, sig.level=0.05,
alternative="greater")$power
plot(n, power)
Jean
On Mon, Jun 16, 2014 at 10:44 AM, Suganthie Jeyaganth
wrote:
> Dear R mailing listers,
> I am try to find the different p
Any Idea to prevent the warning message?
library(XLConnect)
wb <- loadWorkbook("writeWorksheet.xlsx", create = TRUE)
createSheet(wb, name = "CO2")
writeWorksheet(wb, CO2, sheet = "CO2",rownames=c("",1:84))
saveWorkbook(wb)
Warning message:
In names(res)[1] <- colname :
number of items to repla
I am not familiar with the pamr.plotcv function, but in general if it
uses par(mfrow=c(2,1)) to set up for the multiple plots then your
problem with going back to the first plot is that you have lost the
other information (such as user coordinates) needed to add to the 1st
plot. You can see that w
I upgraded to version 3.1.0 of R, still on Win XP, but that did not
solve the problem:
library(spatstat)
data(redwood)
str(redwood)
plot(redwood) # works OK
plot(Kest(redwood)) # works OK
plot(density(redwood)) # produces an "empty" plot
# with just the fi
I made a stupid mistake. My note at the end is irrelevant as is the code
before the library statement.
Instead, the upper bound to your integral must be around .5, not 1.
> integrate(Gauss2F1,lower=0,upper=.5,H=0.7)
0.4954698 with absolute error < 5.5e-15
> integrate(Gauss2F1,lower=0,upper=.51
Hello R-helpers,
I am working through Simon Wood's GAM book and want to specify my own
knot locations (on even tens, i.e. 10, 20, 30, etc.). Then, I want to
compute a GAM on that area, and given the coefficients, reconstruct the
same P-spline that is drawn in plot(my_gam).
I'm failing.
Here
I ususally scan the digest for "surv", so missed your question on the first
round.
You caught predict with a case that I never thought of; I'll look into making
it smarter.
As Peter D said, the clogit function simply sets up a special data set and then calls
coxph, and is based on an identity
#First, correct your code
H <- 0.7;
t <- 1;
s <- seq(0,t,0.1);
x <- 1-t/s;
library(gsl)
Gauss2F1 <- function(H,x){
a <- H-1/2;
b <- 1/2-H;
c <- H+1/2;
# if(x>=0 & x<1){
#hyperg_2F1(a,b,c,x)
# }else{
hyperg_2F1(c-a,b,c,1-1/(1-x))/(1-x)^b
# }
}
#Here
#F <- Gauss2F1(a,b,c,x)
Don't forget the possibility of "hidden files", which will not be
shown using normal file-listing procedures. It may be that R has
stored those files as hidden files.
See, for example:
http://windows.microsoft.com/en-gb/windows/show-hidden-files
http://www.sevenforums.com/tutorials/394-hidden-fi
Dear Martyn and Professor Ripley,
Thank you so much for your help. I used Window's large file search (it was
useful! thank you), but there is no big files detected in C: drive .
Perhaps I will have to reinstall Windows..
Thank you so much for your replies.
Best wishes,
Ayako
___
On Mon, 16 Jun 2014 11:51:44 PM dila radi wrote:
> Dear R users,
>
> I used this codes to find the 1-day annual rainfall maximum
>
> bmFort <- blockmaxxer(Fort, blocks = Fort$year,
which="Prec");bmFort
>
> and the output is given below:
>
> structure(list(obs = c(119, 507, 994, 1398, 1583, 19
Dear all,
I am trying to make a numerical integral of a Hypergeometrical function 2F1
but it is not working. I dont know how to pass the arguments to my
function. Thank you beforehand.
Chuse.
Here is my code:
H <- 0.7;
t <- 1;
s <- seq(0,t,0.1);
x <- 1-t/s;
library(gsl)
Gauss2F1 <- function(H,
Hi Thomas,
Thank you very much for your help. I will try to calculate some summary
statistics and fit observed and simulated data using a time series (as your
examples). With observed data, I would use two values (e.g. 0 and 10),
because I do not have intermediate values.
Thanks!
Javier
2014-06-1
Dear Professor Ripley,
Thank you so much for your quick reply.
I tried 'dianame(tempdir())' and removed several 'Rtemp' and all other files.
Unfortunately it did not changed C: drive space much.
I am really sorry, but could there be other things stored in somewhere in C:
drive ?
I called IT
Hi,
Try
http://social.technet.microsoft.com/wiki/contents/articles/19295.windows-8-how-to-search-for-large-files.aspx
Martyn
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Hiyoshi, Ayako
Sent: 17 June 2014 11:40
To: Prof Brian R
On 17/06/2014 11:39, Hiyoshi, Ayako wrote:
Dear Professor Ripley,
Thank you so much for your quick reply.
I tried 'dianame(tempdir())' and removed several 'Rtemp' and all other files.
Unfortunately it did not changed C: drive space much.
I am really sorry, but could there be other things stor
tempdir() is specific to an R session.
Start up R and run
dirname(tempdir())
and look at that directory. Shut down R, then remove all old
files/folders in that directory, especially those beginning with 'Rtmp'.
An R process tries to clean up after itself but
- it cannot do so if it segfaul
There is a Mtrix Tester on the pages of Charlotte K. Hemelrijk, with
TauKr-Test Partial TAU Kr-TEst , MantelZ, R-Test andK TEst
Is there a similar package on R for that?
(
http://www.rug.nl/research/behavioural-ecology-and-self-organization/people/hemelrijk
scroll down to MatrixTester
Social
On 17.06.2014 10:37, Javier Rodríguez Pérez wrote:
Hi Thomas, Thank you very much for your help. I will try to calculate
some summary statistics and fit observed and simulated data using a
time series (as your examples). With observed data, I would use two
values (e.g. 0 and 10), because I do not
Hello!
I'm trying to generate grid-based landscape mimicking observed data. For
this purpose I'm using simecol and adapting the "CA" (stochastic celular
automaton) model included as example in vignette. I think this example
could be nice given observed data has spatial autocorrelation. My observed
Dear R users,
Hello, I am new to R and confused about my PC's memory space after using R a
while.
My PC is Windows 8, RAM 8GB. I have run some analyses using commands like
"vglm", "aalen(Sruv())", lm()... some regressions.
I also use Stata, and when I tried to run Stata (big file), Stata
Dear R users,
I used this codes to find the 1-day annual rainfall maximum
bmFort <- blockmaxxer(Fort, blocks = Fort$year, which="Prec");bmFort
and the output is given below:
structure(list(obs = c(119, 507, 994, 1398, 1583, 1917), tobs = c(119,
142, 264, 303, 123, 91), month = c(4, 5, 9, 10, 5
Hi,
I have no example at hand, but the usual way could be indeed to
calculate a time series of an adequate statistic (e.g. spatial
statistic) from both, the observed and the simulated data and then to
apply standard model fitting.
Thomas
__
R-help@r-p
40 matches
Mail list logo
