Re: [R] RSQLite: Tracking an error message to its source

2014-01-11 Thread Andrew Hoerner
On Sun, Jan 5, 2014 at 7:12 AM, Hadley Wickham wrote: > > Though I admit that it bothers me that, although the SQLite syntax > > documentation for "CREATE INDEX", here: > > http://www.sqlite.org/lang_createindex.html > > does say the database name is optional, it also says that, if you > include

Re: [R] SE for R-squared

2014-01-11 Thread David Winsemius
On Jan 10, 2014, at 12:32 PM, Troels Ring wrote: > In R package "psychometrics" an estimate of SE of R squared of /sersq <- > sqrt((4*rsq*(1-rsq)^2*(n-k-1)^2)/((n^2-1)*(n+3))) with n sample size, > and k number of parameters if sample size greater than 60 is found. > Does anyone have a formul

Re: [R] Nested for loop help please.

2014-01-11 Thread arun
Hi, I guess you need to replace t() with t.control() or t.sham() output <-vector() for(count in 1:length(pick.a)){  for(count1 in 1:length(pick.d)){  ts1 = sum(t.sham(m.sham,pick.a[count],pick.d[count1])) tc1 <- sum(t.control(m.control,pick.a[count])) output= c(output,ts1/tc1)  }  } m1 <- matri

Re: [R] Levelplots with non-continuous x-axis breaks

2014-01-11 Thread Jim Lemon
On 01/11/2014 09:27 AM, Pachapep wrote: ok, I think I got it. Here is the code I used. Thanks for the help guys.. Hi Pachapep, You can also do it like this: library(plotrix) plot(1:10) # show color.scale(1:50,c(0,1,1),c(0,1,0),c(1,1,0),) gradient.rect(1,4,10,6,c(0,1,1),c(0,1,0),c(1,1,0)) Ji

Re: [R] Problem with barplot

2014-01-11 Thread Jim Lemon
On 01/11/2014 07:41 AM, Clemence Germaine wrote: Hi I’m using bar plot function And I have in my data some zero, but my data are numerics and when i try to use bar plot i have a warning message of « NA was introduced during the automatic conversion » and i can’t see my bar plot. So if somebod

Re: [R] Nested for loop help please.

