search for: datam

Displaying 20 results from an estimated 29 matches for "datam".

Did you mean: data
2010 Aug 17
2
Independent variables omitted in lm and glm
...nsideration. In glm, only one of them was omitted. (attached 2 shots) Please kindly advise whether further info is in need to solve the issue. Also, if similar problems have been encountered, please kindly share your experience. Thank you. Elaine Code rm(list=ls()) library(MuMIn) datam <-read.csv("c:/migration/Mig_ratio_20100817.csv",header=T, row.names=1) dim(datam) datam[1,] # original regression model (16 indep. variables) Mig.lm <-lm(datam$SummerM_ratio~datam$temp_ran+datam$temp_mean+datam$temp_max+datam$temp_min+datam$evi_ran+datam$evi_mean+datam$ev...
2010 Jul 17
2
cca in vegan (formula instead of community matrix data)
...ever, 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 rm(list=ls()) library(vegan) datam <-read.csv("c:/migration/M_R_20100707_winterM.csv",header=T, row.names=1) dim(datam) datam[1,] m.cca <- cca(datam$WinterM_ratio~datam$topo_mean+datam$coast+datam$prec_max +datam$temp_min+datam$evi_min, datam) error in cca.default(d$X, d$Y, d$Z) : All row sums must be >0 i...
2010 Aug 17
2
AIC in MuMIn
...se 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 temp_ran using the code below but failed. Please kindly suggest the potential failure cause. code rm(list=ls()) library(MuMIn) datam <-read.csv("c:/migration/Mig_ratio_20100817.csv",header=T, row.names=1) mig.stds <-glm(datam.std$SummerM_ratio~ datam.std$temp_max++datam.std$evi_mean+datam.std$topo_var+datam.std$topo_mean+datam.std$coast+datam.std$Iso_index_0808,data=datam.std,family=gaussian) target.model &...
2005 Feb 20
2
matrix operations
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 iterating over the entire matrix? Also, does R ha...
2009 Aug 24
1
Saving heatmaps as PDFs
...sing the software "Grab" (Heatmap_CAFvsTNF_run2.tiff). -----R code-------- sample_output <- "stroma_run2" filename <- "exp_limma_results_final_probesets_Pbelowpoint05.txt" # read in data my.data <- read.delim(filename, sep="\t") my.dataM <-my.data[,2:ncol(my.data)] hr <- hclust(as.dist(1-cor(t(my.dataM), method="pearson")), method="complete") hc <- hclust(as.dist(1-cor(my.dataM, method="spearman")), method="complete") # generate heat map library(gplots) if (is.n...
2000 Jan 06
1
nlme
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...
2012 Dec 16
3
averaging X of specific Y (latitude)
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
2017 Dec 14
3
Errors in reading in txt files
...contain 3 elements; only one. > You'll have to fix that line. Delete it, prepend it with a comment character of add enough columns. I definitely don't recommend that. Instead, read ?read.table to learn about the "fill" and "header" arguments. df = read.table("DATAM", header = TRUE, fill = TRUE) will probably work. Best, Ista > > > Berend > >> 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...
2008 Dec 20
1
Help in Lattice!
....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 I use xyplot(Value ~ Visit | factor(ID), data = dataM, typ = 'o',layout = c(3,3),as.table=TRUE); to plot each patient longitudinal observations. How can I add dose information to each panel? e.g., adding "Dose = 30" Best, Haoda
2011 Jul 31
3
export/import matrix
...)) > z <- 10* (1:ncol(y)) > ylim <- range(y) > ylen <-ylim[2] - ylim[1] + 1 > colorlut <- terrain.colors(ylen) > col <- colorlut[y-ylim[1] + 1] > rgl.open() > rgl.surface(x,z,y, color=col, back="lines") Then I export it as write.table(data, file="datam.txt", row.names=TRUE, col.names=TRUE), when I import it back into R again with read.table("datam.txt") it looks like this in R: V1 V2 V3 V4 V5 V6 V7 V8 V9 V10 V11 V12 V13 V14 V15 V16 V17 V18 V19 1 100 100 101 101 101 101 101 100 100 100 101 101 102 102 102 102 103 104 1...
2009 Jan 04
1
POSIXct and chron issues with tz
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 mph2 w/m2 0506010000 0.00 13.6 218.1 36.8 -999 65.1 -999 -999 18.2 0.2 0506010005 0.00 12.9 214.3 36.8 -999 65.5 -999...
2017 Dec 15
2
Errors in reading in txt files
...t 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, > 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 forma...
2009 Jul 10
1
problems with contrast matrix
...LA). All models tested converge but the by-default contrast matrix is really hard to interpret because all the contrasts always involve all the wavelengths. Here is an example of the syntax for lmer and lme, for the complete model: model1<-lmer(PE~Sex+LightIntensity+ordered(WA)+(1|Animal),data=datamic,na.action=na.omit,family="poisson",REML=FALSE) model2<-lme(LA~Sex+LightIntensity+ordered(WA),random=~1|Animal,data=datamic,na.action=na.omit,method="ML") So I built a specific contrast matrix: cmat.r<-matrix(c(-3/6,-2/6,-1/6,+1/6,+2/6,+3/6, # linear effect...
2017 Dec 14
0
Errors in reading in txt files
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 t...
2017 Dec 14
3
Errors in reading in txt files
...2: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 just use the simple script to open it: df = read.table('DATAM', head=T). But it has the error and thus cannot read the file: Error in scan(file = file, what = what, sep = sep, quote = quote, dec = dec, : line 1 did not have 3 elements How to read it with three fixed columns, and how to read the time format in the first column correctly? Thanks for yo...
2017 Dec 14
0
Errors in reading in txt files
...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 just use the simple script to open it: df = read.table('DATAM', head=T). > > But it has the error and thus cannot read the file: > Error in scan(file = file, what = what, sep = sep, quote = quote, dec = > dec, : > line 1 did not have 3 elements > > How to read it with three fixed columns, and how to read the time format in > th...
2007 Jul 19
2
Trend lines on a scatterplot matrix
I'm using pairs() to generate a scatterplot matrix; pairs(~ Fuzzy.gray.white.ratio+Fuzzy.gw.t.score+AgeWhenTested+signal_mean.noise, data=datam,subset=status=="control",main="Controls", labels=c("G/W","Peak Separation","Age","S/N")) How can I add regression lines to the plots?
2006 Apr 08
0
Jupiter Direct Subscription Confirmation Needed
...Update (HTML) DevX Enterprise Update (HTML) DevX Weekly Update (HTML) DevX Java Update (HTML) DevX Open Source Dev Update (HTML) DevX VB Update (HTML) DevX Web Development Update (HTML) DevX Wireless Development Update (HTML) DevX XML Update (HTML) Graphics.com (HTML) Design Firm Management (HTML) Datamation Product Watch (HTML) DRM Watch (HTML) eSecurity Planet (HTML) eCRM Guide (HTML) Executive Tech (HTML) Workstation Planet (HTML) Grid Computing Planet (HTML) Inside ID Weekly (HTML) Instant Messaging Planet (HTML) IT Service Management (HTML) Enterprise Storage Forum (HTML) Optically Networked...
2004 Jul 06
1
R & DataMining
Dear R-user, I wish to know if someone is using R as concern Datamining or KDD (Knowledge Discovery in Database) and if already exists a R package specialized in this kind of analysis. I found this contributes on the R web site: [20] Diego Kuonen. Introduction au data mining avec R : vers la reconqu??te du `knowledge discovery in databases' par les statisti...
2009 Aug 07
0
Bar plots with stacked columns marked with askterisks
...n Data set B i.e. asterisks at the top of the Index1 and Index3 columns. ---- R code ---- title = 'frequency' filename = 'testdata2.txt' sample_output = 'sequences' # read in data my.data <- read.delim(filename, sep="\t") # data from second column onward my.dataM <- my.data[,2:ncol(my.data)] # rownames in first column rownames <- my.data[,1] # define colours for the categories col = c ("red ","orange ","yellow ","green ","violet ","cyan","blue","pink","peachp...