Displaying 20 results from an estimated 6000 matches similar to: "how to compare two dataset with same sampling"
2011 Sep 01
2
how to plot a series of data in a dataframe?
hi
i have a dataframe with the name "obsdata"
       V1 V2 V3 V4    V5    V6    V7    V8     V9   V10   V11   V12   V13
1    1001  3 24 12  24.7  44.4  70.1  49.3   33.7   3.0   6.8   2.7    NA
2    1001  3 25  0  70.1  49.3  33.7 138.2  152.5    NA   4.2   6.9  17.5
3    1001  3 25 12  33.7 187.7 286.5 386.7     NA  16.2  46.0  48.8  43.1
4    1001  3 26  0  88.6 129.4    NA    NA    
2010 Jan 29
2
Create matrix with subset from unlist
Hello all,
I'm trying to create a 2x2 matrix, 32 times after unlist() so that I can 
convert the list to matrix. I've looked through the R archive but 
couldn't find the answer. There is what I've done.
 > f <- system("ls *.txt", intern=TRUE)
 > x <- lapply(f, read.table)
 > x
[[1]]
     V1    V2
1 -27.3  14.4
2  29.0 -38.1
[[2]]
     V1   V2
1  14.4
2009 Jul 31
2
how use the subset?
hi ,everyone
I  want subtract some dataset by subset.
>From the help running  help(subset),  ths information is "*subset(airquality,
Day == 1, select = -Temp)* "
while  I running my script written as "*g1data<-subset(errdata, fac>12) *"
,it is wrong with the error information "*subset.default(newerrdata,
fac>12),can not find fac*"
and g1 in read
2010 Oct 04
3
How To Extract Row from A Data Frame
I have a data frame that looks like this:
> print(df)
  V2     V3     V4     V5     V6     V7     V8     V9    V10    V11    V12
1 FN  8.637 28.890 31.430 31.052 29.878 33.215 32.728 32.187 29.305 31.462
2 FP 19.936 30.284 33.001 35.100 30.238 34.452 35.849 34.185 31.242 35.635
3 TN  0.000 17.190 16.460 21.100 17.960 15.120 17.200 17.190 15.270 15.310
4 TP 22.831 31.246 33.600 35.439 32.073
2012 Mar 15
1
eigenvalues of matrices of partial derivatives with ryacas
Hello,
I am trying to construct two matrices, F and V, composed of partial 
derivatives and then find the eigenvalues of F*Inverse(V).  I have the 
following equations in ryacas notation:
 > library(Ryacas)
 > FIh <- Expr("betah*Sh*Iv")
 > FIv <- Expr("betav*Sv*Ih")
 > VIh <- Expr("(muh + gamma)*Ih")
 > VIv <- Expr("muv*Iv")
I
2005 Oct 13
1
expand.grid problem
Hi all,
I want to make all possible combination from dataset below:
  V1   <- c(0,1,2)
  V2   <- c(0,1)
  V3   <- c(0,1)
  V4   <- c(0,1)
  V5   <- c(0,1)
  V6   <- c(0,1,2,3,4,5,6,7,8,9,10,11,12,13,14,15,16,17,18,19,20)
  V7   <- c(0,1,2,3,4,5,6)
  V8   <- c(0,1)
  V9   <- c(0,1)
  V10  <- c(0,1)
  V11  <- c(0,1)
  V12  <- c(0,1)
  V13  <- c(0,1)
  V14 