2014-01-11 Thread William Dunlap
> > for(count in 1:length(pick.a)){ > + for(count in 1:length(pick.d)){ > + ts=sum(t(m.sham,pick.a[count],pick.d[count])) The variable 'count' in the inner loop is overwriting the value of 'count' set in the outer loop. Use different names in the different loops. for(count.a in seq_along(pick.

Re: [R] Exact logistic regression, data frame issue

2014-01-11 Thread peter dalgaard
On 11 Jan 2014, at 19:09 , Razi Zaidi wrote: > Hi there > > I have the dataframe below saved in an object /tst/ > > agegrp bmigrp pep n > 1 1 a 0 2 > 2 2 a 0 2 > 3 3 a 0 0 > 4 1 b 0 47 > 5

Re: [R] Nested for loop help please.

2014-01-11 Thread Rui Barradas
Hello, You are using the function ?t(), matrix transpose, with more than one argument. What is that supposed to do? The error message says that the other two arguments are not used (because they are illegal). And you don't need to transpose a vector to sum its components. Hope this helps, R

Re: [R] Concatenating rows of a column vector

2014-01-11 Thread arun
Hi, Try ?paste() dat1 <- read.table(text="apkkaakafmfffakkannpaaapkacfaapfd kpaaakfkpkfbfakaaofakapkpppfcgaanfpfakaappffak fkpkfbfakaaofakapkpppfcgaanfpfa aakfkpkfbfakaaofakapkpppfc",sep="",header=FALSE,stringsAsFactors=FALSE) res <-  paste(dat1[,1],collapse="") A.K. So I have a 19x1 col

[R] Nested for loop help please.

2014-01-11 Thread Mathew Nagendran
Hi all I am relatively new to R. I am trying to create a nested for loop but i keep getting an error message (unused argument). Can someone help me find out where I am goign wrong? > m.control=c(1.45,9.40,9.96,4.2,1.86,0.2) > m.sham=c(3.39,23.94,23.62,10.08,2.99,1.09) > > t.control=function(m, a

[R] Exact logistic regression, data frame issue

2014-01-11 Thread Razi Zaidi
Hi there I have the dataframe below saved in an object /tst/ agegrp bmigrp pep n 1 1 a 0 2 2 2 a 0 2 3 3 a 0 0 4 1 b 0 47 5 2 b 0 43 6 3 b 0 31 7

Re: [R] How to show Vignettes download link on CRAN package page

2014-01-11 Thread Prof Brian Ripley
On 11/01/2014 17:29, Hai Qian wrote: Hi, This question might be stupid, but I cannot find the answer anywhere. Some packages, like RUnit ( http://cran.r-project.org/web/packages/RUnit/index.html) has a "Vignettes:" download link on that CRAN page. However, my package PivotalR ( http://cran.r-p

Re: [R] Find the package of a class given classname

2014-01-11 Thread Yuan Luo
This is exactly what I am looking for. Thanks all for your help! Best, Yuan On Sat, Jan 11, 2014 at 1:28 AM, Hervé Pagès wrote: > Hi Yuan, > > > On 01/10/2014 06:49 PM, Yuan Luo wrote: > >> How to find the package of a class given classname? >> For example, there is a class called GAlignments,

[R] How to show Vignettes download link on CRAN package page

2014-01-11 Thread Hai Qian
Hi, This question might be stupid, but I cannot find the answer anywhere. Some packages, like RUnit ( http://cran.r-project.org/web/packages/RUnit/index.html) has a "Vignettes:" download link on that CRAN page. However, my package PivotalR ( http://cran.r-project.org/web/packages/PivotalR/) does

Re: [R] Importing data into R

2014-01-11 Thread arun
Hi, May be this helps: url <- "http://genome.ucsc.edu/cgi-bin/hgc?hgsid=358528009&g=htcDnaNearGene&i=uc003qec.4&c=chr6&l=133562494&r=133853258&o=knownGene&hgSeq.promoter=on&boolshad.hgSeq.promoter=0&hgSeq.promoterSize=10759&hgSeq.utrExon5=on&boolshad.hgSeq.utrExon5=0&hgSeq.cdsExon=on&boolshad.hgSe

Re: [R] Find the maximum value of a column in a matrix

2014-01-11 Thread arun
HI, The input and output dataset looks different (typo?) in "profit" column. dat1 <- read.table(text="Place  Prod  Loctn  profit S1 P1 Loc1    55 S1 P2 Loc1   80 S1 P3 Loc2   70 S1 P1 Loc2   60 S2 P2 Loc1   30    S2 P3 Loc1   40 S2 P1 Loc1  

Re: [R] Selecting individuals to maximize the correlation of two variables

2014-01-11 Thread Enrico Schumann
On Fri, 10 Jan 2014, 汤靖 writes: > Hi, > Maybe it is not directly related to R but sine many are statistical experts > so I post it here for help: > > I have two variables (say x and y) of length n. Now the cor(x,y) is close to > 0. I need to find the subset in {1,.. n} so that the correlation b

Re: [R] looping through 3D array

2014-01-11 Thread arun
Hi Alex, Not sure if this is what you wanted. length(res) #from the previous 'example' using ##indx <- combn(dim(results)[1],2) #[1] 45 mat1 <- matrix(0,10,10)  mat1[lower.tri(mat1)] <- res  mat1[upper.tri(mat1)] <- res A.K. On Saturday, January 11, 2014 12:22 AM, alex padron wro

Re: [R] Problem with barplot

2014-01-11 Thread Marc Girondot
Le 10/01/2014 21:41, Clemence Germaine a écrit : Hi I’m using bar plot function And I have in my data some zero, but my data are numerics and when i try to use bar plot i have a warning message of « NA was introduced during the automatic conversion » and i can’t see my bar plot. So if somebod