search for: hgt

This is no doubt trivial but after searching the help files and the web, I cannot seem to find it. 1) How do I convert 'hgt' into 'HGT' in R? 2) How should I have used the help facilities to find this? At the end of the day, all I want to do is case insensitive string matching... i.e. 'if ("HGT" == 'hgt') print('this should be true')' I tried using some of the regular exp...

...0.134398711,0.167423957) x <- c(0,1,2,3,4,5) y <- c(50, 100, 150, 200, 250, 300, 350) z <- matrix(z, nrow=length(x), ncol=length(y), byrow=TRUE) #persp(x, y, z, theta = 30, phi = 30, expand = 0.5, # box= TRUE, axes= TRUE, ticktype = "detailed", main="Title of plot") hgt <- (z - min(z))/ (max(z) - min(z)) z hgt cols <- grey(hgt) persp(x, y, z, col = cols, theta = 30, phi = 30, expand = 0.5, box= TRUE, axes= TRUE, ticktype = "detailed", main="Title of plot") Thanks, Sander. > version _ platform i386-pc-mingw32 arch...

Dear All, I have been trying to to do this for a few days now. I can generate persp plot OK and can generate it in colour using the command persp(x,y,z,col=terrain.colors(20)) However, I was wondering if it is possible to shade the 3d surface like a contour plot. i.e. black for large z, white for small z, say Thanks in advance Colin

...ng to their z-value. The colors seems to be transposed, i.e. the color of a facet A seems to be used for a facet B which is a (defined?) way off.... An example: nx <- 21 ny <- 21 x <- seq(-1, 1, length = nx) y <- seq(-1, 1, length = ny) z <- outer(x, y, function(x,y) x^2 + y^2) hgt <- z / max(z) persp(x, y, z, col = gray(1 - hgt), theta = 70, phi = 40, expand = 0.4) Now the z-value should define the color, but the color pattern seems to be moved. Thanks, Oliver Niggemann -.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.-.- r-devel mail...

...with it, and if you come up with something, PLEASE let me know- In the meantime, I've got to start fiddling with excel and figuring out how to automate the solver calculation..... Briefly, the point of the program is to approximate the model output from an iterative calculation, Wtmod and Hgtmod, to user-specified endpoints Wt and Hgt, by seeking the optimal values of p, ACT involved in the iterative process. Also, if your interested in recent correspondence that explains the point of the program a bit, and how the function ties in to the iterative process, search the R help forum...

...ll. Since I have replicate observations on the same individuals I was thinking to use the geepack library. The data are organised in a dataframe with the following variables Date = date of sampling, Size = dimensions (mm) Activity duration of activity (min) Water = duration of splashing by waves Hgt = resting eight of each specimen before activity begin Individual = a code indicating the id of the specimen. I have up to 12 replicate observations for individual. Some observation are missing and I organized the data frame to have exactly 12 rows for each specimen, with NAs where there is a mi...

...level3 var4 2501 62 2 level3 var3 2507 58 2 level3 var2 2508 56 2 level3 var1 2513 63 3 level3 var3 3601 73 3 level3 var2 3603 59 3 level3 var1 3609 56 3 level3 var4 3612 61 modela<-lme(height~variety*fertilizer, random=~1|replication) summary(modela) anova(modela) library(multcomp) hgt <- glht(modela,linfct=mcp(fertilizer="Tukey")) summary(hgt) Any body can help me to proceed tukey (or lsd) with lme that would be highly appreciated. Prabhath University of Saskatchewan [[alternative HTML version deleted]]

...r <- 1:1 } r } .midDend <- function (x) if (is.null(mp <- attr(x, "midpoint"))) 0 else mp ############################## 80 Characters Wide ############################## ########## Body follows type <- match.arg(type) leaflab <- match.arg(leaflab) hgt <- attr(x, "height") if (edge.root && is.logical(edge.root)) edge.root <- 0.0625 * if (is.leaf(x)) 1 else hgt mem.x <- .memberDend(x) yTop <- hgt + edge.root if (center) { x1 <- 0.5 x2 <- mem.x + 0.5...

