similar to: Unconsistent behaviour of function cor()

Hi r-users,   I would like to extract the data that match.  Attached is my data: I'm interested in matchind the value in column 'intg' with value in column 'rand_no' > cbind(z=z,intg=dd,rand_no = rr)             z  intg rand_no    [1,]  0.00 0.000   0.001    [2,]  0.01 0.000   0.002    [3,]  0.02 0.000   0.002    [4,]  0.03 0.000   0.003    [5,]  0.04 0.000   0.003    [6,] 

Hi there, I am doing multiple comparisons for data that is not normally distributed. For this purpose I tried both functions kruskal{agricolae} and kruskalmc{pgirmess}. It confuses me that these functions do not yield the same results although they are doing the same thing, don't they? Can anyone tell my why this happens and which function I can trust? kruskalmc() tells me that there are no

Hi, Using the example data from library(gvlma) library(gvlma) data(CarMileageData) CarMileageNew <- CarMileageData[,c(5,6,3)] ?lst1 <- list() ?y <- c("NumGallons", "NumDaysBetw") ?for(i in seq_along(y)){ ?lst1[[i]] <- gvlma(lm(get(y[i])~MilesLastFill,data=CarMileageNew)) ?lst1} pdf("gvlmaplot.pdf") ?lapply(lst1,plot) dev.off() You could also use

Dear R helpers, I have been using the loadings function from the multiv library and I get the typical output (see below). When I try to export these results to a file using a write.table() I get the following error message "Error in as.data.frame.default(x[[i]], optional = TRUE) : can't coerce loadings into a data.frame" Any idea why write.table is doing that and any

Sadly you failed to set your email program to send plain text and the data is corrupted at my end. I also think you need to reduce the size of the data set... the intent here is to increase your understanding, not debug your particular analysis. I will say that I am having a very challenging time understanding what you are trying to accomplish though. What are the equations that you think need

OK, Let's try this again! Here is the reproducible script; it is long because I had to copy the panel dataset here. My question is related to systemfit; I don't know how to get the result for the entire panel. #Reproducible script Empdata<- read.csv("/Users/ngwinuiazenui/Documents/UPLOADemp.csv") View(Empdata) install.packages("systemfit")

Dear R users I have a big matrix like 6021 1188 790 290 1174 1015 1990 6613 6288 100714 6021 1 0.658 0.688 0.474 0.262 0.163 0.137 0.32 0.252 0.206 1188 0.658 1 0.917 0.245 0.331 0.122 0.148 0.194 0.168 0.171 790 0.688 0.917 1 0.243 0.31 0.122 0.15 0.19 0.171 0.174 290 0.474

I have a problem with the textmatrix() function of the LSA package whenever I specify 'removeNumbers=TRUE'. The data for the function are stored in a directory LSAwork which consists of a series of files that houses the text in column form. As long as removeNumbers = FALSE or it is not present the textmatrix function works just fine. The error message I get seems to suggest it is

Hi all I am trying to generate a normal unbalanced data to estimate the coefficients of LM, LMM, GLM, and GLMM and their standard errors. Also, I am trying to estimate the variance components and their standard errors. Further, I am trying to use the likelihood ratio test to test H0: sigma^2_b = 0 (random effects variance component), and the t-test to test H0:mu=0 (intercept of the model Yij = mu

Que tal comunidad, hace dias que estoy en un problema que no puedo resolver. Resulta que debo confeccionar un grafico de barras (columnas) que debe incluir barras de error. La unica forma en que me resulto es con barchart() de la libreria Hmisc, y haciendo un truco un poco complicado para mi nivel de usuario autodidacta. Ahora que ya lo tengo (despues de semanas de busqueda) resulta que me

I currently work on a draft of an aquatic bioassessment. The conditions tested are the following: ER river water T dechlorinated water control 0.5 + 0.5mg / L of malate T + 1 dechlorinated water control + 1g / L of malate T ED dechlorinated water control SED + ER + river water sediment SED ED + sediment + water dechlorinated. It is the result of AChE in muscle (fillet of fish). The production of

Gracias Carlos !!!  ... alguna idea con la libreria lattice ? es solo que tengo que hacer 3 graficos diferentes y el aspecto de los mismos quedara diferente con lattice y ggplot. Voy a ver si me resulta todo en ggplot. Saludos y gracias de nuevo, Eric. On 22/09/18 12:08, Carlos Ortega wrote: > Hola, > > ggplot tiene un /geom /justamente para esto... > >

Hi all, I'm seeing this failure on my PPC G4 box running TOT with llvm-gcc 4.2. Is anyone else seeing this? I'm sure it's related to the byval stuff that's recently gone into LLVM. I'm attaching the output of this command: $ llvm-gcc -emit-llvm -O3 -S -o - -emit-llvm /Users/wendling/llvm/ llvm.src/test/CFrontend/2008-01-25-ByValReadNone.c As you can see in it, there

Dear all, I am attempting to explain patterns of arthropod family richness (count data) using a regression model. It seems to be able to do a pretty good job as an explanatory model (i.e. demonstrating relationships between dependent and independent variables), but it has systematic problems as a predictive model: It is biased high at low observed values of family richness and biased low at

Hello! I have this problem: I want to create a Venn's diagram with three lists of genes'names. The first is all the genome, the second a subset of it comprising all mitochondrial genes, and the third including all genes that correlate with a given gene. This is what I do: > library(gplots) > A<-read.delim("F:/.../mito genes just names.txt") >

Hello R users, Can any one please help me to find a way to build phylogenetic tree by R program from any distance matrix. Suppose I have a data like : MATRIX [1] '1' 0.0 [2] '2' 0.071 0.0 [3] '3' 0.036 0.286 0.0 [4] '4' 0.429 0.75 0.714 0.0 [5] '5' 0.679 0.179 0.214 0.536 0.0 [6] '6' 0.893 0.929 0.964 0.464 0.357 0.0 [7]

Hello, I’m trying to follow the syntax of a script from a journal website. In order to create a regression formula used later in the script, the regression matrix must have column names “X1”, “X2”, etc. I have tried to assign these column names to my matrix ScoutRSM.mat using a for loop, but I don’t know how to interpret the error message. Suggestions? Thanks, Paul

Hi, all, I'm a newbie to ext3 file system, but what a pity if ext3 could not shrink after containing files and subdirectories get deleted. If the ext3 directory could not shrink, then another natural question is: can the deleted directory entries be overwritten by new files/subdirs? The following is an example to detail my question: Suppose a directory named myDir hold 3 files a, b, c.

the following the the lower.tri matrix in a file named luxry.car and i want to read it in R as a lower.tri matrix.how can i do? i have try to use help.search("read"),but no result what i want. 1.000 0.591 1.000 0.356 0.350 1.000

Hi All: Greetings, and best wishes from the festive times here at Kolkata, India -- the time of Durga Puja celebrations. I seek your advice as I try plotting lines for my data. The problem: I have created a dataframe that looks like this (name: myFrame): lowest second third fourth highest significance INAS 0.107 0.115 0.123 0.115 0.166 0.000 MMA 0.091 0.107 0.115