Mike O'Brien mhpob

mgcv

gam

Error in gam.reparam(UrS, sp, grderiv) : NA/NaN/Inf in foreign function call (arg

You're likely trying to take the log of zero.

bam(~s(, bs = 'so'))

The following was adapted from the blog post here.

The key is to use results = "asis". More information can be found here. Basically:

anything output with cat() will be evaluated as RMarkdown.
anything output with print() will be evaluated with R.

for(i in seq(min(station_key$year), max(station_key$year), 1)){

VSCode Remote: Connect to WSL2 from another machine

Do you want to do remote development on your WSL2 container in Visual Studio Code? Read this.

Proxy jump method

On the host set up OpenSSH for Windows
Run wsl --update to make sure you are running the latest WSL
Open WSL and install another SSH server inside WSL with sudo apt-get install openssh-server
Now run sudo systemctl enable --now ssh to automatically start ssh when WSL starts.

These are instructions on how to set up SSH into a WSL system on a remote server. They loosely follow the instructions here https://superuser.com/questions/1622581/ssh-into-wsl-from-another-machine-on-the-network. Assume the Windows Host's name is my.windowshost and your user name is me (so [email protected]) and that your user name on the Windows host's WSL is wsl_me. It does not seem to matter what the WSL host name is: we'll be using localhost for the WSL as we're forwarding traffic through the Windows host (my.windowshost in our example).

Set up OpenSSH on the Windows host (just Google how to do this)

	# Create track
	set.seed(8675309)
	turning_angles <- runif(100, 0, 2*pi)
	runs <- runif(100, 0, 50)


	fake_track <- data.frame(x = 0,
	y = 0,
	t = 1)

	[xml]$Content = Invoke-WebRequest -Uri "https://water.weather.gov/ahps2/hydrograph_to_xml.php?gage=rmdv2&output=xml"

	$smry = $content.site.forecast.datum.primary."#text" \| measure -Minimum

	if($smry.minimum -lt 5.5){
	New-BurntToastNotification -Text "Get ready to go sampling!" ,
	"Richmond-Westham is forecast to be <5.5 ft."
	}

	## This was originally used in review, where Year was nested within Blocks

	library(mgcv)

	# In the BodyWeight data, Rat is nested in Diet with no repeated Rats across Diets
	# Analogous to Year nested in Blocks with no repeated Years across Blocks
	data('BodyWeight')

	# Dummy code levels of Diet as for YearBlock in manuscript
	BodyWeight$Diet1 <- ifelse(BodyWeight$Diet == 1, 0, 1)

	library(dplyr)

	df_north <- tibble(
	week = 1:53,
	N = sample(0:15, 53, replace = T)
	)


	df_south <- tibble(
	week = 1:53,

	git branch -d branch-to-delete
	git fetch -p