similar to: Trend lines on a scatterplot matrix

I want to generate a lattice plot of a multiple linear regression. I'm using the code: xyplot(y ~ x1 + x2 | status, data=datam, xlab="Peak separation",ylab="G/W",main="G/W vs Fuzzy peak separation: Threshold=1.8", groups=Fuzzy.gw.t.score>1.8, subset=(status %in% c("control","patient","sibling")),

Dear List, Some independent variable were missing in calculation using lm and glm (gaussian). (X= Y1+Y2+…..+Y16, Independent number: 16 variable) However, those variables did work well in cor(X, Y) respectively. str(dataframe) was also run to ensure that the variables were all numbers. Moreover, the missing variables were different in lm and glm. In lm, 3 factors were not taken into

Dear List, I tried to do cca based on species data and environmental variables (formula instead of community data). However, there was an error saying row sums must be >0. I searched the previous related messages but found few solutions. Please kindly help and thank you in advance. code This is vegan 1.17-3 Warning message: package 'vegan' was built under R version 2.10.1

Hello, I am using package MuMIn to calculate AIC for a full model with 10 explanatory variables. Thanks in advance in sharing your experience. Q1 In the AIC list of all models, each model is differentiated by model number. Please kindly advise if it is possible to find the corresponding explanatory variable(s) for the model number. Q2 error message I tried to display sub-model with only

In R, I'm imported a data frame of 2,321,123 by 4 called "dataF". I converted the data frame "dataF" to a matrix dataM <- as.matrix(dataF) Does R have an efficient routine to treat the special elements that contain "inf" in them. For example, can you separate the rows that have "inf" elements from the matrix into a separate matrix without

Hello I have a table describing butterfly range traits. It is composed of three columns as below Species name range size (X) latitude of range midpoint (Y) There are 11 kinds of butterflies. Each has its range size, and the latitude of each range midpoint ranges from 9 to 19. I would like to have the average range size of every degree of latitude. For example, the average range

I have some data consisting of multiple trials of an experiment with different values of an independent variable. If I run R> plot(var,result) I get a scatterplot of result versus the independent variable var. If I run R> plot(as.factor(var),result) I get a boxplot of the distribution of result for each value of var. In this plot, each boxplot is labeled by the corresponding value of

Hi, I'm trying to save heatmaps as PDFs. However, the PDF version of the heatmaps (Heatmap_CAFvsTNF_run2.pdf) is blurred when compared to its counterpart, which was saved manually by using the software "Grab" (Heatmap_CAFvsTNF_run2.tiff). -----R code-------- sample_output <- "stroma_run2" filename <-

Among others, datam contains the columns: logconc, tm, dose, subj, bilirubin. None of these are factor variables. The following compartment models work (the first still has not converged after 100 interations): res1 <- nlme(logconc~p2+p3+log(dose/(exp(p1)-exp(p2))* (exp(-exp(p2)*tm)-exp(-exp(p1)*tm))),start=list(fixed=c(5,-2,-0.1)), fixed=list(p1+p2+p3~1),control=list(maxIter=100),

On Thu, Dec 14, 2017 at 1:58 PM, Berend Hasselman <bhh at xs4all.nl> wrote: > >> On 14 Dec 2017, at 19:36, lily li <chocold12 at gmail.com> wrote: >> >> Hi R users, >> >> I have a question about reading from text files. The file has the structure >> below: >> >> Time Column1 Column2 >>

Hi - How can I add different notes in different panels? My data looks like ID Dose Visit Value 1 1 0 0.5 -6.5802e-02 2 1 0 1.0 2.4085e-01 3 1 0 1.5 -2.2907e-01 4 1 0 2.0 2.4074e-01 ... ... ... ... ... 270 45 30 3.0 -8.1316e-01 271 46 30 0.5 -2.0786e-01 272 46 30 1.0 -2.9336e-01 273 46 30 1.5 -5.5657e-01 274 46 30 2.0 -2.1659e-01

Hello I have a problem on keeping the format when I export a matrix file with the write.table() function. When I import the data volcano from rgl package it looks like this in R: > data[1:5,] [,1] [,2] [,3] [,4] [,5] [,6] [,7] [,8] [,9] [,10] [,11] [,12] [,13] [,14] [1,] 100 100 101 101 101 101 101 100 100 100 101 101 102 102 [2,] 101 101 102 102 102 102 102

I use the method, df$Time = as.POSIXct(df$Time), but it has the warning message: Error in as.POSIXlt.character(x, tz, ...) : character string is not in a standard unambiguous format On Thu, Dec 14, 2017 at 1:31 PM, MacQueen, Don <macqueen1 at llnl.gov> wrote: > In addition to which, I would recommend > > df <- read.table("DATAM", header = TRUE, fill = TRUE, >

Hi R users, I have a question about reading from text files. The file has the structure below: Time Column1 Column2 01.01.2001-12:00:00 01.01.2001-24:00:00 12 11 01.02.2001-12:00:00 13 10 01.02.2001-24:00:00 11 12 01.03.2001-12:00:00 15 11 01.03.2001-24:00:00 16 10 ... I

Dear All- I am trying to merge two data files - they have different date formats and different times zones. I need to match up the date/time of the datasets and then invoke a conditional statement, such as: if dataC$mph is >= 12 then keep dataM$co23 for the corresponding time/date stamp. snippets of data files: *dataC.txt* LST in mph Deg DegF DegF2 % volts Deg

In addition to which, I would recommend df <- read.table("DATAM", header = TRUE, fill = TRUE, stringsAsFactors=FALSE) and then converting the Time column to POSIXct date-time values using as.POSIXct() specifying the format using formatting codes found in ?strptime because the times are not in the POSIXct default format. This example might indicate the idea: >

> On 14 Dec 2017, at 19:36, lily li <chocold12 at gmail.com> wrote: > > Hi R users, > > I have a question about reading from text files. The file has the structure > below: > > Time Column1 Column2 > 01.01.2001-12:00:00 This line does not contain 3 elements; only one. You'll have to fix that line. Delete it, prepend it with a

Does anyone know if load balancing and DNAT work well together? I know that load balancing and NAT do not, but what about a simple port forward? I can''t apply Julian Anastasov''s patches, because they don''t work with PPTP patches. :/ Anyhow, a simple: iptables -t nat -A PREROUTING -i eth0 -p tcp --dport xxx -j DNAT --to yyy:xxx iptables -I FORWARD -i eth0 -d yyy -p

I'm more interested in whether or not there are any projects in the works that fuse the tools together to build a better IDE experience. I forget which presentation it was, but one of the developers (I *think* it was the gentleman that presented something on CLANG at the dev's conference 2008) made a bit of a point by saying that LLVM, unlike GCC has been built from the ground up with

How do I add to a trellis plot the best fit line from a robust fit? I can use panel.lm to add a least squares fit, but there is no panel.rlm function. -- Alan S Barnett <asb at mail.nih.gov> NIMH/CBDB