...4/4/2003 17:01:33 2 3 4/6/2003 16:03:07 5 6 4/8/2003 15:03:08 3 7 4/10/2003 14:03:06 4 5 4/2/2003 13:02:00 4 5 4/4/2003 12:14:43 4 3 4/6/2003 11:00:56 3 5 4/8/2003 10:02:06 2 4 4/10/2003 9:02:19 Weather Data DateTime WndSp WndDir Hgt 4/2/2003 17:41:00 8.17 102.86 3462.43 4/2/2003 20:00:00 6.70 106.00 17661.00 4/2/2003 10:41:00 6.18 106.00 22000.00 4/2/2003 11:41:00 5.78 106.00 22000.00 4/2/2003 12:41:00 5.48 104.00 22000.00 4/4/2003 17:53:00 7.96 104.29 6541.00 4/4/2003 20:53:00 6.60 106...

...########################## > > #USER INPUT BELOW > > #Weight at time 0 > Wo<- 9.2 > > #Hg concentration at time 0 (ugHg/g wet weight) > Hgo<- 0.08 > > #Weight at time t > Wt<- 32.2 > > #Hg concentration at time t (ugHg/g wet weight) > Hgt<- 0.110 > > #Prey methylmercury concentration (as constant) > Hgp<- 0.033 > > #Prey caloric value (as constant) > Pc<- 800 > > #Energy density of fish (as constant, calories) > Ef <- 1000 > > #Maturity status, 0=immature, 1=mature > Mat&lt...

...the package to the latest versione (don't ask me what was the previous R version, cause I don't remember). However, now i am in serious trouble. When i try opening the netcdf file, i use a<-retrieve.nc(filename="/Users/P4blo/Desktop/DATA/Raw/Z500/Z500-1979-2010.nc",v.nam="hgt",t.rng=c(1,31)) just as was doing in the past. but know, i obtain not more a 3D matrix but a 4D one! [1] "ordinary" [1] "Attribute time_origin not found" [1] "The chronology is not straight forward: dt= 24 interval span= 275376 data points= 11475" [1] "ncid$...

Ross posted this code: persp(x, y, z, col = cm.colors(10)[floor(9 * hgt + 1)], theta = 35) I hadn't seen the cm.colors function before; it's very nice! I wanted to modify it a bit to allow the start and end colors as well as the max saturation to be specified, and I think I may have turned up a bug in the hsv() function. Here's the evidence of the bug:...

...**BG6oHO4_okvzJRR!72RvWOqz#;2lXvVWA}MijARbU7ucOhs!XJIMq+ zMN|lrR_|yjYp_y9j)JNb{UvKd$oXcOJ#i~H2lMb&YL?Zax=6|1O4IZj)P&)q(=uCm z!p2xD#QKjAv{bS7DzT+XXoO<Dqr^rip>2vaQHgC+Lc<hmj1n8BgyukrRAO@!dw{G{ z9M_OFmD+W^sOP)I=kN~s3e{;Z;rDr83B0ec+!k~Wg*I7<D6z3RI8zDCR9I-Ll88;E zhS?zaX<=_6T&ajHgtT>VrP8&Crksr;6Dq<=Bw0yGHhFBUCN?&E7^{J?nlQF_7%K?l zV-I5mFjf#om4`8(FcKcdd|=EcjHE(l($gfd5d^TOl6eYE3kiD~+2MSUGr(HqeH`J| zc<<vdx5haiVw256?_&YCK31GVgDH=n;l$5YkDuYld^qt_?P0t{7&RWoTflgWFlrSI z7&HPq7Bwzv*XC`hNG<6|{^4y3M@&;ErY*(rB39B-L)<4KE?>m44sGU1MR!v&irVgF z...