node = hou.pwd()
geo = node.geometry()

# specify point group pattern
ptgrp = "0 1"

# create selection object from ptgrp
sel = hou.Selection(geo, hou.geometryType.Points, ptgrp)

# convert to edge selection
sel.convert(geo, hou.geometryType.Edges)

# get edges from selection
edges = sel.edges(geo)

# keep only contained edges
cedges = []
pts = geo.globPoints(ptgrp)
for edge in edges:
    if set(edge.points()).issubset(pts):
        cedges.append(edge)
        
# create groups if needed    
grp = geo.createEdgeGroup("neighborEdges")
grp.add(edges)
grp = geo.createEdgeGroup("containedEdges")
grp.add(cedges)