2008 Sep 05
1
binary order combinations
Dear all!
I have a vector of names
names<-("V1", "V2", "V3",....., "V15")
I could create all possible combinations of these names (of all
lengths) using R:
combos<-lapply(1:15,function(x)
{combn(names,x)
})
I get a list with all possible combinations of elements of 'names'
that looks like this (just the very beginning of it):
[[1]] - the
2009 Apr 04
1
Problem with sample()
Hi,
I'm having a problem using sample() within a function.
Basically I get an error reading:
Error in sample(v, 1, prob = h) : non-positive probability
Can anyone advise me as to the possible origin of this error?
Here is my code
#Discretised Gillespie algorithm function (From SMfSB, D.J. Wilkinson)
gillespied=function (N, T=100, dt=1, ...)
{
        tt=0
	n=T%/%dt
        x=N$M
       
2012 Jul 05
1
Different level set when predicting with e1071's Naive Bayes classifier
Hi!
I'm using the Naive Bayes classifier provided by the e1071 package (
http://cran.r-project.org/web/packages/e1071) and I've noticed that the
predict function has a different behavior when the level set of the columns
used for prediction is different from the ones used for fitting. From
inspecting the predict.naiveBayes I came to the conclusion that this is due
to the conversion of
2015 Aug 22
3
sprintf error: "only 100 arguments allowed"
I'm trying to apply a function defined in the VW R docs, that attemps to
convert a data.table object to Vowpal Wabbit format. In the process i'm
getting the error in printf mentioned in the subject.
The original function is here:
https://github.com/JohnLangford/vowpal_wabbit/blob/master/R/dt2vw.R
Below there is a small example that reproduces the error. The function
works great with
2015 Aug 26
1
sprintf error: "only 100 arguments allowed"
Wouldn't it make sense to have this in the man page?
The 8192-byte limitation for 'fmt' is mentioned but not this one.
Thanks,
H.
On 08/25/2015 02:08 AM, Prof Brian Ripley wrote:
>  From the sources:
>
> #define MAXNARGS 100
> /*               ^^^ not entirely arbitrary, but strongly linked to
> allowing %$1 to %$99 !*/
>
>
>
> On 22/08/2015 04:21, Martin
2010 May 24
2
import data from a csv file
Hi all,
   I have some trouble reading data from a csv file.
   I used command "read.delim("clipboard") to read in the data.
> aalpha.data <- read.delim("clipboard")
> class(aalpha.data)
[1] "data.frame"
> dim(aalpha.data)
[1]  8 25
> colnames(aalpha.data)
 [1] "X"   "V1"  "V2"  "V3"  "V4" 
2006 Jun 19
2
Nested variance-covariance matrix in Multilevel model
Dear R community,
I have trouble implementing a nested variance-covariance matrix in the 
lme function.
The model has two fixed effects called End and logpgc, the response 
variable is the logarithm to base 2 of  Intensity ( log2(Intensity) ) 
and the random effects are called Probe and ProbeNo.
The model has the following nesting structure: A Pixel is nested within 
the ProbeNo,the ProbeNo is
2012 Nov 24
3
How to Label Cases in Scatterplots?
Hi everyone,
i?m trying to graphically display distributions with r and i?m working with
makrodata from the WVS.
the command i?m using is
                                                      
> plot (Makrodata$v11, Makrodata$v12, xlab="Democracy Score Economist",
> ylab= share religious people")
i?m having an additional variable that identifies respectively labels the
2013 Mar 18
2
Confirmatory factor analysis using the sem package. TLI CFI and RMSEA absent from model summary.
Hi R-help,
I am using the sem package to run confirmatory factor analysis (cfa) on some questionnaire data collected from 307 participants. I have been running R-2.15.3 in Windows in conjunction with R studio. The model I am using was developed from exploratory factor analysis of a separate dataset (n=439); it includes 18 items that load onto 3 factors. I have used the sem package documentation
2006 Feb 01
2
sort columns
Hi. 
 
I have a simple (I think) question
 
My dataset have these variables:
 
names(data)
  [1] "v1"   "v2"   "v3"   "v4"   "v5"   "v6"   "v7"   "v8"   "v9"   "v10"
"v11"  "v12"  "v13"  "v14"  "v15"  "v16"  "v17" 
2007 Feb 14
2
[LLVMdev] Linux/ppc backend
Hi Chris,
Chris Lattner wrote:
>> 2) Line 369 of PPCInstrInfo.td, we declare the non-callee saved registers. 
>> However, Linux and Darwin do not have the same set
>> of non-callee saved registers. I don't know how to make the if(isDarwin) test 
>> in here
>>     
>
> Take a look at ARM/ARMRegisterInfo.td for an example of this
I tried to define Defs just
2010 Nov 02
2
count different words in a field
Hi all, 
 
I started to ask this in the other post, but it is off topis...so here it is again.
 
I have a data.frame (created with the helpof this mail list) that looks like this :
 
'data.frame':   22801 obs. of  15 variables:
$ V1 : chr  "HUMUS" "SLABO" "MALO" "SLABO" ...
$ V2 : chr  "IN" "GRANULIRAN"
2005 May 02
1
newbie ifelse matrix question
Hi all,
 
I have time series data in a matrix format 20 rows x 205 columns and have been trying to replace outliers with NA.  My first column contains the outliers threshold (3 Standard deviations) for each row  - here's a bit of the first row
 
   sd3      V1 V2 V3 V4         V5        V6        V7        V8        V9
1  13.03267 1797157 75 84 58  -1.958649  0.048775  2.056198  8.063622 
2010 Oct 13
2
Rmysql - dbWritetable
Dear Rusers,
I am trying to feed my database with data from a file.
But since my file (2010101000.txt) there`s no headers I am facing problem
because the
result data.frame is not separated my columns.
How could I set variables names for each columns in *dbWriteTable*? I have
tried the command below but I don`t know how to fill the *field.type*option.
dbWriteTable(con, "b20101010",