search for: dbh

I ran into a problem with "subset" while working at home that I am not sure is a bug or not. I defined a custom function to take a data frame of tree positions, sizes, and types (containing the columns TRT, COMP, PLOT, X, Y, DBH and CON) and wanted to make a stem map of the plot using different symbols and colors to represent the different types of tree stems. I copied a bit of the code below, the part that gave me troubles: if(dim(subset(mat,CON==1 & DBH>11.43))[1]>0) symbols(mat[mat$CON==1 & mat$DBH&g...

Hey, So I have a scatter plot and I am trying to plot a curve to fit the data based on a Holling Type III functional response. My function is this: nll2<-function(a,b) { conefun<-(a*DBH^2)/(b^2+DBH^2) nlls2<-dnbinom(x=cones ,size=DBH, mu=conefun,log=TRUE) -sum(nlls) } and my plot is this: plot (DBH,cones) DBH is on the x-axis and cones is on the y-axis. How do I get the curve for my function onto the scatterplot? -- View this message in context: http://r.789695.n4.na...

...o plot an equation in two variables to get a feel for sensitivity to its parameters. I've run expand.grid to get made-up vectors of the combinations of the two independent variables, and am trying to plot the output of the dependent, M, against both the dependent in a 3d space. t <- 0:100 DBH <- 10:100 TxDBH<-expand.grid(t,DBH) time<-TxDBH[,1] dbh<-TxDBH[,2] M <- 4000*(1-exp(-time*(1.104-(0.67*0.7)-0.163*log(dbh))^2)) I've tried wireframe(M~t|DBH, data=TxDBH) wireframe(M~time|dbh, data=M) wireframe(M~time|dbh) also the above model specifications using cloud....

Now I have the forest plot data with x, y locations, and I measured the DBH for every indivicuals. The data looks like that: x=runif(100,0,100) y=runif(100,0,100) dbh=runif(100,1,100) rdata=data.frame(x,y,dbh) > rdata[1:5,] x y dbh 1 99.5354145 1.412844 34.10112 2 0.8259361 87.737036 39.12914 3 6.5678613 65.699032 22.55990 4 67.2987881 72.053...

I've been trying to figure out how to prevent a column that is the result of an aggregate function call so that I can use it in further calculations. For example, I would like to aggregate the expf for the data.frame by sp (character) and dbh (double d=rounded to integer) using the command: > st2 <- aggregate( ntrs$expf, by=list(sp=ntrs$sp,dbh=ntrs$dbh), sum ) > st2$expf <- st2$x / 20 > st2$basal.area <- st2$dbh^2 * st2$expf Warning message: "^" not meaningful for factors in: Ops.factor(st2$dbh, 2) > at...

...s', so nls is definitely in the wrong, but the self-start is almost identical to the one for the logarithmic fit, which works perfectly .... I can't see my mistake at all. Can anyone help?? Thanks, Toby Marthews traceson=FALSE;maxint=10000;minstepfactor=2^(-30);tolerance=1e-7 xtxt="DBH in mm";ytxt="Height in m" dbh=c(0.9,1.0,1.1,1.2,4.8,4.9,5.0,5.1,8.9,9.0,9.1,9.2,11.8,11.9,12.0,12.1,14.9,15.1,15.2,15.5) height=c(5.770089,5.154794,4.888847,6.356770,14.849109,14.973146,15.423816,14.865038,21.335985,20.477624,20.915823,21.886004,23.714724,24.182210,23.954929,23.78465...

...ubtitles. I have 9 charts per page (for combination of 3 blocks and 3 treatments). I want to have subtitles for this interaction (e.g. Block A Trt 1, Block A Trt 2, ...) MBE$bt<- interaction(MBE$Block,MBE$trt) par(mfrow=c(3,3)) for(i in unique(MBE$bt)){ ss <- MBE$bt==i plot(MBE$Year[ss], MBE$DBH[ss]) sm<-loess(DBH~Year, data=pMBE[ss,]) x=seq(2004, 2010, by=1) points(x, predict(sm, data.frame("Year"=x)), type="l") } It was possible for the command: MBE$trt.name<- factor(MBE$trt, label=c("Treatment 1", "Treatment 2", "Tr...

...ounce-request at centos.org You can reach the person managing the list at centos-announce-owner at centos.org When replying, please edit your Subject line so it is more specific than "Re: Contents of CentOS-announce digest..." Today's Topics: 1. CEBA:2005-1224-1 CentOS 4 i386 dbh - bugfix update (Extras Only) (Johnny Hughes) 2. CEBA:2005-1224-1 CentOS 4 x86_64 dbh - bugfix update (Extras Only) (Johnny Hughes) ---------------------------------------------------------------------- Message: 1 Date: Sat, 24 Dec 2005 09:54:03 -0600 From: Johnny Hughes <johnn...

Hello, I have a big table with 3 columns and 103918 rows. This is the example,            time             species                 dbh 5 1 4.9377297 575 1 11.64127213 575 1 109.8182438 575 1 8.029809521 5 1 24.32501874 575 1 4.895992119 575 1 11.40567637 575 1 2.795090562 575 1 21.79281837 575 1 52.57476174 575 1 27.7290919 575 1 3.23262083 575 2 19.30612651 575 1 2.956672964 575 1 111.690689 575 1 1...

Dear people of the help list I am drying to analyze my data using the 'lmer' function and I keep having problems. This is the model: > fm1<-lmer(dbh~spec+scheme+(1|Plot),data=d, REML=FALSE). I analyse tree size (dbh) of 3 different species (spec) and 3 planting schemes (scheme). I have 5 plots, which I hope to model as a random factor. (However, the subsequent output is based on some simplified dummy data, which is based on only two plots an...

R gurus, I'm working on data analysis for a small project. My response variable is total vines per tree (median = 0, mean = 1.65, min = 0, max = 24). My predictors are two categorical variables (four sites and four species) and one continuous (tree diameter at breast height (DBH)). The main question I'm attempting to answer is whether or not the species identity of a tree has any effects on the number of vines clinging to the trunk. Given that the response variable is count data, I decided to use Poisson regression, even though I'm not as familiar with it...

...ironment and the script stops after; $AGI->record_file($fname,'gsm',9,-1);... Since the test environment was different to the production environment; I upgraded my test environment from Asterisk V1 and RedHat 9 to V1.2.4 and CentOS. Now in the test environment my script stops after; $dbh ->do("Update CDR SET Recording='0',GSM_File_Size='$size' WHERE GSM_File = '$finame'"); In both instances there are no perl errors and no errors reported in the CLI. In fact the CLI comes back and tells me everything is OK. Is this a known issue? Has someone e...

Dear all, As part of a larger function, I am randomly removing rows from a data frame. The number of removed rows is determmined by a Poisson distribution with a low mean. Sometimes, the random number is 0, and that's when the problem starts: My data frame: > temp occ x y dbh age 801 0 2977.196 3090.225 6 36.0 802 0 2951.892 3083.769 8 40.6 803 0 2919.111 3075.557 8 40.6 804 0 2914.123 3072.700 9 42.9 805 0 2925.353 3074.675 8 40.6 How many rows (nft) shall be removed? > nft<- rpois(1, 2) > nft [1] 2 Ok remove 2 rows: > temp2<- te...

CentOS Errata and Bugfix Advisory 2005:1224-1 CentOS 4 i386 dbh - Bugfix Update (Extras Only) This CEBA is for the Extras repository only, and not for the main CentOS repository. ------------------- The following upstream changes have been rolled into dbh: * Thu Dec 22 2005 Kevin Fenzi <kevin@tummy.com> - 1:1.0.24-3.fc4 - Bump release to fix tagging is...

CentOS Errata and Bugfix Advisory 2005:1224-1 CentOS 4 x86_64 dbh - Bugfix Update (Extras Only) This CEBA is for the Extras repository only, and not for the main CentOS repository. ------------------- The following upstream changes have been rolled into dbh: * Thu Dec 22 2005 Kevin Fenzi <kevin@tummy.com> - 1:1.0.24-3.fc4 - Bump release to fix tagging is...

...lues. - Here is the function rb(this function should return the type I defined in the database): #hosttype.rb require "mysql" module Puppet::Parser::Functions newfunction(:hosttype, :type => :rvalue, :doc => "Gets host type from db.") do |args| dbh = Mysql.real_connect("mysql_server_ip", "username", "password", "puppet") res = dbh.query("SELECT host_type FROM host_extra_info where name=''" + args[0] + "''") # return res.fetch_row puts res.f...

...er? 2. Could someone tell me that how to connect MySQL with Ruby? I have XAMPP server installed and want to connect with ruby. I am using Ruby1.8.7. I tried following code: require ''rubygems'' require ''mysql'' begin # connect to the MySQL server dbh = Mysql.connect("localhost", "root", "root", "status_publish") puts dbh # get server version string and display it puts "Server version: " + dbh.get_server_info rescue Mysql::Error => e puts "Error code: #{e.errno}&quo...

...="intercepts", xlab="average height per site [cm]", xlim=c(20,3020), ylim=c(-2,5), col=coefficient2$color) highlight(coefficient2$intercept ~ coefficient2$average_height, lbls=coefficient2$site_no,col="Red", cex = .6) plot(coefficient2$intercept ~ coefficient2$average_dbh, main="intercepts ::: dbh", ylab="intercepts", xlab="average dbh per site [mm]", xlim=c(-10,360), ylim=c(-2,5),col=coefficient2$color ) highlight(coefficient2$intercept ~ coefficient2$average_dbh, lbls=coefficient2$site_no,col="Red", cex = .6 If I create a...

...re useless." "My guess is that the deviances and coefficients are entirely ok. I'd expect that problems in the general area that Thomas mentions to reveal themselves as a failure to converge." I have this problem with my data. In a GLM, I have 269 zeroes and only 1 one: summary(dbh) Coefficients: Estimate Std. Error z value Pr(>|z|) (Intercept) 0.1659 3.8781 0.043 0.966 dbh -0.5872 0.5320 -1.104 0.270 > drop1(dbh, test = "Chisq") Single term deletions Model: MPext ~ dbh Df Deviance AIC LRT Pr(Chi) <n...

...when a host is removed from the Puppet network forcibly # Alternatives: run puppet on host to be removed with a blank manifest (node blah {}) # use DBI; my $dsn = ''DBI:mysql:database_server_here''; my $db_user = ''root''; my $db_pass = ''xxx''; my $dbh = DBI->connect($dsn, $db_user, $db_pass); my $query; if (($#ARGV + 1) ne 1) { print "Fail.\n"; print "Usage: $ARGV[0] host_name\n"; } else { my $result; print "Delete from hosts table. Affected rows: "; $result = $dbh->do(''delete from ho...