I have the following code for the minimum and maximum of my prediction interval
> y.down=lines(x[x.order], set1.pred[,2][x.order], col=109)
> y.up=lines(x[x.order], set1.pred[,3][x.order], col=109)
domain=min(x):max(x)
polygon(c(domain,rev(domain)),c(y.up,rev(y.down)),col=109)
It doesnt seem t
Dear all,
I would like to add mixed effects in a multinomial model and I am trying
to use MCMCglmm for that.
The main problem I face: my data set consits of a trapping data set,
where the observation at eah trap (1 or 0 for each species) have been
aggregated per traplines. Therefore we have
Yes,
I have a row and column header. I need to save the row and column header.
Here is the example of the data
*Whirr_127.csv*
WHIRR.127 WHIRR.128 WHIRR.129 WHIRR.137 WHIRR.139 WHIRR.140 WHIRR.141
WHIRR.145 WHIRR.146 WHIRR.147
Adrian Cole 0 0 1 0 0 0
You need to learn to use the str() function to debug these kinds of problems.
And you need to learn to provide data with your code (perhaps with the dput
function) so others can reproduce your problem if you want help.
I would guess that you have non-numeric data in your files.
-
Hi,
I have two matrix from two different .csv files.
#load .csv files
a <- as.matrix(read.table("Whirr_127.csv", header=T, sep=",", row.names=1))
b <- as.matrix(read.table("Files_Whirr_127.csv", header=T, sep=",",
row.names=1))
a
b
I managed to do transpose to 'b' without any error.
transpose_ta
Hello
I have data like this
x1 x2 x3 x4 x5
I want to create a matrix similar to a correlation matrix, but with the
difference between the two values, like this
x1 x2 x3 x4 x5
x1 x2-x1x3-x1 x4-x1 x5-x1
x2 x3-x2 x4-x2 x5-x2
x3
I've noticed that ff uses the unzip utility available on it's host OS to
load datasets via ffload. It seems to work fine in linux, but when I try
to use the package in Windows (hence dling Windows unzip utils) I get 2
errors, one telling me that the options aren't being passed to unzip
correctly a
>Le mercredi 12 septembre 2012 à 07:08 -0700, Tim Hesterberg a écrit :
>> One approach is to bootstrap the vector 1:n, where n is the number
>> of individuals, with a function that does:
>> f <- function(vectorOfIndices, theTable) {
>> (1) create a new table with the same dimensions, but with the
The map() function plots polygons. Since several states have noncontiguous
parts, there are actually 63 polygons plotted:
> statepolys <- map("state", namesonly=TRUE, plot=FALSE)
> statepolys
The following commands produce the map after merging the color values in
dState to the list of polygons:
On 12/09/2012 21:52, Ista Zahn wrote:
Hi Andras,
The Sweave manual is at
http://www.statistik.lmu.de/~leisch/Sweave/Sweave-manual.pdf. There
you can learn about the echo option any many other wonderful things.
That's for R 2.7.1: there have been quite a few changes since and a
current manual
Many thanks Paul, that is a nifty way to do it.
Regards
Ben
Paul Murrell wrote:
Hi
On 07/09/12 09:35, Benjamin Tyner wrote:
Update: seems one way to skin this cat is to add
gp = gpar(fontfamily="mono")
to the viewport() call itself. If anyone has any suggestions for a
robust way to extra
Even this one-liner shows Michigan with 2 colors.
map("state", col=1:3, fill=TRUE)
This link might help:
http://tolstoy.newcastle.edu.au/R/devel/06/08/6354.html
Kevin
On Wed, Sep 12, 2012 at 10:38 AM, Josh Browning
wrote:
> I'm having issues when trying to do a filled State plot using the maps
I have not found an approach to time series analyses of data that include
left-censored observations in the resources I have here. If there is an
approach please point me to some references.
Thanks,
Rich
__
R-help@r-project.org mailing list
https://
This is close, but it does quote the header names, but does produce
the same dataframe when read back in:
> RAW.API <-
> structure("id,event_arm,name,dob,pushed_text,pushed_calc,complete\n\"01\",\"event_1_arm_1\",\"John\",\"1979-05-01\",\"\",\"\",2\n\"01\",\"event_2_arm_1\",\"John\",\"2012-09-02\
Dear R experts,
I'm reading data from an online database via API and it gets delivered in this
messy comma separated structure,
> RAW.API <-
> structure("id,event_arm,name,dob,pushed_text,pushed_calc,complete\n\"01\",\"event_1_arm_1\",\"John\",\"1979-05-01\",\"\",\"\",2\n\"01\",\"event_2_arm_1\
Thanks you all,
With your recomendations I built my solution as follows:
:> writeLines(c("uno dos tres", "cuatro cinco", "seis"), "tempfile.txt")
:> o <- pipe("wc < tempfile.txt", open="r")
:> readLines(o)
:[1] " 3 6 31"
best regards,
--Sergio.
__
On 12-09-12 4:52 PM, Ista Zahn wrote:
Hi Andras,
The Sweave manual is at
http://www.statistik.lmu.de/~leisch/Sweave/Sweave-manual.pdf. There
you can learn about the echo option any many other wonderful things.
It is also available within R as
vignette("Sweave")
in current versions.
Duncan M
If you want to solve it _numerically_ there are many options within R:
if you want to solve symbolically, you'll probably need to interface
to some symbolics software, e.g., Ryacas.
Cheers,
Michael
On Wed, Sep 12, 2012 at 11:08 PM, Bert Gunter wrote:
> Well, then ...
>
> I've cc'ed this to R-Hel
Well, then ...
I've cc'ed this to R-Help so that you can get other opinions to that I
express below.
R is for data analysis, statistics, and graphing data. R **can** do
(natively) some symbolic mathematics/algebra; and there are R packages
(Ryacas for one) that provide interfaces to symbolic math
It is my normal practice to install R libraries without root. Just use your own
library directory instead of the system library.
---
Jeff NewmillerThe . . Go Live...
DCN:Basic
Dear Samantha
I suggest that you directly use the systemfit() command rather than
the wrapper function zelig(). You can impose parameter restrictions
with the arguments "restrict.matrix", "restrict.rhs", and
"restrict.regMat". You can use argument "restrict.matrix" to specify
the restrictions nume
Hi Andras,
The Sweave manual is at
http://www.statistik.lmu.de/~leisch/Sweave/Sweave-manual.pdf. There
you can learn about the echo option any many other wonderful things.
Best,
Ista
On Wed, Sep 12, 2012 at 1:43 PM, Andras Farkas wrote:
> Dear All,
>
> I am working with a Sweave document to be
``commandArgs`` is the function in base R that lets you access arguments
passed to Rscript, does not allow any fancy parsing of options but if you
are not able to install optparse it is better than nothing.
- Trevor
On Wed, Sep 12, 2012 at 1:32 PM, Tom Roche wrote:
>
> https://stat.ethz.ch/pipe
Homework? We don't do homework here.
-- Bert
On Wed, Sep 12, 2012 at 11:56 AM, pial das wrote:
> Hello,
>
> I am brand new in this site as well as in R. Please accept my early apology
> if it's not the right place to ask such a question.
>
> I have the following equation:
>
> b=[2(1-a)(1+2a)^0.
Hello,
I am brand new in this site as well as in R. Please accept my early apology if
it's not the right place to ask such a question.
I have the following equation:
b=[2(1-a)(1+2a)^0.5]/(1+3a)
How can I solve the above equation for 'a' using R?
Thanks in advance.
KD
[[alternative H
Hello all,
I am following some standard code from Zelig manual when using a SUR
(Seemingly Unrelated Regression Model) to constrain parameters across
equations.
Please see code below:
setwd("C:/Research/Economics/SUR_FX/Model")
# Seemingly Unrelated Regression
# Load our library.
library(Zelig)
Hi,
If you want to just count the words in R, try this:
vec1<-c("uno dos tres", "cuatro cinco", "seis")
#get individual word count within quotes
res1<-unlist(lapply(strsplit(vec1, " "),length))
res1
#[1] 3 2 1
#get whole word count
length(unlist(strsplit(vec1, " ")))
#[1] 6
- Original Me
Dear All,
I am working with a Sweave document to be converted into PDF using Rstudio. It
seems to me that my R code will also show up after conversion, which I would
like not to happen. Is there a way to specify a command that I could place on
the beggining and at the end of the R code that wo
Hello Amelie,
I don't have an answer to your question, but I just wanted to point out
this page I noticed recently (
http://hlplab.wordpress.com/2009/05/07/multinomial-random-effects-models-in-r/),
which might be helpful.
I'm also interested in figuring out how to do a multinomial glmm, so if you
I'm having issues when trying to do a filled State plot using the maps()
package. It seems that the color values aren't being plotted correctly,
and I'm not sure why...
I've started debugging it, but I'm stuck. I've attached the code that I've
got so far. The for loop plots states 1:i. What I'
I'm trying to use a data base from SPSS and get all data
representing the true data. I would like to use a variable as expansion
data. In our data base we have the variable with the factor expansion,
but we have not idea how to set it in R.
Any idea?
Thanks in advance!
José
[[alternati
https://stat.ethz.ch/pipermail/r-help/2012-September/323551.html (slightly
edited)
>> how to structure an R file such that it can be both
>> 1. used as a script via, e.g., (from OS commandline)
>> $ Rscript foo.r bar=baz
>> 2. imported and called as a function via, e.g. (from R commandline)
>
I encountered this problem when I was learning nonlinear regression with R :
there is this function "plotfit" in package "nlstools",which "displays
a superimposed plot of the dependent variable versus one the
independent variables together with the fitted model.",but when I
execute the example give
Hello again
I tried to do everything how described on the page:
http://rwiki.sciviews.org/doku.php?id=getting-started:installation:android
So, i pushed the three files in the folder /sdcard/gcc
Then i tried to extract the files with an android terminal with the
commands:
cd /sdcard/gcc
tar xjf
suddenly a few days ago it takes about a minute for local help to come up .
i type :
?hist
and get
starting httpd help server ... done
but only after a minute does
http://127.0.0.1:15686/library/graphics/html/hist.html
appear .
Lee De Cola, PhD, MCP.
DATA to Insight
ldec...@comcast.n
thanks all for the prompt answer.
Yes, I need to go through the R tutorial rather than learning a snippet of
codes from googling.
On Wed, Sep 12, 2012 at 11:32 AM, Rui Barradas wrote:
> Hello,
>
> The input must be a matrix, not a list (or its special case data.frame).
>
> Var <- read.table(tex
Hi all,
Sorry about posting a really novice question.
I was able to run rcorr after converting the list to a matrix by your help.
I'm though wondering if there is any way to find out an exact p value as
the output only gave me 0 for P value as shown below.
I've added options(digits=10), which doesn
Le mercredi 12 septembre 2012 à 07:08 -0700, Tim Hesterberg a écrit :
> One approach is to bootstrap the vector 1:n, where n is the number
> of individuals, with a function that does:
> f <- function(vectorOfIndices, theTable) {
> (1) create a new table with the same dimensions, but with the coun
Hello,
You could post to R-Help, maybe it will of use to others. And since it's
a follow-up I'll post it there.
To get more (or less) digits use the respective argument to print().
Var <-
structure(list(D.Prime = c(0.17234, 0.14399, 0.14626, 0.035147,
0.058957), T.statistics = c(4.9268, 2.892
On Wed, Sep 12, 2012 at 12:40:54PM -0500, Julio Sergio Santana wrote:
> Hi,
> I'm trying to use pipes in R. By now, I could launch the linux command "wc"
> (to count words from a text), but
> I don't know how to capture the results, say in a vector of chars...
> Here is the R code I'm trying:
>
>
On 12/09/2012 1:40 PM, Julio Sergio Santana wrote:
Hi,
I'm trying to use pipes in R. By now, I could launch the linux command "wc"
(to count words from a text), but
I don't know how to capture the results, say in a vector of chars...
Here is the R code I'm trying:
:> f <- pipe("wc", open="w")
:>
No, this is not what I want, "wc" is just an example. What I want is to
shut a process to do a task outside R, and to
get its results back in R
On Wed, Sep 12, 2012 at 1:07 PM, arun wrote:
> Hi,
>
> If you want to just count the words in R, try this:
> vec1<-c("uno dos tres", "cuatro cinco", "
Hi,
I'm trying to use pipes in R. By now, I could launch the linux command "wc"
(to count words from a text), but
I don't know how to capture the results, say in a vector of chars...
Here is the R code I'm trying:
:> f <- pipe("wc", open="w")
:> writeLines(c("uno dos tres", "cuatro cinco", "seis")
(Sorry about the double post; it seems the last time I posted this it
didn't show up correctly because of some character encoding issue.)
Hello everyone,
I am testing a model in which I have a two-level factor (let's call it
First [1, 2]) nested under a four-level factor (let's call it Second [A,
Thanks for including an example that could be copied and pasted.
However TkPredict and Predict.Plot both need at least one numeric (not
factor) predictor. So I added another column called x that was just
1:5 to the sample data frame and included it in the model.
Here are a couple of approaches us
Hello,
The input must be a matrix, not a list (or its special case data.frame).
Var <- read.table(text="
D.PrimeT.statistics
11.7234e-01 4.926800
21.4399e-01 2.892000
31.4626e-01 2.642800
43.5147e-02 1.112400
55.8957e-02 2.723700
", header=TRUE)
# l
On Wed, Sep 12, 2012 at 9:02 AM, Geophagus wrote:
> Hi @ all and thanks a lot for your advises.
> So far I don't know the function dput. Now I read something about it and I
> have the following result:
>
>> TEST
> year value state
> 1 200715FL
> 2 200816FL
> 3 200814TX
>
Have you read (relevant sections) of "An Introduction to R" to gain
some basic understanding of how R works? If not, please do so before
further posting.
Briefly, Var is a data frame and rcorr wants a matrix.
-- Bert
On Wed, Sep 12, 2012 at 8:04 AM, Jason Love wrote:
> Hello,
> I'd like to tes
On Wed, Sep 12, 2012 at 4:04 PM, Jason Love wrote:
> Hello,
> I'd like to test a significance of two variables in their correlation using
> rcorr, which gave me an error of format incompatibility.
> Below are the lines that I typed in the R window and let me know if anyone
> knows how to resolve t
Hello,
I'd like to test a significance of two variables in their correlation using
rcorr, which gave me an error of format incompatibility.
Below are the lines that I typed in the R window and let me know if anyone
knows how to resolve this.
Var=read.csv("03apr10ab_corr_matrix_in_overlaps.csv",hea
On Wed, Sep 12, 2012 at 1:24 PM, Simon Fry wrote:
> Hello again
>
> I tried to do everything how described on the page:
>
> http://rwiki.sciviews.org/doku.php?id=getting-started:installation:android
>
> So, i pushed the three files in the folder /sdcard/gcc
> Then i tried to extract the files with
On Wed, Sep 12, 2012 at 8:36 AM, Vincy Pyne wrote:
> Dear Mr Weylandt and R helpers,
>
> Thanks a lot for your suggestion. Unfortunately the "return" statement in
> my original R code returns me different results which are obtained after
> processing the function I have constructed.
>
> My requi
Hello everyone,
I am testing a model in which I have a two-level factor (let's call it
First [1, 2]) nested under a four-level factor (let's call it Second [A, B,
C, D]). I have used the following model to get coefficients representing
whether, for each level of Second, there is a significant diff
Thanks for the reply. That was the problem. My getwd() was
"C:/Windows/system32". I changed and now it works!
Best regards,
João Fadista, Ph.D.
Post Doc
Lund University Diabetes Centre
CRC, Malmö University Hospital
Jan Waldenströms gata 35, building 60, level 13
SE-205 02 Malmö, Sweden
Tel: +46
One approach is to bootstrap the vector 1:n, where n is the number
of individuals, with a function that does:
f <- function(vectorOfIndices, theTable) {
(1) create a new table with the same dimensions, but with the counts
in the table based on vectorOfIndices.
(2) Calculate the statistics of
I don't think you can with the current version of ggplot2.
There appears to be development work on the issue but It's beyond my level of
knowledge to implement the developmental programming althought it looks good.
https://groups.google.com/forum/?fromgroups=#!topic/ggplot2/mBdq01BxoKE seems
Are you sure that you have write permission to the directory that R is trying
to write to?
John Kane
Kingston ON Canada
> -Original Message-
> From: joao.fadi...@med.lu.se
> Sent: Wed, 12 Sep 2012 12:37:42 +
> To: r-help@r-project.org
> Subject: [R] SNPRelate package error
>
> Dear
Dear all,
I am using the R package SNPRelate but I found an error when I run the
following command. Do you know what might be the problem? Thanks in advance.
> vcf.fn <- system.file("extdata", "sequence.vcf", package="SNPRelate")
> snpgdsVCF2GDS(vcf.fn, "test.gds")
Start snpgdsVCF2GDS ...
Hello,
Two corrections:
1. To convert from Celsius to Farenheit it's degC*1.8 + 32, not degC*1.4
+ 32
2. The axis tick marks should be in the same _places_, but with
different _labels_.
degC<-seq(-40,120,by=20)
plot(degC,1:9)
axis(3, at=degC, labels=degC*1.8+32) # here
mtext("degF",3,line=2
Dear Meghan,
you can do this with rma() in metafor, but you will have to set up the loop
yourself. Here is an example. First, I need to simulate some data:
library(metafor)
library(MASS)
k <- 40
On 09/12/2012 09:29 AM, Jonas Stein wrote:
Hi,
how can i plot two different x axis in a ggplot2 qplot?
I want to plot Farenheit and Celsius in one diagram.
x1:Farenheit x2:Celsius
Hi Jason,
Here's a basic example of how to do this sort of thing.
degC<-seq(-40,120,by=20)
plot(degC,1:9)
axis(3,
Hi
and please do not use HTML messages, they get scrambled during processing.
Plain text is preferred.
Regards
Petr
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
> project.org] On Behalf Of PIKAL Petr
> Sent: Wednesday, September 12, 2012 10:49 AM
Hi
This is the error I get with your code and this prevents us to give you
reasonable help.
Error in eval.quoted(.variables, data) : object 'My_data' not found
Petr
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
> project.org] On Behalf Of Vincy Py
Hi
It is good to cc to R help. It would be even better if you provide toy data to
illustrate what do you have and what do you want. Sometimes you will find
solution when preparing such data yourself without need for posting to the
rhelp.
I am still not sure what do you really want to achieve s
Hi @ all and thanks a lot for your advises.
So far I don't know the function dput. Now I read something about it and I
have the following result:
> TEST
year value state
1 200715FL
2 200816FL
3 200814TX
4 200813TN
5 200918FL
6 200915OH
7 201019
Dear all,
I would like to add mixed effects in a multinomial model and I am trying
to use MCMCglmm for that.
The main problem I face: my data set consits of a trapping data set,
where the observation at eah trap (1 or 0 for each species) have been
aggregated per traplines. Therefore we have
Dear Mr Weylandt and R helpers,
Thanks a lot for your suggestion. Unfortunately the "return" statement in my
original R code returns me different results which are obtained after
processing the function I have constructed.
My requirement for storing the product-wise random numbers is just a pa
On 11/09/2012 13:09, Duncan Murdoch wrote:
On 12-09-10 9:21 PM, Peter Dunn wrote:
Hi all
I've been getting some emails from the R package maintainers that I need
to update some code in a CRAN packge that uses FORTRAN, to comply with
(not so recent) changes. I've been a little busy...
I'm havi
Hi,
how can i plot two different x axis in a ggplot2 qplot?
I want to plot Farenheit and Celsius in one diagram.
x1:Farenheit x2:Celsius
kind regards,
--
Jonas Stein
__
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
P
Hi
> -Original Message-
> From: r-help-boun...@r-project.org [mailto:r-help-bounces@r-
> project.org] On Behalf Of Vincy Pyne
> Sent: Wednesday, September 12, 2012 8:51 AM
> To: r-help@r-project.org
> Subject: [R] How to append the random no.s for different variables in
> the same data.fra
On Wed, Sep 12, 2012 at 7:51 AM, Vincy Pyne wrote:
> Dear R helpers,
>
> (At the outset I sincerely apologize if I have not put forward my following
> query properly, though I have tried to do so.)
>
>
> Following is a curtailed part of my R - code where I am trying to generate
> say 100 random
Dear R helpers,
(At the outset I sincerely apologize if I have not put forward my following
query properly, though I have tried to do so.)
Following is a curtailed part of my R - code where I am trying to generate say
100 random no.s for each of the products under consideration.
library(plyr
72 matches
Mail list logo