Dear Mr/ms
My name is dwi nabila. I am a postgraduate student of Statistics Department
from University of Padjadjaran, Indonesia. Now I am at fourth semester, and I
will write a thesis to complete my study. My adviser give me advice to write
about FPCR. I have gcm data(x) , where there are 64 v
Hi, I could not open the file.
Can you load the data into R and then use the dput function?
The first link deals with dput. The latter 2 links are for dealing with
NAs in data.
http://www.inside-r.org/r-doc/base/dput
dput {base}
http://r.789695.n4.nabble.com/subset-without-removing-NAs-td348201
Or use an offset
lm( y ~ x+offset(x), data = dat)
The offset gives x a coefficient of 1, so the coefficient of x in this
model is the difference between the coefficient of x in the model without
an offset and 1 -- the thing you want.
-thomas
On Wed, May 1, 2013 at 2:54 PM, Paul Johnson wro
It is easy to construct your own test. I test against null of 0 first so I
can be sure I match the right result from summary.lm.
## get the standard error
seofb <- sqrt(diag(vcov(lm1)))
## calculate t. Replace 0 by your null
myt <- (coef(lm1) - 0)/seofb
mypval <- 2*pt(abs(myt), lower.tail = FALSE,
Greetings to r-help land.
I've run into some program crashes and I've traced them back to methods()
behavior
after the package gdata is loaded. I provide now a minimal re-producible
example. This seems bugish to me. How about you?
dat <- data.frame(x = rnorm(100), y = rnorm(100))
lm1 <- lm(y ~ x
Hello,
I am work with a linear regression model:
y=ax+b with the function of lm.
y= observed migration distance of butterflies
x= predicted migration distance of butterflies
Usually the result will show
if the linear term a is significantly different from zero based on the
p-value.
Now I
Hi All,
i've problem to understand how to work with R to generate a hierarchical
clustering
my data are in a csv and looks like :
idcode,count,temp,sal,depth_m,subs
16001,136,4.308,32.828,63.46,47
16001,109,4.31,32.829,63.09,49
16001,107,4.302,32.822,62.54,47
16001,87,4.318,32.834,62.54,48
16002
On Tue, Apr 30, 2013 at 5:12 PM, Dorte Klerke wrote:
>
> Hi I have been trying for days now to read the contents of this xml file into
> R. I tried the simpler solution xmlToDataFrame, which worked on another
> file, but somehow it doesn't work for this file.
>
> How is this possible?
How is t
Looks like my post might get some problems, so I re-wrote my question plus some
new one...
I noticed that the points on the biplot are not exactly the same as the
predicted values.
Another relevant question: should I expect that all the vector points have the
same length if chose parameters
Hi,
You could also do:
v<- read.table(text="
Name Year_1_value Year_2_value Year_3_value
A 1 2 3
B 2 7 19
C 3 4 2
D 10 7 6
E 4 4 5
F NA 3 6
",sep="",header=TRUE,stringsAsFactors=FALSE)
names(v)[-1]<-gsub("(.*\\d+)_.*$","\\1",names(v)[-1])
v2<- v
v1<-reshape(v,direction="long",varying=2:4,sep="_")
1. Read "an Introduction to R" or other R tutorial to learn how R works.
2. You apparently wish to apply a function, f, to each row of a data
frame or matrix to classify it as growing, declining, etc. Only you
know what that function should look like. Write it.
3. Apply it using ?apply or perhaps
Here is one way to, for each row in the data.frame v, regress the numbers in
columns 2 through 4 on the numbers 1 through 3, storing only the slopes, and
then creating a column saying if the slope is greater than zero or not.
> v[,"Beta"] <- vapply(seq_len(nrow(v)),
Folks,
This is probably a "help me google this properly, please"-type
of question.
In TIBCO Spotfire, there is a procedure called "line
similarity". I use this to determine which observations show a growing, stable
or declining pattern... sort of like a mini-reg
Tena koe Jie
Try
x <- ' 12"ab"34 '
plot(1:10, xlab=x)
HTH ...
Peter Alspach
-Original Message-
From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org] On
Behalf Of Jie
Sent: Wednesday, 1 May 2013 6:46 a.m.
To: r-help@r-project.org
Subject: [R] Quote as element of a ve
Hi,
Try this:
dat1<-read.table(text="
ID X0 X1 X2 X3 X4 X5 X6 X7 X8 X9 X10 X11 X12 X13 X14 X15
1 5184 0 0 1 0 0 0 0 0 0 1 0 0 0
0 0 0
2 6884 0 0 1 0 0 1 0 0 0 0 0 0 0
0 0
y = c("1", "2", '"', "a", "b", '"', "3", "4")
y
#[1] "1" "2" "\"" "a" "b" "\"" "3" "4"
A.K.
- Original Message -
From: Jie
To: "r-help@r-project.org"
Cc:
Sent: Tuesday, April 30, 2013 2:46 PM
Subject: [R] Quote as element of a vector/list
Dear All,
I would like to store
Dear All,
I would like to store quote as part of an vector.
For instance, I would like to get an character object as
x = " 12"ab"34 "
or
y = c("1", "2", """, "a", "b", """, "3", "4")
Is that possible? Thank you.
Best wishes,
Jie
[[alternative HTML version deleted]]
Another alternative might be to post the the r-sig-mixed-models list
gls() is part of nlme.
-- Bert
On Tue, Apr 30, 2013 at 11:12 AM, Bert Gunter wrote:
> Actually, Jeff, it may not be.
>
> The OP has likely overfit -- as a guess, Depth is essentially
> continuous and as.factor() produces too m
Actually, Jeff, it may not be.
The OP has likely overfit -- as a guess, Depth is essentially
continuous and as.factor() produces too many levels. I would recommend
that the OP find a local statistical consultant to help him, as he
seems out of his depth, statistically. Failing that, he should
pro
Hi I have been trying for days now to read the contents of this xml file into
R. I tried the simpler solution xmlToDataFrame, which worked on another
file, but somehow it doesn't work for this file.
I have tried different other solutions, but I just can't seem to get it
quite right.
doc <-
Dear Rxperts,
Is there a simpler way to generate multipanel grouped individual profile
plots? All individuals of a group within a panel have the same color. As
of now I am using lattice::xyplot to get the desired effect. Please feel
free to suggest other ideas.
Also, I am trying to create a gener
I noticed that the points on the biplot are not exactly the same as the
predicted values.
Could any body give me a hint about why?
Thanks.
To: Jim Lemon
Cc: "r-help@r-project.org"
Sent: Tuesday, April 30, 2013 9:51 AM
Subject: Re: [R] biplot for principa
Hi Olvia,
Welcome to the R-help list.
Are you seriously proposing to add error bars to a stacked barplot? Why?
Do you mean a dodged barplot?
My first thought is that there probably are several better ways of presenting
your data than stacked bar charts--I sometimes think that anything other t
That is unlikely in the absence of a reproducible example [1].
[1]
http://stackoverflow.com/questions/5963269/how-to-make-a-great-r-reproducible-example
---
Jeff NewmillerThe . . Go
Hello everyone,
I'm trying to perform some trend analysis using the package "Kendall" and for
some observations I get the following error "WARNING Error exit, tauk2, IFAULT
= 12"
Could anyone tell me how to create a logical vector of observations having its
element equal to TRUE if this error h
On Apr 30, 2013, at 6:00 AM, arun wrote:
> Hi,
> May be this helps:
> funcName<- function(df1, x){
> whatCol=df1[[x]]
> print("Got it")
> print(whatCol)
> }
>
> funcName(df,"ColA")
> #[1] "Got it"
> #[1] 1 2 3 4 5
> funcName(df,"ColB")
> #[1] "Got it"
> #[1] A B C D E
> #Levels: A B C D E
Hi,
I was wondering if there is a package/function for Panel Granger
non-causality tests? I am interested in Toda-Yamamoto procedure in panel
data setting.
Thank you,
--
View this message in context:
http://r.789695.n4.nabble.com/Panel-Granger-Causality-Tests-tp4665834.html
Sent from the R h
I am trying to find the percentage of the parameters explaining the bacterial
community composition. I have one data matrix with relative abundance of OTUs
and one with environmental parameters. I used varpart in vegan package but the
values in the venn diagram is bigger than 100% in total.How i
Hi there,
I've been battling with an extension of this in my own data: getting
appropriate error bars once data is stacked in a bar graph.
(original question:
http://r.789695.n4.nabble.com/ggplot2-se-variable-in-geom-errorbar-s-limits-
td3311176.html). It wouldn't let me reply to that thread.
A
Dear R users,
For the purpose of validating a prediction model using validate() from the rms
package, I am running into some trouble with using the fastbw() function
breaking up natural groups of variables.
Is there any way I can specify to keep certain variable together? In
particular, if int
Question about package Coxme:
I develop a cox model that includes a variable "treatment" with 3 levels (A, B,
C):
> model_alea_int <- coxme(Surv(delai, status) ~ (1|trt)+ strata(center) , data)
I am surprised that the output given in R is 3 coefficients for random effects
whereas only 2 dummy
Hi there
I am a Masters student at the University of Stellenbosch. I have been using R
to analyze the data, using the GLS model, of one of my experiments.
The problem that I am having is that whenever I run my model using:
fit.glsmodel1<-gls(Number~as.factor(Season)+as.factor(Depth)+as.factor(Ori
-Original Message-
I have a dataset which for the sake of simplicity has two endpoints. We would like to test
if two different end-points have the same eventual meaning. To try and take an example
that people might understand better:
Lets assume we had a group of subjects who all rece
very helpful!! Thanks a lot.
From: Jim Lemon
Cc: "r-help@r-project.org"
Sent: Monday, April 29, 2013 6:53 PM
Subject: Re: [R] biplot for principal componens analysis
On 04/30/2013 08:24 AM, capricy gao wrote:
>
>
> I did a PCA for my data which has a dim
This comes up regularly. Type "?print.survfit" and look at the comments there under
"value".
Terry T.
- begin included message
Hi,
I'm not sure if this is the proper way to ask questions, sorry if not. But
here's my problem:
I'm trying to do a bootstrap estimate of th
Thanks to you all for your replies. I didn't realize bigmemory is only
available in Unix environments - when I saw
> install.packages('bigmemory')
Installing package into C:/Users/BenC/Documents/R/win-library/3.0
(as lib is unspecified)
--- Please select a CRAN mirror for use in this session -
Hi,
Not sure if this is what you meant.
dat1<- read.table(text="
id t scores scores2 scores3
2 0 1.2 1.4 1.9
2 2 2.3 2.5 2.2
2 3 3.
a) This looks like homework. The Posting Guide clearly indicates that this list
is not for homework help.
b) This is a statistics theory question that happens to use R, not an R
question that happens to be about statistics. Also off-topic per the Posting
Guide... there are other forums for stat
Hi,
May be this helps:
funcName<- function(df1, x){
whatCol=df1[[x]]
print("Got it")
print(whatCol)
}
funcName(df,"ColA")
#[1] "Got it"
#[1] 1 2 3 4 5
funcName(df,"ColB")
#[1] "Got it"
#[1] A B C D E
#Levels: A B C D E
A.K.
>I am trying to extract the 2nd column from a dataframe using a
On Tue, Apr 30, 2013 at 10:40 AM, Fabio Berzaghi wrote:
> Dear All,
>
> Is it possible to add a percentage completion bar to R? I find it
> frustrating when running long calculations that I don't know at what point
> the process is. It would be very helpful so I am not guessing if I should
> wait
On 30/04/2013 10:40, Fabio Berzaghi wrote:
Dear All,
Is it possible to add a percentage completion bar to R? I find it
frustrating when running long calculations that I don't know at what
point the process is. It would be very helpful so I am not guessing if I
should wait a few minutes or hours.
Dear All,
Is it possible to add a percentage completion bar to R? I find it
frustrating when running long calculations that I don't know at what
point the process is. It would be very helpful so I am not guessing if I
should wait a few minutes or hours.
__
Hi,
I am using Extrafont package to install more fonts for my graphs. My
primary graphic tool is ggplot2. I seem to have problem installing the
package, but could not pinpoint where it is.
I try to follow the instruction here:
https://github.com/wch/extrafont
I guess (but am not sure) th
Thanks a lot Joshua.
Regards,
Indrajit
On Tue, Apr 30, 2013 at 1:13 PM, Joshua Wiley wrote:
> Hi Indrajit,
>
> In your first SAS code, change to type=un. cs imposes the (somewhat
> dubious) assumption that the variance of both the intercept and slope are
> equal. If you are using lme4, all ran
Hi all,
I was running the adf test in R.
CODE 1:
adf.test(data$LOSS)
Augmented Dickey-Fuller Test
data: data$LOSS
Dickey-Fuller = -1.9864, Lag order = 2, p-value = 0.5775
alternative hypothesis: stationary
CODE 2:
adf.test(diff(diff(data$LOSS)))
Augmented Dickey-Fuller Te
Hi Indrajit,
In your first SAS code, change to type=un. cs imposes the (somewhat
dubious) assumption that the variance of both the intercept and slope are
equal. If you are using lme4, all random effects in a single block (e.g.,
(1 + month | batch) the 1 = intercept and month = random slope) wil
Hi David,
Sorry if I have explain me wrong!
Mylist<-list(A=data.frame, B=data.frame, C=data.frame, D=data.frame)
I know that this is a malformed command, only was a fast way to descrive
mylist.
Thank you!
-Mensaje original-
De: David Winsemius [mailto:dwinsem...@comcast.net]
Enviado el:
Wow! It's true!!!
Many thanks!!!
montserrat
-Mensaje original-
De: MacQueen, Don [mailto:macque...@llnl.gov]
Enviado el: lunes, 29 de abril de 2013 22:39
Para: De Castro Pascual, Montserrat; r-help@r-project.org
Asunto: Re: [R] rbinding some elements from a list and obtain another list
Hi All,
I am trying to shift from running mixed models in SAS using PROC MIXED
to using lme4 package in R. In trying to match the coefficients of R
output to that of SAS output, I came across this problem.
The dataset I am using is this one:
http://support.sas.com/documentation/cdl/en/statug/630
BINGO!!!
Many thanks arun!
You are a genious!!!
montserrat
-Mensaje original-
De: arun [mailto:smartpink...@yahoo.com]
Enviado el: lunes, 29 de abril de 2013 18:33
Para: De Castro Pascual, Montserrat
CC: R help
Asunto: Re: [R] rbinding some elements from a list and obtain another list
Hi
On 30/04/2013 07:44, Anupam Tyagi wrote:
R-3.0.0 prints an error message at startup. I uninstalled a previous
version and installed this on Windows 8. Uninstall of previous version
deleted all previously installed packages and I forgot to keep a list
of them. How do I correct this error?
"
Error
51 matches
Mail list logo
